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CHAPTER  1 


A  SYNTHESIS  TECHNIQUE  FOR  UNCERTAIN  MUI TIPLE  INPfT-iXTPUT 
FEED3ACK  SYSTEMS 


1.1  Introduction. 

The  flight  control  problem  is  basically  a  multiple  input-output 
(mio)  feedback  problem  with  an  acertain  plant.  In  this  report,  the 
lateral  and  longitudinal  systems  are  treated  separarelv.  The  lateral 
system  is  taken  as  three  by  three  with  rudder,  ailerons  and  canard  as 
the  plant  control  input  variables.  The  output  variables  vary  with  the 
mode,  for  example  sideslip  (•!)  ,  roll  ($)  and  side  force  (A^)  in  one 
case.  The  longitudinal  system  has  elevator  and  flap  as  the  plant  inputs 
and  the  outputs  are  usually  normal  acceleration  (A^)  and  pitch  (■?)  . 
In  this  work  the  linearized  plant  equations  are  used,  so  different 
flight  conditions  give  different  coefficients  in  the  linear  equations, 
thereby  giving  a  set  of  plants  P  •  {P}  .  The  fact  that  P  is  only 
known  to  be  a  member  of  a  set  P  (but  it  is  not  known  which  member), 
constitutes  uncertainty.  Of  course  much  of  this  uncertainty  could  be 
removed  by  "scheduling",  but  in  any  realistic  control  system,  there 
always  remains  some  uncertainty.  Nevertheless  one  is  interested  in  a 
feedback  design  which  does  not  use  scheduling,  even  if  only  as  an 
option.  Also,  it  is  desirable  to  know  the  trade-off  between  the 
extent  of  scheduling  and  the  resulting  economy  achieved  in  other  system 
factors.  Hence,  it  is  important  to  have  a  synthesis  technique  for 
uncertain  mio  systems,  which  is  flexible  enough  to  include  scheduling 
of  any  desired  extent. 
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The  difficulty  inherent  in  developing  a  synthesis  theory  for 
uncertain  mio  systems,  can  be  appreciated  by  considering  the  3  *  3 
output-feedback  structure  in  Fig.  1.1,  with  3  *  3  matrices: 

P  =  [p^ji  of  plant;  G  *  [g^ ]  ,  F  =  [f^]  of  compensation  functions. 
We  even  let  G  *  [g^]  be  diagonal  and  consider  only  the  transfer 
function  relating  output  (1)  to  input  (1): 

T11  “(pll  fll  81  +  P12  f21  g2  +  p13  f 31  g3^(1  +  p22  s2)(1  +  p3383) 

-  P23  P32  g3^  "  (p21  fll  81  +  p22  f21  82  +  p23  f31  g3) 

+  E P 1 2  §2^  ^  33  S3)  —  P ^2  ^3_3  ^2  §2^ 


^P31fll8l  *  P32f2182  +  p33f3183^P23p12g283  ~  P13g3^1+P22g2)  ^ 

U  +  Pu  gjMU  +  P22g2)(1  +  P3383>  ~  P23P328283] 

-  Pi2glfpl282(1  +  P33g3)  ~  P32P13g2g3^  +  P13gl ^P12P238283  "  p1383(1+k282 

(1.1a) 

There  are  nine  such  functions,  each  one  a9  complicated  as  this. 

Note  that  the  variables  normally  refer  to  Laplace  transforms  or  transfer 

functions.  Functions  of  time  will  be  explicitly  expressed,  e.g.  f(t)  . 

Suppose  that  there  is  significant  uncertainty  in  the  functions 

2 

p^  ,  for  example  p^  -  (A^s  +  Bij^Eijs  +  Jijs  +  Fij^  ’  with 
intervals  of  uncertainty,  for  (ij)  *(11)  say  Au  6  [1,2]  , 

Bu  6  [2,4]  ,  En  €  [1,2.5]  ,  Jn  £  [-1,3]  ,  Fn  €  [-2,1]  ,  or 
each  may  be  a  function  of  other  parameters  with  large  uncertain 
intervals  (as  in  the  flight  control  problem).  The  problem  seems 
hopeless,  unless  one  is  thoroughly  familiar  with  the  system,  with 
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considerable  cut  and  try  experience  so  that  he  has  a  good  initial 
starting  point. 

A  synthesis  technique  has  been  developed  fl-4  1  which  considerably 

reduces  the  complexity  of  the  above  problem.  Since  it  has  beer,  presented 

i 

;in  detail  in  the  references,  it  is  only  summarized  in  the  next 
sections. 

jl. 2  Summary  of  MIO  synthesis  technique 

Fig.  1.2  applies  for  a  3*3  system  and  consists  of  nine  single- 

|loop  structures.  Note  that  in  the  first  row  the  elements  inside  the 

|loops  are  G^(s)  ,  Q^(s)  the  same  for  each  of  the  three  structures. 

But  the  disturbances  d^(s)  are  different.  Also,  the  prefilters 

are  different.  This  situation  is  repeated  for  each  of  the  three  rows. 

The  F  are  the  elements  of  the  F  matrix  in  Fig.  1.1,  while  the 
uv 

are  the  diagonal  elements  of  the  G  matrix,  which  for  the  moment 
is  taken  as  a  diagonal  matrix  (see  Sec.  1.4).  The  Quv  in  Fig.  1. 
are  the  inverses  of  the  elements  of  the  P  ^  matrix,  P  1  *  l^^uv^ 
where  P  is  the  plant  matrix.  Because  of  the  plant  parameter 
uncertainty,  there  is  a  set  P  *  {p}  of  plant  matrices,  which 
generates  a  set  of  Q.  matrices.  Since  y  *  PS  in  Fig.  1.1, 

<5  -  P-1y  ,  so  -  yL/Qu+  ’  6  2  *  yl/Q21  +  y2/Q22  +  y3/Q2! 

etc.  The  following  is  a  physical  interpretation  of  (for  example)  Qi2  : 

Adjust  the  inputs  6^  ,  <$2  ,  so  that  y^  ■  y^  *  0  .  Then  the 

.resulting  ratio  of  transforms  y2/6^(s)  is  equal  to  Q^2  ,  y2/62  “  Q 22  ’ 

jy^/f^  *  Q^2  •  In  practice,  one  simply  inverts  P  in  order  to  obtain 

the  Q 
^uv 


The  elements  b  in  Fig.  1.2  are  real  positive  functions  of 
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frequency  to  ,  obtained  as  follows.  The  3  >•  3  system  in  Fig.  1.1  has 
nine  closed-loop  transfer  functions 

y  (s) 

T  (s)  =  -}L~—  •  : 

ov  r  (s) 

v 

In  many  problems,  for  u  ^  v  one  wants  to  be  small,  ideally  zero, 

which  is  certainly  possible  at  one  combination  of  plant  parameter  values. 
But  if  there  is  uncertainty,  this  cannot  in  general  be  achieved  for  all 
P  €  P  ,  so  one  sets  a  lower  bound  which  can  be  a  function  of  w 


T  (jto)  <  b  (ui)  ,  u  ^  v 
uv  uv 


(1.2) 


to  be  achieved  for  all  P  €  P  .  Such  a  system  is  called  "basically 

noninteracting"  BNIC.  In  this  report  the  designs  are  all  desired  to  be 

BNIC.  Usually  one  wants  some  sensible  response  for  the  diagonal  elements 

T^Cs)  in  Fig.  1.1,  so  in  minimum-phase  systems  [3],  specifications  on 

|T  (jui)  |  suffice  of  the  form 
uu 


a  (w)  5  T  ( iw)  $  b  (w) 
uu  uu  uu 


(1.3) 


See  Chapter  2  for  examples  of  b  ,  a  of  (1.2,  1.3).  Thus,  the 

uv  uu 

disturbances  d  in  Fig.  1.2  are  derived  from  the  Q  .  of  the  plant 
uv  ui 

and  the  specifications  b.  on  the  closed-loop  response  functions. 

im 

The  problem  in  Fig.  1.2,  is  to  choose  the  ,  F _  so  as  to 


satisfy 

the  specifications 

on  the 

T  of 

ij 

Fig.  1.1;  for  example 

te 

choose 

F11  ’ 

so 

that 

yll 

satisfies 

the  specification  (1.3) 

$ 

with  u 

-  1  . 

Choose 

F12 

•  G1 

so  that 

y^2  satisfies  (1.2)  with 

*.  *  • 

u  -  1  , 

v  =  2 

,  etc. 

The 

output 

*12  in 

Fig.  1.2,  is  the  sum  of 

two  components:  one  due  to  the  unit  input,  the  second  due  to  d^  • 
As  it  is  required  (Eq.  1.2)  that  -  b12^  ^or  al*  ^  *  ic 


Is  clearly  best  to  let  F^  *  which  we  can  control,  be  zero.  Similarly 

let  F  *  0  for  all  u  4  v  .  Since  G,  is  the  same  for  all  3  structures 
uv  1 

in  the  first  row  of  Fig.  1.2,  it  must  be  "good  enough"  (have  large  enough 

.rain  over  large  enough  frequency  range)  to  '■'.v;. the  "feedback  needs"  of  the 

worst  of  the  three  at  each  a  value.  If  it  is  possible  to  find  F. .  , 

ij 

which  solve  these  nine  single-loop  problems,  then  it  is  guaranteed  that 

these  same  F . .  ,  G.  when  used  as  elements  of  the  F  ,  G  matrices  in 
ij  i 

Fig.  1.1,  are  satisfactory  for  the  original  mio  problem  [1].  Thus,  the 
3  *  3  mio  problem  has  been  converted  into  three  loop  problems  (choice 
of  G^)  and  nine  prefilter  problems  (choice  of  F_). 


1.2.1  Derivation  of  Table  1.1. 

Table  1.1  lists  the  demands  on  the  three  loop  transmission  functions 
G^  so  as  to  solve  the  nine  single- loop  problems  of  Fig.  1.2, 
for  the  BNIC  case.  They  were  obtained  as  follows.  In  Fig.  1.2  each 
diagonal  y  element  has  two  components,  one  due  to  the  unit  input, 
the  second  due  to  d^(s).  Let  their  respective  outputs  be 


y  (s) 
'uu 


t  (s) 
uu 


T  (s)  +  T  (s) 
uu  duu 


(1. 4a-d) 


d  Q 
uu  uu 


duu  1+L 


G  Q 
u  uu 


From  (1.3)  it  is  necessary  that 


v  (jw)  €  [a  (oj)  ,  b  (uj)  ]  . 

•uu  uu  uu 


(1.5) 


3ut  ,  v  (j-i)i  =  i  t  (jui)  +  Ci-j)  ;  ,  and  since  their  relative  phases 

'  uu  uu  duu  r 

are  not  known,  the  harsher  condition  is  used: 


t  ( jto)  +  t  (jcu)  £  [a  , b  ] 
uu  1  duu  '  uu  uu 


(1.6) 


Suppose  |t  ( j  to )  6  [a'  (oj),b'  (to)  ] 

uu  uu  uu 


(1.7) 


which  must  be  chosen  such  that. 


b  +  r  $  b  ,  a  -  t  5a, 
uu  duv  uu  uu  duu  uu 


[b  -  a  +  2t  ]  =  b  -  a 

uu  uu  duu  max  uu  uu 


(1. 8a- c) 


For  example,  suppose  | Yuu^ J 2) |  €  [.9.1.1]  is  required.  This 

allowance  of  1.1-.9  =  .2  variation  in  |T  (j2)|  can  be  split  up  in  an 

infinitude  of  ways  between  t  (j2)  and  t.  (12):  e.g.  t,  =  .05  , 

uu  duu  °  duu 


b  =  1.05,  a  ■  .95  satisfying  (1.8c).  But  tj 
uu  uu  y  °  duu 


.07  ,  b  =  1. 
uu 


a^u  ■  .97  would  also  be  satisfactory,  etc.  What  is  the  optimum  choice? 

To  satisfy  (1.7),  a  minimum  value  of  | L ( j to )  |  is  needed,  and  a  minimum 

value  of  | L(  1u>)  I  is  needed  to  achieve  It.  |  j  its  allowed  value. 

1  duu1 

Hence,  the  optimum  division  is  when  these  two  minimum  required  | L j 
values  are  the  same. 


In  table  1.1,  the  (1.7)  obligation  on  L  is  denoted  as  an  A  t 

u  uu 


of  the  form 


1.2.2  Optimization  and  trade-off 

There  are  interactions  betweeen  the  entries  of  Table  1.1,  via  the 


b^v  .  For  example,  suppose  b2^  is  decreased  from  its  original  value  of 
-3  -3 

say  10  to  6  <  10  ,  which  certainly  doesn't  violate  the 

_3 

requirement  ^  •  This  helps  ,  because  |l  +  |  now 

may  be  less  than  before.  It  helps  for  the  same  reason,  but  clearly 


makes  it  harder  on  .  But  suppose  imposes  a  tougher 

obligation  on  than  does  (this  is  denoted  by  ^  5 

~  D22  means  that  each  imposes  the  same  requirement  on  L^)  • 

Then  it  doesn't  matter  if  b2^  is  decreased  to  some  extent,  thereby 
increasing  the  obligation  on  -  Note  that  this  interaction  is 

always  between  the  elements  of  the  same  column  in  Table  1.1.  Thus  b^  , 
b^  ^  ,  b ^  appear  only  in  the  first  column;  b^  >  ^2  ’  ^32  on-*-y  in 
the  second  column  etc.  Similarly,  decrease  of  b^  helps  (A  ^,B  j.)  , 


but  makes  it  harder  on  . 

In  the  above,  it  was  possible  to  help  at  no  expense  to  the 

other  L's.  Equilibrium  is  defined  as  when  this  is  no  longer  possible. 
The  following  has  been  proven  [2]: 

Th eo rem ;  A  necessary  and  sufficient  condition  for  equilibrium  is 

when  all  the  L  are  dominated  by  the  members  of  the  same  column, 
u 

For  example,  dominates  ,  (A22,  B22^  dominates  L?  , 

D^2  dominates  •  This  does  not  preclude  one  or  more  of  ~  , 

( ^  ,  3^^)  ^21 

Trade-off 

After  equilibrium  has  been  reached,  it  is  still  possible  and  may 

be  desirable  to  trade-off,  i.e.  sacrifice  one  or  more  L  for  others. 

u 

For  example,  if  the  sensor  of  y2  is  noisier  than  the  others,  or 

Q22  has  higher-order  modes  in  a  rather  low  <u  range,  it  is  important  to 

reduce  |L2(jm) |  rapidly  vs.  w  .  This  can  be  done  by  reducing  b^  (°r 

if  ^23  */*  ^A22’  B22^  »  °21  at  Che  exPense  of  °13»  (A33»  B33>  , 

thereby  imposing  greater  demands  on  (if  dominates  L^) 

and/or  on  ,  if  (A^,  B33)  dominates  .  Of  course,  if  so, 

eventually  one  or  more  of  (A^,  B22>  ,  D21  ~  D23  say  the  former. 

Then  b22  can  be  decreased  but  only  until  it  approaches  a22  in  value, 

whereas  any  b  (u  ^  v)  can  be  reduced  by  any  amount. 

Equilibrium  and  trade-off  give  the  designer  valuable  flexibility, 
which  he  can  exploit,  according  to  the  circumstances  of  his  specific 


system. 
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1.3  Necessary  conditions  for  use  of  of  MIO  synthesis  technique 

Certain  conditions  must  be  satisfied  in  order  that  the  above 

synthesis  technique  may  be  used.  These  conditions  are  not  difficult 

to  derive,  by  simply  examining  the  operations  involved  in  executing 

the  design  techniques.  Normally  the  plant  matrix  P  is  available,  so 

its  inverse  P  ^  must  exist,  in  order  to  evaluate  the  Q  from 

uv 

P  ^  =  [1/Quvl  •  This  means  that  det  P  must  not  be  indentically 

zero  for  any  P  €  F  .  This  is  the  much-touted  "controllability"  condition. 

If  dec  P  =  0  it  means  that  the  plant  outputs  y^  are  not  independently 

realizable,  for  that  specific  P  . 

The  other  conditions  are  due  to  the  need  that  the  single-loop 

problems  of  Fig.  1.2  be  solvable,  for  only  then  can  the  F  ,  G  be 

u  u 

found.  The  uncertainty  in  the  Q  must  be  compatible  with  the 

specifications  on  the  T  (jw)  .  If  Q  (s)  is  minimum-phase,  then 

uu  uu 

very  general  uncertainty  and  T  performance  tolerances  can  be 

handled  (see  [1]  for  details).  But  if  Q  (s)  is  nonminimum-phase 

uu 

(has  zeros  in  the  right  half  complex  plane) ,  then  the  T^u  tolerances 
must  be  compatible  with  the  limitations  on  the  realizeable  bandwidth 
of  Lu(jw)  [4],  due  to  the  right  half-plane  zeros  of  Q^ts)  ■ 

Right  half-plane  poles  of  Q^u  can  be  handled  [2 ] . 

1.3.1  A  high-frequency  condition 

There  is  an  important  high-frequency  condition  on  the  Q^v  which 

must  be  satisfied.  It  will  be  derived  for  the  2x2  case,  and  only 

stated  for  the  3  *  3  case.  At  high  frequencies,  the  a  of  Fig. 

uu 

1.2,  Table  1.1  and  Eq.  (1.3)  are  so  small  that  they  can  be  taken  as 

zero.  All  the  A  of  Table  1.1  then  disappear  and  the  B  have  tne 
uu  r  uu 


In  the  2*2 


same  form  as  the  D  ,  with  b  replacing  the  t 
uv  uu 

case  the  inequalities  of  Table  1.1  then  are: 


1  4. 

T  | 

-1  <  LlI 

Q12 

b12 

Q12| 

L1! 

k  |  i 

'QU 

’  b 

QH! 

|-1  e  b21 

Q21  j 

b22 

lQ2ll 

1  + 

L2 

1  '  ba 

Q22 1 

b  12 

lQ22  i 

In  (1.11a),  one  might  as  well  let  bi]/b2i  =  ^12^22  sance  I1  + 

must  satisfy  both  inequalities.  This  can't  hurt  and  may  help,  because 

suppose  b12//b',2  >  bll^b21  '  Reducin8  b^/b^  doesn’c  hurt  buC 

may  help  (b)  ,  if  b^/b^  <  b21^bll  ’  F°r  the  same  reason  set 
b21^bll  *  b22^b 12  in  ^ •  These  Cwo  results  are  identical: 

bUb22  "  b12b21  •  “d  8lve 

n  +  t,rlii  +  L,rl  t  fo-ffi1-)  •  u.u) 

i  i  ^Hy22  1 


Now  we  know  that  sooner  or  later  in  w  ,  any  practical  |l  +  L  |  <  1 

unless  Lu  is  restricted  to  |arg  L|  5  */2  which  is  impossible  in  any 

practical  system,  (This  impractical  condition  is  blithely  assumed  in  the 

"optimal  quadratic  regulator,",  so  popular  in  Modem  Control  Theory). 

If  we  play  safe  and  want  to  allow  |l  +  L  |  <1  for  all  w  larger  than 

some  w  ,  then  (1.12)  requires  that  there  exist  a  to  such  that 
u  um 


1  P11P22 |  _  Q12Q21 
I  P12P21 1  Q11Q22 


>  1  for  all  P  £  P  and  all  to  >  u 


(1.13) 


The  analogous  condition  in  the  3*3  case  is  more  complicated  in  its 
derivation  [1],  so  is  only  stated  here: 


12 


PL1P22P33*  >  IP1LP23P32I  +  'P12P21P33'  +  ^P12P23P3l‘ 


+  PPP  +  PPP  l 
1  132231 1  !  13  21  32 1  ’ 


(1.14a,b) 


P 

uv 


1/Q 


uv 


The  above  are  not  necessary  conditions  because  in  (1.12)  we  could 
stagger  the  w  regions  in  which  j 1  +  |  ,  il  +  Lvl  are  less  than 

unity,  i.e.  have  L  in  the  second  or  third  quadrants  when  L  is  in  the 
first  or  fourth  and  vice  versa.  3ut  it  simplifies  matters  if  the 
condition  is  satisfied.  By  renumbering  the  terminals  at  the  plant  input 
(but  not  at  the  output)  the  new  p^  is  old  p^2  ,  (p^  ew  =  ^pll^old 

etc.  And  if  this  renumbering  is  done  at  the  plant  output. 


(p  .)  ■  (p-,)  i  .  etc.  In  both  case  (1.13)  becomes 

11  new  21  old 


(  Q12q21  \  a  (  qilQ22  \ 

V  QU^22  4iew  ^  Q12q21  4ld 


(1.15) 


Of  course  renumbering  of  either  the  input  or  the  output  terminals  of  the 
plant  changes  the  expressions  (1.14,1.15).  This  may  be  tried  if  in  a 
specific  numbering,  the  inequality  is  not  satisfied  VP  €  P  .  It  is  worth 
noting  that  Rosenbrock  [5]  has  a  condition  for  his  synthesis  technique, 
which  is  harder  to  satisfy  than  (1.13)  or  (1.15),  even  though  his  synthesis 
technique  does  not  cope  with  plant  parameter  uncertainty.  His  condition  is 
that  the  diagonal  element  IpuuI  be  large  than  the  sum  of  all  the  other 
elements  in  the  same  row  (or  of  the  same  column) ,  and  this  must  be  so 
for  all  j  6  [0,”)  . 
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Finally,  suppose  that  in  the  2x2  case  neither  (1.13)  nor  (1.15) 
is  satisfied  for  all  w  >  some  ,  precluding  the  use  of  this  synthesis 
technique  in  precisely  its  present  form.  This  is  the  situation  in  some 
of  the  moc.'S  of  our  aircraft  flight  control  problem,  no  matter  what  the 
numbering  of  the  plant  terminals,  i.e.  (1.14)  is  not  satisfied.  Never¬ 
theless,  all  is  far  from  lost.  In  the  2x2  case  one  uses  this 
synthesis  technique  for  designing  one  of  the  (say  u  =  1) .  So 

fhe  design  of  ,  F^  is  that  of  single-loop  system  design.  Then  one 
uses  the  exact  system  equations  to  design  ,  F,  .  This  second  part  is 
not  difficult  because  it  too  becomes  a  single-loop  type  of  problem,  when 
is  known: 


P11P21F1  L1 


C  -  hs?l 
22  L20+W  ’ 


21 


p110  (1+L1} 
L20+  P 


(1. 16a-d) 


+  L.) 
Pi  1  P'10  1 


11  11  22 


10 


(1  +  Lj) 


P11P22"P12P21 


^20  *  ^2^220  a  nomina^  1°°P  transmission  at  some  €  P  chosen 
as  the  nominal  plant.  If  I^q  is  large  enough  then  from  (1.16a) 

T22  *  F2  ‘  In  large  enough  I^q  also  guarantees  | T2 1  f  can 

be  made  as  small  as  desired  over  any  finite  u>  range.  Since  F^  and 

are  known,  it  is  simple  enough  to  choose  to  satisfy  the  tolerances 

on  T21  ,  T22  •  this  way,  the  restriction  (1.13)  can  avoided,  and  .in 
a  similar  manner  restriction  (1.14)  for  the  3X3  system  (see  Sec.  2.4). 

The  above  situation  in  the  YF16CCV  system  motivated  us  to  research  this 
problem  and  develop  an  improvement  of  the  mio  synthesis  theory.  This  is 
presented  in  Appendix  1.  It  is  a  good  example  of  fundamental  research 
motivated  by  practical  problems. 
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A  shortcoming:  use  of  nondiagonal  G 

A  good  criterion  of  any  feedback  design  technique  is  that  the 

amount  of  feedback  should  be  related  to  the  amount  of  parameter  un--e  rt.i  in  tv  and 

co  narrowness  of  die  tolerances.  Ln  partieuLir,  i:'  there  is  n>>  ur.cc  rt  .  ir.t 

there  is  no  need  for  any  feedback  at  all,  and  the  design  technique  should 

emerge  with  this  result.  As  presented  so  far,  the  above  design  technique 

fails  this  test  because  of  the  existence  of  the  B..  ,  D..  elements  in 

ij 

Table  1.1.  This  can  be  corrected  by  inserting  a  precompensator  H  in 
front  of  the  plant,  whose  function  is  to  diagonalize  the  new  modified 
plant  PH  ,  as  much  as  possible.  If  there  is  no  uncertainty  it  can  be 
done  precisely  and  then  Q„  =  00  for  all  i  4  j  ,  and  in  Table  1.1 
the  all  give  |  1  +  L^|  1  $  ®  with  no  resulting  obligation 

on  the  L 

u 

When  there  is  parameter  uncertainty,  the  elements  of  H  can  be 
chosen  to  minimize  the  B^  ,  factors  in  the  design  problem.  Let 

the  new  modified  Q  elements  be  denoted  by  w„  ,  the  old  (due  to  P  alone) 
by  .  For  the  2  *  2  case 

[1/v^]  -  (PH)"1  *  H"1n/QiJ ]  -  I 


'!u/i  + 

X  — ) 

lii/ 

X  ^ 

QllA 

1  <>21  > 

1  ^22 

V22  /  , 

^21  \ 

V22  / 

^22 

q2lV 

**22  V 

\  ■  ■  ■ 
1  Q12 

(1.17a-d) 


[ v i j 1  =  H~  .  Xl  “  'Jl2/Vll  ’  X2=V21/V; 


In  the  low  and  medium  ut  range  where  |L  t  is  >1  significantly 

i  ^22 

1 1+L  !  saa  | L  i .  so  D,  o  for  example  becomes  ciosely  |  G,  |  5  - - 

u  u  12  i  °12W12 


Hence  it  is  desirable  to  choose  A  to  minimize  max  .  . 

1  P  1  Q11  Q22 

This  is  also  the  best  that  can  be  done  for  (A^,B^)  in  this  u>  range. 
In  the  higher  range  where  i L  I  ~  1  or  less  (l.lla.b)  apnlv  so  then 


one  seeks  to  minimze 


max  | w , , /w 


11'  12 1 


max 


Q11  ^21 


-L  +  -J: 
QL2  Q22 


(1.18) 


if  the  right  side  in  (1.11a)  really  is  the  constraint  on  .  Often 
the  right  side  is  so  large  that  instead  one  arbitrarily  sets  a  gain 


margin  type  of  constraint  like  | 1  +  L, 


-1 


<  G  db  in  the  higher  -o 
m 


range. 


1.5  Summary 

We  have  presented  the  highlights  of  a  multiple  input-output  (mio) 
synthesis  technique  for  uncertain  feedback  system.  Its  main  attractive 
feature  is  its  reduction  of  the  highly  complex  uncertain  mio  system  into 
a  number  of  single  loop  uncertain  systems.  Under  very  general  conditions 
the  solutions  of  the  latter  are  guaranteed  to  be  satisfactory  for  the 
original  mio  problem. 

The  technique  was  applied  very  successfully  to  four  linear  time 
invariant  mio  problems  with  highly  uncertain  plants  [1-  3  ]  .  The  most 
complex  was  3x3  with  a  mixture  of  BNIC  and  more  complex  performance 
tolerances.  Also,  it  was  extended  to  a  class  of  nonlinear  uncertain  mio 
plants,  including  a  2  x  2  highly  nonlinear,  uncertain  and  interacting 
example  [4],  However,  all  of  these  were  fabricated  academic  examples. 
The  balance  of  this  report  applies  this  technique  for  the  first  time, 
to  a  real  life  practical  design  problem. 


CHAPTER  2 


DESIGN  FOR  DIRECT  SIDE  FORCE  MODE 


The  synthesis  technique  of  Chapter  1  is  applied  to  the  design  of 


the  "Direct  Side  Force  (A  )  "  mode  of  a  fighter  CCV  aircraft.  Direct 

y 


side  force  (Ay)  is  commanded  and  the  objective  is  to  achieve  specified 


Ay  step  response  tolerances,  while  maintaing  specified  small  sideslip 


(3)  and  roll  ($)  .  The  control  surfaces  are  <5  (rudder),  <5  (ailerons) 

r  a 


<5^  (i~Tnard).  The  definition  and  units  of  the  symbols  are  all  listed  in 


the  List  of  Symbols  at  the  beginning  of  the  report. 


2.1  Plant  equations 
The  plant  equations  are 


sx”Ax  +  B<^ 
x  “  [p,r,3] ' 


1  -  U,<J>,Ay] 


2.  ■  C  x  +  D  ^ 

i  -  [6  6  ,6  ] ' 

—  r  a  c 
p  ■  i  ,  r  *  i 


yaw  rate 


(2 . la-g) 


qSb  $ 


2V 


2  VC 


16 


'Ip 


'Hr 


b 

2  VC 


n6 


np 


nr 


C  +E 

yp 


2mV 


2  VC 


yr  qSb 


i£ 


(2.2) 


V 

s  V 


i 


p* 

[l' 

y' .  . 

kv- 


£ 


’¥•  '.  »  .■  '.  *  .■'.■•  H  'A'  .'■"  ■>  'V*  ' 
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r 


^ior  ^2.da  ^£.6c 


B  =  q  S  ^  I  Cn,r  Cn,a  V..  i 


(2.1) 


C  .  C  .  C  . 

yOr  yoa  yoc 


i>  = 


b/I> 

0 

0 


0 

b/I 


0 

0 

1/mV 


Y  2 

f»  _  /  i  41  2  mV 

E  =  (a+  7"  > 


*v. 

0  0  1 

o 

o 

L_ 

0 

X*'  c  ■ 

1/s  0  0 

.  D  * 

0  0 

0 

s 

C31  C32  c33 

_  Cy<$r  Cy6a 

Cy6c_ 

(c  c  c  )  =  b<^E  (C  C  —  c  ) 
^  31,C32,C33;  2Vm  yp  ’  yr’  b  yBJ 


Note  that  the  third  row  of  D  is  precisely  the  product  of  V 
and  the  third  row  of  B  . 

From  (2.1a)  x=(sl  -  A)  ^B  5_  , 


(2.4a,b) 


u. 

m 


(2 . 5a,b) 

(2.6a) 


[C(sI-A)-iB+D]<5  =  P  8 


(2.6b) 


Our  technique  requires  evaluation  of  P 


-1 


U/Qyl 


Let 


R  =  DB_1 


0  0  0 

0  0  0 

0  0  V 


(2.6c) 


from  (2.3,  2.5b),  so 


[l/Q..]  =  P"1  =  { [C(s I-A) _1  +  DB~^]B }_1  =  B_1 [ R  +  C(sl  -  A)”1]' 
=  3_1v  [R(sI-A)+C](sI-A)_1rL  =  B_1(sl-A)  [R(sl-A)+C]_1  . 


(2.7) 


It  is  readily  found  that 


Cs+R(s  I-As) 


0  Os 

1  0  0 

-(as+YJV  Vs  Vs‘ 


(2.8a) 


and  its  inverse  is 


(as+Y  )/s 
<9 


(2.8b) 


P  1  *  [1/Q^j]  is  obtained  by  premultiplying  (2.8b)  by  B-1(sI-A)  .  This 
was  programmed  on  the  digital  computer. 


2.2  Flight  Conditions 

Ten  flight  conditions  were  considered  in  pairs  of  two.  Conditions 
(1,6)  are  quite  similar  to  each  other,  so  are  (2,7),  ...  (5,10)  in 
having  the  same  mach,  velocity  etc,  differing  only  in  the  Ci5c  values, 
due  to  the  odd  cases  having  v.c.  =  15,  the  even  cases  v.c.  =*  25  .  The 
parameter  values  for  these  conditions  are  listed  in  Table  2.2.  The 


Quv  transfer  functions  for  the  10  cases  are  listed  in  Table  2.3  in  the 
form  =  k.ir(s+p^)  .  The  value  of  k  is  listed  first,  followed  by  the 

p,  .  At  the  end  in  brackets  is  the  zero  frequency  value  of  Q  ^  .  For 


example  in  Case  1,  =  .0723(s+.907+j  8.2)  (s+.C)07-j  8.2)  , 

whose  value  at  s  =  0  ,  is  4.92  .  *  -.0909(s+6. 13) (s-6. 10)  whose 

value  at  s  =  0  is  3.40  ,  etc. 

The  ten  flight  conditions  give  it  each  frequency  ten  different  values 

of  Q  (jw)  .  The  sets  {Q  (iuj)}  ,  denoted  as  "templates"  Tp  (oi)  , 
uv  uu  uu 

are  needed  in  the  synthesis  technique.  A  number  of  these  are  shown  in 
Figs.  2.1-3,  in  Nichols  chart  coordinates.  A  missing  flight  condition 
number  in  these  figures,  means  that  its  Q^(jo>)  value  is  very  close  to 
that  for  another  flight  condition.  The  star  refers  to  Case  10. 


Table  2.2 


Fighter  CCV ,  Trimmed  Aero  Dace,  01/18/79 

Area.  ?  =  J'^O.OOO  ft'  ;  Mae,  e  =  10.9  ft.  :  .-'pm.  h  =  '9.0  ft. 

I  _  =  12.606  ,  I  =  12.824  ,  Lood  factor  =  1.0, 

a  h  xy 

CG  location  35.300  %  Mac.  at  W.L.  93,200  .  I  =  3100  ,  I  =  47000 

,  XX  yy  ’ 

:  =  53300  ,  I  =354  slug  ft'  . 

zz  xz 


Flight 

1 

2 

3 

4 

5 

Condit ions 

6 

7 

8 

9 

10 

Mach 

.6 

.8 

.9 

.6 

.9 

Altitude  (ft) 

0 

20000 

0 

30000 

30000 

Impact  P,  Qc 

582.98 

510.38 

1462.84 

173.48 

435.30 

Dynamic  P,q 

533.22 

435.97 

1199.75 

158.67 

357.01 

(lb/ft2) 

P  /P 

1.00 

.460 

1.00 

.298 

.291 

s  o 

.69 

895.41 

Q  /P 
c  s 

Vel  (ft/sec) 

.276 

669.84 

.524 

828.32 

.691 

1004.76 

.276 

596.94 

TRIMMED  LATERAL  DIRECTION  DATA 


a  (degrees) 

i.  890 

2.145 

.901 

6.129 

2.470 

103C  (1/deg) 

n6 

2.291 

2.759 

2.169 

3.463 

2.947 

10^(1 /deg) 

-1.836 

-1.999 

-1. 721 

-2.705 

-2.005 

102C  (1/ deg) 

y  p 

-2.030 

-2.201 

-2 . 100 

-2.203 

-2.220 

102c  (1/rad) 

np 

-1.949 

-2.658 

.9888 

-2.709 

-2 .  780 

10(C  ( 1/ rad) 

IP 

-2.908 

-3.148 

-3.180 

-2.790 

-3.203 

10C  (1/  rad) 

yp 

1.054 

1.329 

.5486 

1.495 

1.227 

10(C  (1/deg) 

nr 

-3.834 

-3.665 

-3.562 

-3.863 

-3.837 

_  VI  _ 


Con'c  Table  2. 

2 

M 

10C ^  r( 1/ deg) 

1.  752 

2.109 

1.608 

2.324 

IOC  (1/deg) 
yr 

6.910 

6.  .396 

6.434 

7.218 

:  -  •' 

103C  -  (i/deg) 

-1.217 

-1. loo 

-  .0989 

-1.-65 

a  -  r 

B 

10^Clor(1/de8) 

3.  720 

3.090 

2.057 

3.229 

>V- 

103Cy6r( 1/deg) 

2.347 

2.147 

1.240 

2.873 

>  *  . 

103C  (1/deg) 

a  a 

-1.008 

-1.082 

-1.067 

-.9541 

a 

“^tS, 

-1.012 

-1.060 

-  .9443 

-1.089 

103c  . 

2.343 

2.567 

2.508 

2.352 

<  *  • 

yoa 

3 

10  C  . 

1.310 

1.314 

1.306 

1.190 

no  c 

1.010 

i.ooa 

1.027 

.8761 

l°ict6c 

2.091 

2.253 

3.058 

1.290 

V' 

2.104 

2.199 

3.075 

1.293 

i 

lo3cy6c 

2.066 

2.027 

2.207 

1.396 

» ■  * 

1.800 

1.770 

1.912 

1.248 

2.136 
7.  351 
-  1 .  i  >8 
3.306 
2.080 
-1. 107 
-1.082 
2.671 
1.310 
.9977 
2.816 
2.727 
1.965 
1.720 


Tr.msfrr  functions  of  Q  for  10  flight  conditions 


Con't  Table  2. 


.077;  9. 79,-10.0(8.22)  .0310;  1. 39  0648(-2 . 97F-03)  .860E-04;  425(.0366) 


It  was  found  that  condition  (1.14)  was  not  satisfied  for  all  P  €  P  , 


no  matter  what  plant  terminal  numbering  w.j.s  tried.  The  modification 
described  in  Sec.  1.  5,  was  therefore  used:  The  technique  is  used  in  its 
original  form  to  derive  and  L?  .  Then  the  exact  equation  is  used  to 

derive  •  First,  the  bounds  on  |T^(jw)|  are  needed.  Examination  of 
the  simulation  data  supplied  by  the  contractor  (A^  mode  response  -  Fig.  2.4), 
suggests  modelling  T33(s)  as  a  second-order  function  with  a  pair  of  poles 
at  approximately  -5.  This  gives  a  first  cut  at  iT^Cj-))  •  Then,  curves 
of  larger  and  smaller  values  are  tried  in  order  to  widen  the  acceptable  set. 

In  this  way  the  bounds  shown  in  Fig.  2.5  were  obtained.  It  is  required 
(recall  1.3)),  that 


a33(oi)  <  | T33 ( jus)  |  $  b33(w)  .  (2.9) 

This  approach  may  be  criticized  as  rather  crude  for  translating  time-domain 
into  uj-domain  tolerances,  but  we  found  it  highly  reliable  in  many  design 
examples  for  example  [8,9],  and  it  is  easy  to  execute.  Reasons  for  working 
in  the  frequency  domain  rather  than  in  state-space,  have  been  given  in 
[10,11]. 

In  this  mode  *  r3  *  0  ,  so  the  only  command  input  is  r3  of 
Fig.  1.1,  and  in  Table  1.1  only  D^3  ,  D23  ,  (A33  ,  B33)  are  applicable. 

Consider  D^3  written  as 
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maximum  of  one  degree  steady-state  value  of  6  ,  then 

b^(0)  *  1/ (57. 3)  (32.2)  =  (.5)10  ^  .  We  can  prescribe  this  value  for 

all  a)  .  Comparing  this  value  to  b^  in  Fig.  2.5,  clearlv  b^>>  b,^  * 
at  least  from  .=0  to  ~  20  .  But  if  it  will  help  in  the  design, 
can  be  decreased  right  down  to  *  This  will  make  it  harder  on 

the  design  of  ,  because  the  tolerances  on  I  are  Chen 

narrower,  but  it  is  later  found  that  the  demands  on  ,  L ^  are  much 
harder  than  those  on  . 

Returning  to  (2.10)  since  b^  >>  ^23  *-n  c^e  l°w  frequency  range, 

b^^  J  will  dominate  in  the  numerator  unless  iQ^/Q^I  >:> 

We  therefore  sketch  in  Fig.  2.6a,b  these  two  ratios  from  the  data  available 
in  Table  2.3.  Cases  1-5  suffice,  because  Cases  6-10  are  pairwise  very 
similar  to  Cases  1-5.  Clearly,  ^^IQ^/Q^I  dominates  at  least  up  to 
m  *  10  ,  so  in  this  w  range,  we  need 


1  +  L,  1=  |1  +  G.Q 


b33  Q11 

b13  Q13 


(2.11) 


which  is  ~  66db  +  (-50  db)  *  16db  for  w  6  [0,6]  .  It  helps  to  approximate 

|l  +  G^Q^I  by  |g^Q^jJ  ,  because  (2.11)  then  becomes 


|G1(»U  b33/b13  lQL3l 


(2.12) 


and  the  bounds  on  G^  are  easily  determined  in  this  range.  The  value 
(1500)  ^  was  chosen  for  b^  »  corresponding  to  1.5°  max  8  for 


A  of  lg  .  The  maximum  value  of  |Q, 
y  0  '  1 


.0375  for  Case  4,  giving 


55  for  IG, I  .  .  Case  9  was  chosen  as  the  nominal  plant  case  Q,~  , 

1  min  10 

with  the  nominal  loop  transmission  L^q  =  G^  •  It  is  convenient  to 

have  the  bounds  on  L^q  rather  than  on  G^  ,  because  in  the  higher  w 


In  any  case,  the  procedure 


range  one  can  no  longer  take  |l  +  L^|  =  . 

for  finding  the  bounds  on  L^q  in  the  low  w  range  is  fairly  clear 
from  the  above. 


As 


increases,  b  decreases  (see  Fig.  2.5)  and  so  does  'Q  /<) 

5  ’  U  I  <  I  J 


(see  Fig.  2.6b).  The  first  term 


23 


13 


:11  I 


— —  I  in  (2.10)  is  never  large,  so 


12 


the  demand  on  |l  +  L^|  in  (2.10)  steadily  decreases.  Eventually  |l  +  L^| 


may  be  quite  small  which  would  permit  -*■  -1  .  But  it  is  sensible  to 


impose  reasonable  stability  margins  on  the  loop  by  itself.  Physically, 

this  would  correspond  to  the  case  y^  =  d>  =  0  ,  y ^  =  0  and  the  only 


feedback  is  from  S  because  then  in  5 


.-1 


y_  (see  Eq.  (2.5)),  5^  =  y 


So  in  the  higher  w  range  bounds  were  found  on  such  that 


,-l 


|1  +  LjJ  £  2db  .  The  resulting  bounds  on  are  shown  in  the  Nichols 


shart  in  Fig.  2.7.  The  next  step  is  to  find  a  L^s)  *  such  that  L  Q(jw) 


satisfies  the  bounds  of  Fig.  2.7.  This  was  done  simply  by  cut  and  try. 


The  Bode  plots  of  L1Q(ju)  are  shown  in  Fig.  2.8,  and  L1Q(jw)  is  also 


shown  in  Fig.  2.7. 


The  resulting  | ( jto>  |  =  |  L,  0/Q 


10,x110' 


is  sketched  in  Fig.  2.13,  and 


it  is  seen  it  must  be  greater  than  one,  over  an  extremely  large  w  range. 
The  reason  for  this  is  evident  from  Fig.  2.7.  The  intrinsic  demands 


on  due  t0  c^e  specifications,  become  negligible  at  w  *  25.  However 


Liq(  jtu)  |  cannot  decrease  very  rapidly  vs  ,  because  its  phase  lag  is 


not  allowed  to  be  more  than  ~  130 


And  this  constraint  holds  until 


| Li© ( jn* )  |  has  decreased  to  about  -24  db  .  This  is  due  to  the  large  spread 


in  the  high  frequency  values  of  (see  Fig.  2.1).  The  present  design 


is  clearly  overly  conservative.  At  u>  =  8 


'10 1 


is  about  3  db  larger 


than  necessary.  From  15  to  -24  db  levels  there  is  a  change  of  39  db ,  at 


130 


an  average  slope  of  -  x  12  =  -8.7  db/octave,  requiring  39/8.7  =  4.5 


,4.5 


octaves,  giving  (8)  2  »  180  for  the  frequency  at  which 


'10  1 


-24  db 


-  29 


The  decrease  of  |l^q|  can  he  more  rapid  thereafter,  but  again  we  have 
been  rather  conservative. 

There  is  a  trade-off  between  the  complexity  of  the  compensation 
function  (the  number  of  its  poles  and  zeros),  and  the  b tndwidth  of  Gj  . 
It  is  clearly  possible  to  decreas  G^  more  rapidly  in  the  mid-frequency 
range.  Also  more  far-off  poles  can  be  inserted  in  the  higher  w  range. 
All  this  requires  more  careful  shaping  of  L^Cjui)  ,  with  more  poles  and 
zeros.  This  will  be  done  in  later  designs. 

It  is  also  noted  that  a  significant  reasona  for  the  large  bandwidth 
of  ,  is  the  large  gain  margin  (about  25  db)  needed  for  ,  because 

of  the  large  spread  in  the  high-frequency  value  of  (see  Fig.  2.1). 

This  spread  can  be  reduced  by  high-frequency  scheduling,  which  is  done  in 
Sec.  2.8.  The  analytic  expressions  for  the  L^^s)  ,  G^s)  etc  are  given 
in  Appendix  1  of  this  chapter. 

We  next  turn  to  the  design  of  ^20^'  From  Table  1.1,  the  only 


(since  only 

A  is 

y 

being  commanded)  is 

l>23  of  Table  1.1, 

Q22 

1  _1_  u 

<M 

CM 

O' 

13  q21 

1  +  b33 

1  q23  1 

( 0  TON 

5 

b23 

.ij; 

The  bounds  on  I^qO^)  are  obtained  in  the  same  manner  as  those  on  L^q  , 
with  the  result  shown  in  Fig.  2,9.  This  figure  has  precisely  the  same 
features  as  Fig.  2.7  and  for  precisely  the  same  reasons  -  especially  the 
even  larger  gain  margin  needed,  due  to  the  larger  spread  in  the  high- 
frequency  gain  factor  of  (see  Fig.  2.2).  The  I^gCju)  chosen  is 

shown  in  Fig.  2.10,  and  is  also  rather  conservative.  It  should  be  noted 
that  as  Q220  has  a  r*8ht  half-plane  pole  at  .0283  (see  Table  1.2, 

Case  9),  it  is  imperative  that  L^q  has  this  same  pole;  otherwise  there 


is  attempted  right  half-plane  pole  cancellation,  which  would  give  a 
dipole  in  the  right  half-plane,  and  the  system  would  then  be  unstable. 


2.4  Design  of  third  (A  )  Loop 

V 

It  is  recalled  that  condition  (1.14)  could  not  be  satisfied  in  this 
problem,  so  the  design  of  the  third  loop  is  performed  by  using  the  exact 
equation.  In  Fig.  1.1, 

T  =  (I  +  L)-1LF  ,  L  =  PG  .  (2.14) 

We  are  interested  in  the  actual  expression  for 


l3f3[(i+l1)(i+l2)-y12] 
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(1+L^  (i+l2)  (i+l3)-[y12(i+l3)+y13(i+l2)+y23(i+l1)  ]+r 


(2.15) 


.  Qliq22 
12  Q12Q21 


qiiq33 


Q22Q33 


13  Q13Q31  ’  23  Q23Q32 


(qi 9Q9 


12^23^31  ^21^32^13) 


Q11  q22  Q33  '  QoiQ-joQ! 


(2 . 16a-d) 


The  above  was  obtained  by  expanding  the  matrices  in  (2.14).  For  purpose 
of  design  of  L3  ,  it  is  desirable  to  rewirte  (2.15)  in  the  following  form: 


F33  W* 


33  L 


30 


L30  *  G3  Q330 


(2.17a,b) 


T(ju)) 


Cy13(1+L2)  +  y23(1+L1)"I'J  ^ 
(l+LjKl+i^)  -  y12  j1 


(2.18) 


Of  course,  <p  is  a  function  of  w  and  of  the  plant  parameters.  There  is  a 
different  \\>  value  for  each  of  the  ten  flight  conditions. 

The  procedure  is  to  find  the  templates  Tp  (w)  ,  of  ijj(jw)  i.e. 

To  (u)  ”  {'Kit))}  .  Gone  templates  are  ohot/n  in  Fig.  2.11.  It  is  recuircd 
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that 

a33(uJ)  S  |l33(jw)|  <  b33(^)  (2.19) 

of  Fis.  2.5.  The  function  of  I  is  cn  guarantee  that  over  the  range 

•  >  f  : 

b33 

$  log-22  (2.20) 

a33 

The  Nichols  chart  is  the  perfect  design  tool  for  this  purpose  because 

it  contains  loci  of  constant  |L/(1+L) |  .  Here  L3q ^ ^  takes  the  place 

of  L  .  As  there  is  no  uncertainty  in  L3Q  ,  the  template  (set  {L^/'f}) 
of  L-jq/Y  is  si®ply  the  template  of  1/Y  shifted  by  the  value  of 

L30(jw)  •  One  positions  the  template  of  1/Y  on  the  Nichols  chart  until 
(2.20)  is  barely  satisfied,  giving  a  bound  on  L3q  .  The  design  procedure 
is  in  fact  precisely  the  same  as  in  single-loop  design,  described  in  detail 
in  [3], 

The  bounds  on  L3q  thus  obtained  are  shown  in  Fig.  2.11,  together 
with  the  L3Q(»  chosen.  A  Bode  plot  of  L3q  is  shown  in  Fig.  2.12. 

There  is  here  too  the  same  phenomenon  as  in  the  design  of  L  n  ,  L.n 
(Figs.  2.7,  2.9).  In  the  higher  ui  range,  an  additional  constraint  of  the 
form  1 1  +  l>3 1  1  <  some  specified  value  becomes  more  demanding  on  L3 
than  (2.19).  Also,  the  large  spread  in  the  value  of  the  high-frequency  gain 
factor  of  Q33  (see  Fig.  2.3),  forces  | L3q ( jou)  |  to  be  decreased 
slowly  until  |L3q|  ~  -30  db  .  However,  since  |Q33q(ju)  |  is  considerably 

larger  thatn  | Q ,  |  and  I ^22  ®  ’  Cbe  demand  on  ^3  =  ^30^330  ^  aS 
considerably  less  than  on  ,  G3  ,  as  seen  in  the  Bode  plots  of  the 
(asymptotic)  G3  ,  ,  G3(jw)  in  Figs.  2.13,  and  2.14. 


A  log 
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Proper  choice  of  only  guarantees  that  the  change  in  I T 3 j ( .1  '■**)  I 

is  within  the  amount  allowed  by  (2.20).  The  role  of  in  Fig.  1.1, 

is  to  position  | T ^ ^ ( .i us )  |  within  the  tolerances  of  (2.19).  For  example 
if  .1^(10)  =  -23.5  db  .  b  j  j  ( 10)  =  -*.5,  and  if  .  (L  /  • )  /  |  l-(  L  /  • )  j 
in  practice  ranges  from  -10  db  to  +  2  db  at  u>  =  10  ,  then  it  is 
necessary  that,  -13.5  £  j ( j  10 )  [  s  -10.5  db  .  Bounds  in  |F^(jui)  |  are 
obtained  in  this  manner  and  then,  by  cut  and  try,  F^(s)  was  selected 
which  satisfied  these  bounds  (see  Appendix  1  for  F^(s)). 

Note  that  T^Cjw)  =  F^(ju)  over  the  significant  frequency  range, 
because  the  bandwidth  of  L^(jcj)  is  so  much  larger  than  that  of  T^O'-o)  • 
This  happens  in  systems  in  which  large  parameter  uncertainty  demands  large 
loop  bandwidths. 


tf 


2.5  Simulation  results.  Design  1 

The  system  was  simulated  on  the  digital  computer  for  Ig  step 


command  of  .  The  responses  are  shown  in  Figs.  2.15a-c  in  degrees 


for  2,  and  g's  for  A^  .  The  Ay  responses  for  the  five  cases  were 


very  similar,  so  result  in  almost  a  single  curve  in  Fig.  2.15c.  The 
other  five  cases  (6,  7,..., 10)  are  pairwise  very  similar  to  these  first 


five  and  are  not  shown.  The  control  surface  deflections  5  ,5  ,5 

r  a  c 


in  degrees  are  shown  in  Fig.  2.16a-c.  The  control  surface  rates  are 


shown  in  Figs.  2.17a-c.  The  actual  control  surface  limits  are  respect¬ 


ively  i30°  ,  ±20°  ,  ±25°  for  5,6  5  ,  so  Case  3  is  well 

r  a  c 


below  saturation,  while  Case  1  barely  causes  canard  saturation  (25.7  ) 


The  other  canard  values  are  33°  ,  120°  ,  54°  for  Cases  2,  4,  5 


respectively,  so  that  to  avoid  this  amplitude  limiting,  A^  commanded 


must  be  reduced  to  . 76  ,  .21,  .46  g's  respectively.  In  all  cases 


-  -  A.  -r-  £-  -1  ■ 
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Che  canard  is  the  limiting  surface. 

The  race  limits  are  120,  56,  100  deg/sec  for  6,5.6 

r  a  c 

whereas  Fig.  2.17  shows  Chat  in  Case  A,  rhe  canard  rate  is  j  semt-wha r 
greacer  limiting  factor  than  the  canard  amplitude,  limiting  A  to 
.17  g  for  this  Case  4.  This  can  be  improved  by  slowing  up  T^  .  If 
the  step  response  of  is  slowed  up  by  X  (the  time  scale  in  Figs.  2.15 

multiplied  by  X  ,  then  the  6’s  will  be  reduced  by  this  factor  very  closely, 

because  as  noted  in  Sec.  2.4,  ^33^“^  *  F^jw)  over  the  significant 

range  of  ^3^^)  •  This  slowing  up  by  a  factor  X  is  achieved  by 
multiplying  each  of  the  poles  by  X  .  The  reason  is  that 

with  its  much  smaller  bandwidth  than  the  ,  dominates  the  system 
response  to  commands  (see  Sec.  2.10). 


2.6.  Responses  under  amplitude  and  rate  limiting 

The  control  surfaces  amplitude  and  rate  limits  were  incorporated  into 
the  system  model,  in  both  cases  as  simple  algebraic  nonlinearities.  If 
the  rate  saturates  first,  the  amplitude  can  continue  to  increase  until 
it  becomes  saturated,  and  vice  versa.  The  commanded  input  was  in 

each  case  made  so  high  as  to  force  the  system  to  saturate.  The  resulting 
5^  responses  are  shown  in  Figs.  2.l8a-c,  and  the  6^  responses  in 
Figs.  2.19a-c.  In  all  cases  the  canard  rate  was  the  first  to  saturate  and 
the  canard  amplitude  eventually  also  saturated.  Figs.  2.20a-c  show  the 
6»  ,  Ay  responses  under  these  very  hard  saturated  conditions.  In 

all  cases,  the  responses  are  very  nice,  with  little  overshoot.  Note  that 
in  all  cases  except  3,  3.1  g's  were  c snmanded  and  9.3  g  for  Case  3,  so 


as  to  deliberately  drive  the  system  into  very  hard  saturation. 


2.7  Design  two 


Several  more  designs  were  made,  primarily  in  order  to  decrease  the 

bandwidths  of  G^  ,  G.?  .  In  design  2  some  sacrifice  was  made  in  rhe  final 

values  of  „•  and  :  when  A  is  commanded.  Inis  is  eviuenr  from  the 

v 

smaller  values  of  G^(0)  ,  G^CO)  in  this  design  as  compared  to  Design  1. 
Also,  an  additional  two  far-off  poles  were  assigned  to  G^(s)  ,  G^(s)  in 

order  that  they  may  reach  zero  db  magnitude  at  a  considerably  smaller 
frequency  than  in  the  first  design  -  compare  Designs  1,  2  in  Figs.  2.13-4. 

Otherwise,  the  design  procedure  is  identical  to  that  of  Design  1,  and  is 

therefore  not  described.  Because  of  the  faster  reduction  of  |l  |  ,  |L0q| 
a  greater  burden  was  imposed  on  ,  so  in  Fig.  2.13  has  a  larger 

bandwidth  for  Design  2,  than  in  Design  1.  But  this  is  acceptable,  because 
the  demand  on  is  so  much  less  than  on  ,  G2  .  The  Bo<*e  pl°ts  °f 

Liqi  L20  ,  L3Q  are  shown  in  Fig.  2.21a-c.  The  computer  simulation  results 
are  shown  in  Fig.  2.22  It  is  seen  that  the  6,  «J>  responses  are 

as  expected  somewhat  poorer  than  those  in  Figs.  2.15-2.17  for  Design  1. 


2.8  A  design  with  scheduling 

One  of  the  reasons  for  the  large  G^  ,  G2  bandwidths  of  Design  1, 

noted  in  Sec.  2.3,  was  the  large  gain  margins  needed  (Figs.  2.7,  2.9)  for 

L^q  >  ^2q  .  This  was  due  to  the  large  uncertainty  in  the  high-frequency 

gain  factors  of  ,  Q22  (Figs.  2.1,  2.2).  A  significant  decrease  in 

the  G^  ,  G2  bandwidths  should  be  possible,  if  there  can  be  a  significant 

decrease  in  this  high-frequency  uncertainty.  We  therefore  examine  Eqs. 

2.1-8  for  lim  Q  .  The  results  are  (letting  q  =  Q  ^  ) : 

uv  uv  uv 
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Thus  all  Che  Q^v  are  inversely  proportional  to  q  as  s  -*■  00 
Hence,  the  part  of  the  uncertainty  due  to  the  q  variations  could  be 
eliminated,  if  q  can  be  accurately  measured  rapidly  enough.  We  assume 
that  q  can  be  measured  with  t  20%  error,  i.e.  the  measured  value  may 
in  fact  be  between  .8  and  1.2  cimes  the  correct  value.  The  sero  u)  values 


-36- 


of  the  Q  are  independent  of  q  ,  so  the  measured  q  are  to  be  used 
only  in  the  higher  u>  range.  Accordingly,  we  insert  in  front  of  the 
plant  a  diagonal  scheduled  compensate *■,  each  of  whose  elements  i< 


J(s)  = 


1  +  s  q 


Ts  1200 


(2.23) 


where  q  is  the  measured  value  of  q  .  One  should  reallv  experiment 
m  ■ 

with  the  value  of  .  But  our  primary  purpose  here  is  to  demonstrate 

how  well-known  techniques,  such  as  scheduling,  can  be  incorporated  into 
our  design  procedure. 

This  incoporation  of  scheduling  into  our  design  technique  is  very 
simple.  Let  the  new  plant  be  P^  =  PJ  (s)  ,  so  the  new 


n 


(2.26) 


with  J(s)  a  function  of  q  and  therefore  of  each  flight  condition.  The 

procedure  is  now  exactly  the  same  as  before,  but  Q  (functions  of  the 

flight  condition)  replace  Q  .  But  now  there  are  thirty  cases  because 

for  each  flight  condition  q^  in  (2.25)  may  be  .8  q,  q  or  1.2  q. 

Some  of  the  templates  of  Quun  are  shown  in  Figs.  2.23-2.25,  for 

,  ^22n  *  Comparing  with  the  templates  of  Q  (old)  in  Figs.  2.1-3, 

it  is  seen  that  there  is  little  difference  at  low  ^  values,  but  with 

significant  changes  of  the  templates  at  higher  u>  values.  The  new 

template  Q  at  large  u  (100)  is  actually  larger  in  size  that  the 

unscheduled  Q, ,  (100)  .  However,  the  nominal  (Case  9)  Q, ,  is  near 
11  1  In 

the  top  of  its  template,  rather  than  the  bottom  as  it  was  for  the 
unscheduled  case.  So  the  resulting  bound  on  L ^  at  high  w  is  not 
so  far  down,  we  don't  need  as  large  as  gain  margin  and  \  G  (j  100) 


for 


the  scheduled  design  is  much  smaller  than  for  Design  1  or  2.  In  finding 

these  templates  of  Q  ,  three  different  values  were  assigned  to  q 
r  uun  m 

of  (2.25)  for  each  flight  condition.  These  were  q  =  .  8q  ,  q  and  L.2  q. 

m 

The  bounds  on  trhe  Nichols  chart  analogous  t»>  those  in  the  rirsr 

design  (Figs.  2.7,  2.9,  2.11)  are  shown  in  Figs.  2.26a-c,  which  also 

include  L^qOuj)  .  Bode  plots  of  the  L  are  shown  in  figs.  2.27a-c 

Digital  simulation  results  are  shown  in  Figs.  2.28a-c.  For  each  flight 

condition  three  runs  are  shown,  one  for  the  case  q  =  . 8q  ,  the  second 

m 

if  happended  to  be  q  ,  and  the  third  for  q^  =  1.2  q. 

Figs.  2.13,  2.14  show  that  a  scheduled  design  can  be  more  economical 
in  G  bandwidths  required.  We  have  probably  overestimated  the  measure¬ 
ment  error  i.e.  the  allowance  £  [0.8q,  1.2q]  is  probably  much  more  than 
the  actual  error.  Also  experimentation  in  the  value  of  used  in  (2.25) 

will  lead  to  better  results,  as  would  a  more  sophisticated  J(s)  function 
of  (2.25).  In  any  case,  our  main  point  was  to  demonstrate  how  easily 
scheduling  may  be  incorporated  into  this  new  design  technique.  Also  how 
the  design  technique  reveals  the  frequency  ranges  in  which  scheduling  is 
worth  implementing. 


2.9  Nondiagonal  G 

The  three  designs  presented  so  far,  all  have  the  shortcoming  discussed 
in  Sec.  1.4.  The  are  the  products  of  a  design  technique  which  requires 
loop  transmissions,  even  when  there  is  no  uncertainty.  The  remedy  was 
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Tn  our  problem  only  A  is  commanded  so  there  is  need  onlv  for 

y 

h  .  .  The  others  of  diagonal  elemenrs  can  be  made  Zt*rn.  In 

L  >  *  —  - 

•4ener.it ,  the  »- 1  i-men  C  s  of  H  may  he  functions  or  r  requeue  v .  hue  ,(s  rnr 

primary  purpose  is  to  demonstrate  the  new  technique,  we  shall  concentrate 


on  u)  =  0  ,  and  indicate  the  procedure  which  can  be  followed  for  u>  in 
general. 


Recall  from  Sec.  1.3,  (Eqs.  2.10,  2.13)  that  the  bounds  on  , 

are  determined  by  the  requirements  that 


!i  +  llI  i 


23  |  Q, 


1  +  L, 


33  I  Q. 


g22  |  222  | 

’  13  Q21  +  b33  |  Q23  | 


(2.27a,b) 


Recall  also  that  it  was  shown  in  Sec.  1.3  that  the  second  terms  in  the 

above  involving  ,  Q23  easily  dominated  the  right  sides  of  (2.27a,b) 

up  to  about  10  'V*  15  rps.  Hence,  it  is  desirable  to  make  the  new  , 

Q_.  very  large.  Let  the  new  P  be  P  =  PH  ,  so  the  new 
z  o  n 

^1/Quvn^  H  P  =  H  t/Quv^  •  LeC  hll  =  h22  h23  =  1  ’  hl3  =  a’ 

h23  =  6  ,  the  others  zero  and  then 


Q13n  *  (P13  "  aP33)  ’  Q23n  =  [P23  "  SP33^ 


( 2 . 28a-c) 


If  the  b^^  ,  b^^  values  were  the  same  for  all  flight  conditions,  then 
a  reasonable  objective,  at  any  uj  ,  is  to  choose  a  ,  S  to  minimize 


max  I  P^3  -  ctP  ^  ^  I  over  various  flight  conditions.  However,  in 
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practi re  ciie  maximum  A^  commanded  is  generally  different.  Hence,  if 
sav  a  constant  maximum  3  angLe  is  tolerable  independent l v  of  A 

y 

c  .'raroanced  .  then  b,  ^  sh*>ti  l.i  be  inverse  I  v  nronort  ien.i  i  t  ■■■  A  .  ■'mm.spd.  -  ; . 
Til  is  was  assumed  here,  so  an  estimate  was  needed  or  che  maximum  A 

y 

commanded  for  each  flight  condition.  The  estimate  was  made  approximately 
on  the  basis  of  the  A^  needed  to  achieve  amplitude  saturation  in 
earlier  designs  with  lg  the  maximum  value: 


Cases  1,6  -  lg  ;  Cases  2,7  -  .6g  ;  3,8  -  lg 

•'-,9  -  .25g  ;  5,10  -  .5g  . 


(2.29) 


Accordingly,  we  seek  to  minfmax  |  W^(P^  ~  ^ )  |  ]  ,  where 

*  ^33^13  *  wh°se  values  are  now:  2000  ,  1200  ,  2000  ,  500  ,  1000 
respectively  for  the  cases  (1,6),  (2,7) ,...,  (5,10)  .  The  result  is 
a  *.571  ,  giving  5.64  for  the  min.  max.  value,  which  is  15  db  needed 
for  G^  at  u>  *  0  ,  instead  of  55  (=  34.8  db)  neeeded  previously. 
Similarly,  the  optimum  8  is  found  to  be  .347  ,  with  only  12  db  needed 
for  at  a)  =  0  .  In  principle,  this  process  can  be  continued  as 

a  function  of  frequency,  but  was  not  done  here  because  there  is  not 
much  change  in  the  above  parameters  up  to  w  =  10  .  So  a  matrix  of 
constants  is  used  here  for  H  . 
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Some  of  these  new  Q, ,  ,  templates  are  shown  in  Figs.  2.29a,b. 

1  In  22n 

They  are  somewhat  larger  than  those  of  the  diagonal  G  design,  in 
Figs.  2.1.  2.2.  Bur  the  much  greater  reduction  of  off-diagonal 
elements  in  (2.28a,b)  leads  to  the  nondiagonal  design  being  significantly 
more  economical  than  the  diagonal.  However,  this  point  should  not  be 
overlooked.  The  choice  of  H  should  consider  the  effect  on  the  uncertain¬ 
ty  rangers  of  the  new  Quun  •  A  more  sophisticated  optimization  technique 
is  needed  for  this  purpose.  The  new  bounds  on  L^q  ,  L^q  are  shown  in 
Figs.  2.30a,  b  and  Bode  polts  in  Figs.  2.31  a,  b.  The  new 
values  are  shown  in  Figs.  2.13  ,  2.14  where  it  is  seen  that  they  are 
significantly  more  economical  than  those  in  previous  designs. 

Simulation  results  are  shown  in  Figs.  2.32a-c:  B  ,  0  ,  A  in 

y 

Fig.  2.32a  ;  5  ,  6  ,6  in  Fig.  2.32  b  and  6,6  ,6 

r  a  c  °  r  a  c 

in  Fig.  2.32c.  In  these  simulation  results,  we  deliberately  commanded 

o 

sufficient  A  so  as  to  obtain  the  maximum  available  25  canard  value, 

y 

except  for  Cases  3,  8.  The  commanded  A^  values  are  given  in  brackets 
in  the  A^  responses  in  Fig.  2.32a.  Exceedingly  small  B  ,  0  values 
are  achieved  in  this  design,  part  of  the  reason  being  the  smaller  A^ 
being  commanded  for  6  of  the  10  cases.  There  is  no  rate  limiting  in 
all  of  these  cases.  The  A  commanded  are  (in  g's): 

y 

Case  1,  lg  ;  2,  .8g  ;  3,  lg;  4,  .23g  ;  5,  .5g  ;  6,  . 73g  ; 


7,  . 6g  ;  8,  lg  ;  9,  .17g  ;  10,  .35g  . 


(2.31) 


From  Fig.  2.13,  2.14,  it  is  seen  that  this  is  the  most  economical 
design  so  far.  It  could,  of  course,  be  made  even  more  economical  by 
incorporating  "scheduling"as  well. 
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2.10  Trade-off  between  speed  of  response  and  rate  limitin 


APPENDIX  2.  I 


COMPENSATION  FUNCTIONS  FOR  VARIOUS  DESIGNS 

1  (Sects . 2 . 3-5) 


_ (260)1Q4(s+7) (s+60) _ 

(s2  +  .706s  +  27. 3) (s  +  25) [ s 2  +  400s  +  (400)2] 

(52) IQ4 (s+7) (s+60) _ 

(s  +  25) [s2  +  400s  +  (400) 2 ] 

_ (5)lQ5(s+10)  (s+180) _ 

(s  -  . 0283) (s  +  l.l)(s  +  60) [s2  +  280s  +  (350)2] 

-(4.08)  1Q4(s  +  180)  (s+10) 

(s  +  60) [s2  +  280s  +  (350) 2 ] 

28, 100(s2-H2s-H0Q)  (s+200) _ 

(s2  +  4s  +  6.25) (s  +  25) (s  +  60) (s  +  600) 

(22.5)  (s-HOO)  (s2+12s+1Q0)  (sH-200) _ 


(s  +  4s  +  6. 25)  (s  +  25)  (s  +  60)  (s  +  600) 


100 _ 

(s  +  5)(s  +  20) 
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Design  2 


-5.4(1  +  asl(ljjg) _ 

,  +  S_  J,  l  S  \  f  S  V  ]" 

1  40  )|_1  \15dJ  \i50J  J 


34.9  (l+.ls)  1  + 


^  ,  +_s_+^02]2 

v1  40 J  j.  1  350  \,350y  J 


8  (l+.Ols) 


V1  +  15o) 


1.®.  (  -\  a.  !z _ A  ( i  _i_  s  \ 

1  6  V  1  30 y  V1  60C  J 


3  (s  +  5)(s  +  20) 


Scheduling  design 


L10  * 


(4 . 14) 105(s2+12s+100) (s+75) 

(s2+. 706s+27. 3) (s2+45s+625) (s2+136s+(170) 2]  2 


0.84  x  105  (  32+12s+ 100) ( s+75 )( 1+s / sQ  ) 
1  (s2+45s+625) (s2+136s+170)2]2 


„  (1.029) 10^3+5) (s+30) (s+60) _ 

20  (s-.0283) (s2+48s+1600) [s2+170s+(170) 2] 2(s+l .1) 


-(4.195)10  (1+. 2s) (s+30) (s+60) ( l+s/sQ) 


2  (s2+48s+1600)  [s2+170s+('170)2]2 


(4.8)104(s+25) 


L 


30 


(s+3) (s+IO) (s+400) 


38.h(s+35) (s+100) 
(s+3) (s+10) (s+400) 


F 


3 


100 

(s+5) (s+20) 2 


L. 

1 


G.Q..(old) 


1 _ 

(1  +  s/7.56) 


Nondiagonal  design 


„  _ (93.6)1Q6(s+2.3)  (s+14) 

10  (s+6) (s2+0. 296s+5 .4) [s2+83s+(83) 2 ] 2 

(44.6)106(s+2.3) (s+14) 

1  (s+6) [s2+83s+(83)2]2 

=  _ (50)1Q8(s+4)  (s+40) _ 

0  (s-.0264)(s+. 957) (s+10) [s2+200s+(200)2]2 

„  -5.2  x  108(s+4) (s+40) 

2  (s+10) [s2+200s  +  (200) 2] 2 


G 


3 


1.2 


1  +  1.67s 


F 


3 


100 


(s+5) (s+20) 


oj  ( RP3 ) 

O  e 


12  5  10  20  50  100  200  50' 


fig-  2 • 6a , b  Bode  plots  (asymptotic)  of  | Q  /Q  |  |Q  /Q 

ii  12  ll'  13 


es 


omparison 


to  (RPS) 


w(RPS) 


Fig.  2.15 


Fig.  2.16 


Fig.  2.17 


Fig.  2.15-2.17  Computer  simulations  of  design  1,  for  A  command  of  lg. 

y 

(15)  8  ,<t  (deg . )  ,  (16)  5  (degrees),  (17)  A  (g's). 


Fig.  219 


>  - 


•ig.  2.20 
Figs.  2.18-2.20 


Fig.  2.19 
* 

Responses  due  to  very  large  A  commands  causing  saturation 
9.3g  for  case  3,  3.1g  for  others. 
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Fig.  2.24  Templates  of  Q7~ 


scheduled  design 


OB 
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scheduled  design 


Fig.  2 . 28a-c .  Scheduled  design,  simulation  results  for  A  command  =  Ig. 

(a)  8 ,4>(deg.)  >  A^(g's)  (b)  ^(degrees)  y  (c)  ^(deg. /sec .) 


diagonal  compensation  design. 
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CASE  3 


(cont.)  Simulation  results  for  design  l  (6  <j>  a  ) 

’  ’  v 


with  F  -  F  -  36/ (s+3) 
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Fig.  2.33  (cont.) 

Simulation  results  for 

1  (-,:.Av)  with  F3  =  F 

=  36/(s+3) (s+12) 


design 

3  b 


Fig.  2.34  Simulation  results  for  design  1  (6.»5.)  with  F  *  F 


36/ (s+3) (s+12) 


DEGRE 


Fig.  2.34  (cont.)  Simulation  results  for  design  1  (6  ,5^)  with  F^*  =  36/ (s+3) (s+12) 
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CHAPTER  3 


VERTICAL  TRANSLATION  MODE  (c^) 

C.:. ipc^r  3  is  devoted  to  a  longitudinal  node  in  which  verr:  .. 

velocity  (Vz>  and  pitch  angle  (?)  are  the  output  variables  of 
interest.  The  plant  inputs  (5^  control  surfaces)  are  the  elevator 
and  flaps  (Fig.  3.1). 


3. 1  Plant  equations 

The  plant  equations  are  gotten  from  [12,  pp.  292-8],  modifi  .or 
inclusion  of  flaps.  The  parameters  are  defined  in 


y  *  [u,  w,  a] 1  ,  6  =  [6g,  5f  ]  • 
Ay  *  B6  ,  z  *  [a,  0]  -  My 


(3.1a,b] 
(3. 2a, b) 


X 

w 

g  cos  a 

s-Z 

w 

-Vs 

-M 

s(s-M  ) 

w 

q 

_CD6 

-CDf 

m 

m 

~^L3 

^Lf 

m 

m 

cCM5 

CCMf 

I 

y  • 

I 

y 

(2.3) 


(3.4) 


0 


3 


12.6s  -s V 


I 


! 

4 

I 

\ 

b 

' 


M 


0 


0 


1 


Since  P  is  defined  by  z  =  P6  ,  we  have 

P  =  MA_1B  =  [P. . ] 

J-J 


-  -pSVCD/m  , 

X  = 
w 

pSV(CL-CDct)/2m 

=  -PSV2Cri  ./2m 

D6 

xf = 

-  SV2CDf/2m  , 

=  pSc2C  / 41  , 

Md  y 

PSVcC  / 21 

Ma  y 

»  -pSVC  /m 

L 

Z  =* 
w 

-PSV(CLa4CD)/2m 

-  -pSV2CL5/2m  ,  Zf  =*  -pSV2CLf/2m  , 

“  0SVc2(V4ly  *  M5  *  PSv2cCM6/2Iy  * 

*  -0ScCL./4m  ,  Mf  -  SV2cCMf/2ly 

Z  »  -  SVcCT  /4m 

q  Lq 


11 


12 


21 


pSV 


2  dec  A 

,2 


pSV 


2  dec  A 


P$V2 
2  dec  A 


n?u(l)  s4  +  Pu(2)s3  +  pu(3)s2 
[p12(l)s4  +  p12(2)s3  +  p12(3)s2 

[  P21(l)s2  +  p21(2)s  +  p21(3)  ] 


'22 


2 

PSV 

2  dec  A 


[p22(l)s2  + 


p22(2)s  +  P22(3)] 


z  =  My  -  MA  ]B 


+  Pu(4)s] 
+  p12(4)s] 


Pll(1)  =  b21  +  12 • 6b3i 


pu<2)  -  -Ku»21-b21(Xu+«,)-12.6b21.32-n.6b31<VZw) 


p  O)  =  b  ,  ( a  >1  +12.6.1  J  )+b  iX  M-*-i.2.oa.  X  rj,  V  - 
11  11  21  q  21  32  21  u  q  32  u  j 2 

+b-,(12.6X  Z  +VZ  -12 . 6a  0a  ) 

31  u  w  w  12  21 


Pll(4)  *  -bHa2ia32V“b2ia32XuV+b31(a2lai3“VXuZw+ai2a21V) 


t,  3 . 9a-d) 


P12(l)  =  b^,2+12.6b32 


p12<2)  ~  "bi2a21~b22(Xu+Mq)_12'6b22a32"12‘6b32(Xu+Zw) 


Pi2(3)  -  b12(a21Mq+12.6a21a32)+b22(XuMq+12.6a32Xu+a32V)+ 
+b32(12-6XuVVZW-12-6a12a21) 


(3 .10a-d) 


P12(4)  =  -bi2a21a32V"b22a32XuV+b32(a21a13-VXuZw+a12a21V) 


P21(1)  ’  b31 

P21  (2)  =  -b2la32-b31(W 


p  (3)  -  bua2ia32+b21a32Xu+  b31(XuZw"a12a21) 


(3 .lla-c) 


Q ( 1 )  -  p11(1)p22(1)-p12(1)p21(1) 

Q(-)  =  P11(2)p22(l)+plia)p:2(2)-p12(2)p21(l)-p,  ,(l)Poi 
Q (3)  =  p11(3)p22(l)+p11(2)p22(2)+p11(l)p22(3)-p12(3)p21 

“P12^2^P21^2^"P12^1^P21^3^ 

Q(4)  -  p11(4)p22(l)+Pll(3)p22(2)+p11C2)p22(3)-p12(4)p21 

-pl2(3)p2i(2)-pi2(2)p2i(3) 

Q(5)  =  p11(4)p22(2)+pu(3)p22(3)-P12(4)p21(2)-p12(3)p21 

Q(6)  «  pi;L(4)p22(3)-p12(4)p21(3) 

QQ  =  Q(1)s6+Q(2)s5+Q(3)s4+Q(4)s3+Q(5)s2+Q(6)s 


3.2  Flight  conditions 

As  in  Chapter  2,  there  are  five  different  combinations  of  mach 

and  altitude.  The  aerodynamic  coefficients  C,  .  ,  Cw.  ,  C„.  of  the 

Li  ’  Mi  Di 

Tables  below  are  functions  of  the  elevator  and  flap  angles,  and  were 
supplied  for  five  different  values  of  elevator  angle  and  for  two  of 
flap  angle.  Since  these  surfaces  work  generally  in  opposite  directions, 
only  such  combinations  were  used,  giving  a  total  of  five  cases  for 
each  of  the  five  flight  conditions  -  a  total  of  25  different  sets  of 
plant  parameter  values.  Some  of  the  data  common  to  both  the  lateral 
and  longitudinal  modes,  has  already  been  listed  in  Tab  lie  2.2,  so  is 


not  repeated  in  Table  3.1. 


.(1/ 


(1/ rad. )  -3.93490  -4.15044  4.96202  -3.94221  -5.10697 


1002  .00007  .00072 


cases  a-e  are  due  to  the  following  5  combinations  of  Elev.  and  Flap 


3 . 3  Transfer  funct i ons  for  25  flight  conditions 


The  fransfer  functions  for  the  25  flight  conditions  of  Sec.  3. 2, are 


in  Appendix  3.1.  Those  of  the  elements  ptJ  of  P  ■  [p,.]  ,  as  well 

ij  iJ 


_  1 


as  chose  or  Cne  inverse  i 


f  1  1  are  p  reser.  cev. 

1J 


e .  . 

iJ 


refers  to  the  high-frequency  gain  factor.  Thus  p  -*■  k^/s  "  as 

s  -*•  “>  ,  and  "gain"  refers  to  k.  .  for  p..  or  Q, .  ,  as  the 

iJ  ij  ij 

case  may  be. 


3.4  Problem  statement  and  response  specifications 

t 

The  original  problem  statement  was  "Vertical  velocity  V  ■  /  A  dt 

z  0  n 

is  commanded.  This  is  to  be  achieved  with  very  small  effect  on  the 
pitch  (9)  angle  -would  like  0.1°  max.,  but  .5°  may  have  to  be 
accepted.  The  elevator  and  flap  are  needed  for  this  purpose.  For 
test  results  and  responses  see  "a2  mode  response"  (Fig.  3.2b  here)  — 
try  for  steady  state  achieved  in  <4  seconds,  hopefully  in  2  seconds. 

Re  maximum  values,  see  how  much  can  be  done.  About  25ft/sec.  rate  of 
climb  was  about  the  maximum  achieved  in  test  results,  but  more  would 
be  desirable". 

Suppose  we  allow  .165°  for  9  due  to  25  ft/sec.  Then  b^  of 
Chapter  1  (the  equivalent  of  Table  1.1  for  n  *  2,  and  of  Eq.  1.2)  is 

b21  *  57TT^25  ’  (1*15)10"4  ("78<7db)  •  (3.16a) 

To  find  («u*bu).  suppose  we  use  a  critically  damped  second-order 
model,  and  ask  for  90%  response  in  2.8  seconds.  From  the  response  curves 
at  4*1  e.g.  [13, p.  200],  this  corresponds  to  2.8  *  4  ,  so  u)n  *  1 
We  shall  use  w  ■  1.4  ,  i.e.  the  T  (s)  model  is 


(3.16b) 


Tn(s) 


s2+2(1.4)s+(1.4)2 


and  shall  use  |T^(jo>)|  =  b^(“J)  ,  because  preliminary  examination 

indicates  that  as  in  Chapter  2,  iL^Cj-'.h  (here,  ct  ,L  (j^);  there) 
will  be  very  large,  so  A  T^^(joj)  due  to  varying  flight  conditions  will 
be  snail  over  the  important  <u  range.  Accordingly,  b^(ui)  ■  ! T  < jai>  | 
is  sketched  in  Fig.  3.2a. 


3.5  Bounds  and  design  of  I^gCjw) 

The  equivalent  (for  n  =  2)  of  V  of  Table  1.1  is  used,  written 


Q22  I 


|l  +  L  | 

2  b21  1  Q21 


G2Q22 


(3.17a) 


In  the  low  id  range,  |l2I  >>  1  so  write  |l  +  L2 1  =  | L2 1  =  I G2Q22 1  » 

so  (3.17a)  can  be  written 


Ig2(j«)|  * 


“li 

b2iq21(Ju) 


(.86)10  b. 


(3.17b) 


Clearly,  at  any  given  u  ,  min  | Q2i I  "*®eded,  so  we  prepare  asymptotic 

plots  of  |Q2^(jiu)  |  .  This  was  done  only  in  cases  la  to  5a,  shown  in  Fig.  3.3. 
Me  also  want  to  check  when  the  approximation  |l  +  G2Q22|  *  |g2Q22I  is  no 
longer  valid,  so  similarly  prepare  asymptotic  plots  of  I Q22  ^  *  shcwti  in 

Fig.  3.4.  Case  4a  was  chosen  as  the  nominal  case.  With  these  sketches,  it 
was  relatively  easy  to  prepare  Table  3.2  based  only  on  Cases  la-5a. 


Calculation  of  bounds  on  Lon(j“>) 
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For  j  HO  ,  we  use  the  constraint  |l  +  L,q  \  1  <  5  db  ,  and 
use  the  templates  of  to  obtain  the  bounds  on  L^Cjuj)  ,  which  are 

shown  in  the  Nichols  Chart  in  Fig.  3.5.  Some  templates  of  Q9^  needed 
for  tills  p-;n>’ae,  arc  shown  in  Fig.  3 .  :  ■  ■_ 

half-plane  pole,  which  must  be  present  in  L9q  ,  so  Arg  L9q  =  -180° 
at  w  *  0  .  By  cut  and  try  a  l^gCju)  which  satisfies  the  bounds 
of  Fig.  3.5,  was  found.  Its  Bode  plot  is  shown  in  Fig.  3.’,  and  it  is 
also  shown  in  Fig.  3.5.  The  resulting  | G2 ( I  is  sketched  in  Fig.  3.8. 
The  dashed  lines  show  the  greater  reduction  of  ^(jto)!  vs  w  that 
may  be  achieved  (should  this  be  important)  by  adding  a  zero  at  -400  and 
a  pair  of  complex  poles  with  ?  =  .3  and  to  =  800  . 


3.6  Design  of  ,  F 


The  exact  equation  for  T^s)  is  used  to  determine  L1Q(s)  .  It  is 


F.LH1+L-) 

T  (3)  -  - =— i _ _ _ 

lr  '  (i+l1)(i+l2)-y 


Q11Q22 

Q12Q21 


G  Q 
i  Hi 


(3. 18a-c) 


This  can  be  manipulated  into 
* 

F1  L1  * 

Tii(s)  *  -Hr  •  l  -  - 

1  1+L,  i 


1  - 


1+L. 


(3. 19a, b) 


Since  (3.19a)  is  exactly  the  expression  for  a  single  loop  system  with 
plant 


'll 


‘11 


(3.20) 


1  - 


1+L. 


single-loop  synthesis  theory  [3]  may  be  used  to  find  the  bounds  on 


L10  "  G1Q110 
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1  '  l+i:o 


Q110  Q220 
Q120Q210 


( 3  . : !  ; .  b 


To  do  such  a  design,  the  templates  of  are  needed,  and  are 

k 

available  because  L^g  is  known.  Some  of  these  templates  of  Q  are 

shown  in  Fig.  3.9.  The  small  circle  marks  the  nominal  Q  q  ‘  or<^er 

* 

to  get  an  idea  of  how  large  in  =  G^Q^  must  be ,  in  order  to 

k 

handle  the  problem,  we  sketch  Q  ^ju)  in  Fig.  3.10,  and 

see  that  it  is  very  large  over  the  bandwidth  of  T^Qu)  .  A  Nyquist 

* 

sketch  of  Q11q  makes  it  obvious  that  it  has  a  net  right  half-plane 
pole?  excess  of  right  half-plane  poles  over  right  half-plane  zeros  =  1  , 


so  Q110(s)  was  factored  on  the  computer  and  found  to  have  a  pole  at 
3.94  (due  to  Q22Q) 

* 

There  is  no  need  to  find  bounds  on  L^q  ,  because  it  is  so  easy  to 
•  *  *  . 

have  |L^/(1+L^)|  almost  invariant  over  the  significant  ui  range  of 
T^(jw)  .  Rather  G^(jw)  is  simply  chosen  to  achieve  reasonable 
stability  margins  for  T^  over  the  various  flight  conditions.  The 

k 

correspondingly  L^g(ju))  chosen  is  sketched  in  Fig.  3.1i.  The  locus 

*  k 

of  the  .templates  of  *  G^Q^Cjw)  ,  is  shown  in  Fig.  3.12,  with  the 


triangles  marking  L^q  . 


Finally  F^(a)  is  chosen  to  achieve  T^(j^)  of  Fig.  3.3. 
The  compensation  functions  for  the  above  design,  are: 


•(1  +5^)0+  TTr2_  > 
3o  l^OO 


'■j  ,  ( s  J  — 


(1  +  70  )(1  +  30Q  )(1  +  3000 


G2(s)  = 


23°(1+  A-)(l+!)(l  +  8§) 


<1  +  26><1+i>[1+B0  +  (^ 


_\2, 

V350 J  j 


Fl(s) 


1  + 


2s 


+  (_s_r 

1.4  \l.Uj 


(3. 22a-c) 


3.7  Simulation  results 

Computer  simulations  of  the  design  are  shown  in  Fig.  3.13  for  commanded 
V ^  of  25  ft/sec.  and  it  is  seen  that  the  specif ications  of  Sec.  3.4  are  remark¬ 
ably  well  satisfied.  Also,  for  this  command  of  25  ft/sec.,  the  saturation 
value  of  6^(20°)  is  exceeded  moderately  only  for  Case  3  (mach  .9;  alt.  0  ft). 

The  rate  saturation  value  5^(56° /sec)  is  exceeded  only  for  Case  4  (mach  .6; 
alt.  30,000  ft.) 

Fig.  3.14  shows  the  faster  responses  achieved  by  using  the  faster  prefilter 
F.  »  1.4/(s  +  1.4).  In  both  cases,  the  results  are  highly  satisfactory. 


APPENDIX  3 . L 


Trans  ter  Functions  tor  25  Flight  Conditions 


"Gain"  refers  to  high-frequency  gain  factor. 


CiSr  1'  U,H>  '1-4';S  Ai-J-i 


•  i  .  j :  ■  =  -27  i* 

'  '  -  P  ■ 

I  “  ^  V  w  1  *  W  •  t  -  /  ) 

(“  .//-"..I  i.  ?  u  ) 

0  L  l  -  S  Jr  Pi: 

(  -  0  i  **  J4;  +C1(  CiC  ) 

(-Jj 155t-ul.-0, 7  opE -01  ) 

Gil-.  P  1  =  -I?*. 

7  -  c  ■  -  a  ;i  F  p  :  2 
<  -C,  ♦'JGc  +C1  ,  0:0  ) 

CjC  ) 

p  .  l  .  s  ■  p  _  p :  2  _ 

(*Ji4>r.*ul  i  Ci  u_  __  ) 

I  *  ji  vi  I 

GM\  P'l  = 

'  ~  P  J'. 

<-J,l?7E>Cl,  C  :  0  ) 

P‘L2S  JF  P»L 

("Oo^O^F+ulf  v- 1  C  ) 

< -Oi lPPc-Gl »-C»  755E-JI  ) 
Gilj  P22  =  -i , GO 

Z.5.'?S  G F  °i2 

(—0»244»c4Gl*  0  3  0  ) 

PTLES  0 F  Pli 

( -0. 40^5+Gl ,  0,0  ) 

(-0, lo5~-Cl,-C. 755E-OI ) 


(  C  ,  c  **  C  ;  —  w 1 

(  7  .  C  ,  >  ,  u  ) 

(  0.  142c  ♦Ci.  »  0.  C  ) 

(-0.  15*>t:-Cl  ,  G.7^t-Cl) 

(  0«3p25+C1,  0 • C  ) 

(  v, «  C  »  0 .  0  ) 

(  J  •  1 4  >  c  ♦  0 1  »  J  •  0  ) 

(-0. 15rc-:i ,  0.  755--C1) 

(-0, ?Go2-Cl ,  0*  C  ) 

(  0.  :4<tE  +  Cl  ,  O.C  ) 

(-0.1p5c-C1,  0.755E-C1) 

(-0.272E-C1,  O.C  ) 

<  0,  142E  +  Q1 ,  O.C  ) 

(-0. lhhE-Cl ,  C.755E-C1) 


GAIN  Oil  *  0 . 1166+04 

011(G)  =  C, 

ZEROES  CF  -111 
(  ”0.  21  Ifc-Cl  «  C.C  ) 


POLES  CF  Qil 
(-0.244E  +  C1 ,  0.0 


GAIN  012 
012(01  * 


-12.1 

-20.5 


ZEROES  OF  012 
(  -0.  21  IE-01  *  0.0 

PCLES  CF  Cl 2 
(-0.480E  +  C1,  C.O 
< -0. 253E-C1 ,  C.O 


GAIN  C21  * 
021(0)  * 


-1*6. 

C. 


ZEROES  OF  C21 
(-0.211E-01,  C.C 

POLES  OF  021 
(  -0.  157E  +  C1  *  C.O 

GAIN  Q22  ■  8. 

022  (  0  )  *  3.  70 

ZEROES  OF  022 
(-0.211E-01,  C.O 

POLES  OF  Q22 
(-0.829E+Q1,  0.0 
(-0.277E-C1,  C.O 


(  0.0 


.  o.; 


(-0.272E-C1,  O.C 


(  0. 352E+C1 »  O.C 


(  O.o  ,  0.0 

(-0. 3C6E-01 «  0.0 


|  (  0. 685E  +  C 1  *  0.0 


CASE  NO.  IB 


(A,B>  HEANS  A+JB 


GAIN  »ii  *  -384. 

ZEROES  OF  Pil 

(-C.  52PE  +  C1  ,  C.  0  ) 

(— Q«26dc  —  Cl*  C«  j 

(  J.686E+C1, 
(C.O  * 

C.C  > 

0.0  ) 

°_LES  OF  Pil 
(-0.4046  +  C1,  C.C  ) 

(-0. 1555-01,-0. 755E-01 1 

(  0.142E+01, 
(-0. 155E-C1, 

0  .  C  ) 

0.  75 5E-  Cl  > 

GAIN  012  *  -184. 

ZEROES  OF  P 1 2 
( -0. 480E  +C1 1  0.0  ) 

(-0.253E-01,  0.0  ) 

(  0.352E+C1, 
(0.0  , 

0.0  ) 

0.0  > 

PCLES  OF  P12 
( -0.  404E  +  01  *  0.0  I 

(-0.  155E-Cl,-0.  7556-01) 

(  0.142E+C1, 
(-0.155E-01 , 

O.C  > 

0. 755E-C1 ) 

GAIN  P21  =  -20.7 

ZEROES  OF  P21 

(-0.  157E  +  C1  *  C.O  ) 

(-0.298E-C1, 

0.0  ) 

PCLES  OF  P21 
(-0.4046  +  C1,  C.O  ) 

<-0.1556-Ci,-C. 7556-01) 

(  0.142E+C1, 
(-0. 155E-C1 , 

0.0  ) 
0.755E-C1) 

GAIN  P22  *  -1.90 

ZEROES  OF  P22 
(-0.2446+01,  C.O  1 

(-0. 272E-01 ♦ 

O.C  ) 

POLES  OF  P22 
(-0.404E  +  Q1,  C.O  ) 

( -0. 155E-01 , -0. 755E-01 ) 

(  0.1426+C1, 
(-0.1556-01, 

O.C  ) 

0. 755E-C1 ) 
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1 

L, 

GAIN  oil  a  0 • 16 IE  +04 

& 

5UUI  *  C. 

&>> 

*  V,  •  « 

ZEROES  3 F  Oil 
(-0. 2265-01  ,  C.O 

) 

(0.0  , 

O.C 

a 

PCLES  OP  Oil 
(-0.244E+01  •  0.0 

) 

(-0.2726-C1, 

O.C 

* .  * 

GAIN  Q12  *  -16.7 

Q12( 0  )  »  -27.6 

m 

ZEROES  OF  01 2 
(-0.226S-C1,  C.O 

) 

•  . 

PCLES  OF  Q12 
(  -0.  48CE  +  C1  ,  0.0 
( -0.  2536-C1 »  C.O 

) 

) 

(  0. 352E+CL * 

0.0 

1  /• 

gain  021  *  -148. 

■  r*»» 

021(0)  a  c. 

,  -  .  » 

ZEROES  3F  02 1 
(-0.226E-C1,  C.O 

) 

(  0.0  , 

0.0 

kY.* 

POLES  OF  Q21 
( -0. 1576+01  *  0.0 

) 

(-0.2S8E-01. 

O.C 

SB 

m 

GAIN  02 2  •  8.00 

022(0)  »  3.70 

ZEROES  OF  022 
(  -0.  226E-C1 .  C.O 

) 

i 

PCLES  OF  022 
(-0.  82<3t  +  Cl ,  C.O 
(-0.268E-01.  C.O 

) 

) 

(  0.686E+G1 t 

o 

. 

o 

CASE  NG.  1C 


(  A  ,  8  )  MEANS  A+JB 


GAIN  PH  a  -442. 
ZEROES  OF  PH 


(-0.B2-5+C1,  C.O  ) 

(-0.2555-01,  O.C  ) 

(  0.632S+C1, 
(0.0  , 

O.C 

c.c 

PCLES  OF  PH 
( -0.  404c  +  C1 ,  C.Q  1 

(-0 .1556-C1, -0.755E-01  ) 

<  0.142E+C1, 
(-0. 155E-CI , 

O.C 

0. 755E-C1 

GAIN  P12  *  -207. 

ZERCES  OF  P12 
(-0.464E+Q1,  0.0  ) 

(-0.272E-C1,  C.  0  » 

(  0. 336E+C1 , 

(  0.  0  , 

C.C 

0.0 

PCLES  OF  P12 
(-0.4046+01,  C.O  ) 

(-0. 155E-C1 , -0. 755E-01) 

<  0. 142E+01 , 
(-0. 155E-C1 , 

0.0 

0. 755E-C1 

GAIN  P21  *  -23.5 

ZERCES  OF  P21 
(-0.1585+01,  0.0  1 

(-0.287E-C1, 

0.0 

PCLES  OF  P21 
(-0.404E+01,  0.0  I 

(-0. 155E-C1 , -0. 755E-01) 

(  0.142E+C1, 
(-0. 155E-C1 , 

O.C 

0. 755E-C1 

GAIN  P 22  -  -1.85 

ZEROES  OF  P22 
(-0.261E  +  0I,  0.0  1 

(-0.290E-01, 

0.0 

PCLES  OF  P22 
(-0.404E+01,  C.O  ) 

(-0. 155E-C1 ,-0. 755E-01 1 

(  0. 142E+C1, 
(-0. 155E-C1 , 

O.C 

0. 755E-01 

GAIN  Qli 
011(0)  * 


0.219E+G* 


C. 


2 : 3Cc 5  J-  311 
(  -C.  2  945  -  Cl  ?  C.G 

) 

( 

0.0  » 

0.0 

PCLcS  OF  Gil 
(-0.261S  +  01  ,  C.O 

) 

( 

-0.290E-C1, 

o.c 

GAIN  012  *  -19.5 

C12(0)  *  -29.9 

ZEROES  OF  G12 
(— O.  294c-Clf  C.O 

) 

POLES  OF  012 
(-0.464E  +  C1,  C.O 

) 

( 

0. 336E+C 1 t 

0.0 

( -0. 272S-C1 »  C.O 

) 

GAIN  021  *  -172. 

021(0)  *  C. 

ZEROES  OF  Q2 1 
(-C.294E-C1,  C.O 

) 

( 

0,0  t 

0.0 

POLES  OF  021 
(  -0.  158E+G1  ♦  0.0 

1 

< 

-0. 2876-C1, 

o.c 

GAIN  022  «  9.14 

022(0)  *  4.C9 

ZEBCES  OF  G22 
(-0.294E-01,  C.0 

) 

POLES  OF  022 
(-0.824E+C1 .  C.O 

) 

( 

0.682E+C1 » 

o.c 

(-o.25ae-cif  c. o 

) 

case  NC.  10 


(A,3»  MEANS  A+JB 


G-'IN  PH  »  -33S, 


Z6  =  CES  OF  Oil 
< -C.  325E  +  C1  ,  C.O  ) 

1  “0*454:” 0 1  *  0*  C  ) 

(  0. 6 36  E  +  C I . 

(  o.c  , 

O.C  ) 

O.C  ) 

?CLE  S  3F  Pil 
( -0.4046+C1,  C.O  ) 

<-0. 155E-01 ,-0. 755E-01 ) 

(  0. 142E  +  C1 . 
(-0. 155E-C1  , 

O.C  ) 

0.755E-C1) 

GAIN  P 12  »  -207. 

ZEROES  OF  P12 
(-0.A64E+C1,  0.0  ) 

C  -0.  2  72E-C1  .  0.0  I 

(  0.336E+C1, 

(  0.  C  . 

0.0  ) 

0.0  ) 

PCLES  OF  P12 
( -0,  4046  +  01  .  C.O  ) 

(-0. 155E-C1 » -0. 755E-01 ) 

(  0.142E+01, 
(-0. 155E-C1  , 

0.0  > 

0.755E-01I 

GAIN  P21  *  -18.3 

ZEROES  OF  P21 
(-0.158E+C1,  C.O  ) 

(-0.281E-C1, 

0.0  ) 

PCLES  OF  P21 
(—0.404c +01.  C.O  | 

(-0.155E-C1.-0. 755E-01 ) 

(  0. 142E+C1, 
(-0. 155E-C1, 

O.C  ) 

0. 755E-C1 ) 

GAIN  P22  *  -1.85 

ZEROES  OF  P22 

(  -0.  26  ic  +  01 .  0.0  1 

(-0.290E-01, 

O.C  ) 

ROLES  OF  P22 
(  — 0.  4046  +  C1 .  0.0  ) 

(-0.155E-C1  . -0. 755E-01 ) 

(  0.142E+01, 
(-0. 155E-C1 . 

O.C  ) 

0.755E-C1 ) 
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GAIN  Oil  a  0.171E+0* 

0  i  1  (  0  )  »  C. 

Zr^CeS  Dr  311 
(  —  0*  3  04;  —  Cl  *  C.  0 

) 

( 

0.0  , 

o 

• 

o 

PILES  Or  QH 
(  -0.261E+C1 ,  0.0 

) 

( 

-0.290E-C1 , 

o 

• 

o 

GAIN  012  ■  -15.2 

012(0  )  *  -22.6 

ZEROES  OF  Q12 
(-0.304E-G1,  C.O 

) 

PCLeS  OF  Q12 
(-0.464E+C1,  C.O 

) 

( 

0.336E+01  , 

o 

• 

o 

(  -0.  272E-C1 1  0.0 

) 

GAIN  Q21  »  -173. 

021(0)  »  C. 

ZEROES  OF  Q21 
(-0.304E-G1.  C.O 

) 

( 

0.0  , 

o 

• 

o 

POLES  OF  021 
( -0.  1 58E+G1 «  0.0 

) 

( 

-0.281E-C1, 

o 

• 

o 

GAIN  022  »  9.2 

022(0)  »  4.09 

■ 

ZEROES  OF  Q22 
(-0.304E-CI,  C.O 

I 

POLES  OF  022 
(-C.828E  +  G1,  C.O 

) 

( 

0.636E+C1 » 

o 

• 

o 

( -0. 252E-01 >  C.O 

) 
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CASS  NC.  It 

(A, 8)  MEANS 

A  +  JB 

GAIN  PH  *  -268. 

Z:»GES  OF  PH 
( -C.  82A5  +  C1 ,  C.O  ) 

(-C.247E-G1,  C.O  ) 

<  0.6336+C1, 
(0.0  , 

0.0 

0.0 

POLE S  OF  Pll 
( -0.  **046  +  C 1  ,  C.  C  ) 

( -G.  1 555-01, -0.  75  5E-01  ) 

(  0.142E+C1, 
(-0.155E-C1 , 

o.c 

0. 755E-C1 

GAIN  P12  *  -2C7. 

ZEROES  OF  012 
(-0.4o4E  +  Cl»  0.0  ) 

(-0.272E-C1,  C.O  I 

(  0.336E+C1, 
(0.0  , 

O.C 

0.0 

POLES  OF  P12 
(  -0.  404E+C1  ,  C.O  ) 

<-0.15  5t-Cl,-C.  755E-01 ) 

(  0.1425+01, 
(-0.1555-C1 , 

0.0 

0. 755E-C1 

GAIN  P21  *  -14.3 

ZEROES  OF  P21 
(-0.153E  +  C1,  0.0  ) 

(-0.2765-C1, 

0.0 

POLES  OF  P21 
(  -0.404c  ♦01,  C.O  > 

( -0. 155E-01 , -0. 755E-01) 

(  0. 142E+C1, 
(-0.155E-C1, 

0.0 

0. 755E-01 

GAIN  P22  »  -1.85 

ZEROES  OF  P22 

C  -0.  26  IE  >01 ,  C.O  I 

(-0.29QE-QI, 

0.0 

POLES  OF  P22 
(-0.4045+01,  0.0  I 

(-0.155E-CI ,-0. 755E-01 ) 

(  0.142E+01, 
(-0. 155E-01 , 

0.0 

0.755E-C1 

GAIN  Oil  * 
01 1  (  0  »  * 


0.133E+C4 


C. 


ZEROES  IF  211 
(  -0.  3 1  IE  -Cl  *  C.O 

) 

< 

0.0  , 

0.0 

PCLES  OF  211 
(-0.261E+CI,  C.O 

) 

( 

-0. 290E-C1  • 

o.c 

GAIN  Q12  »  -11.9 

Q12 (0  )  *  -17. A 

ZEROES  OF  Q12 
(-0.311E-C1,  C.C 

) 

PCLES  OF  212 
(-0.4646+01,  C.O 

1 

( 

0. 336E+C1  , 

0.0 

(-0.272E-01,  0.0 

) 

GAIN  Q21  *  -172. 

221(0)  3  C. 

ZEROES  OF  22 1 
(-0.  311E-C1 1  C.O 

) 

( 

0.0  , 

0.0 

PCLES  OF  Q21 
(-0.  158E  +  C1  •  0.0 

) 

I 

-0.276E-C1. 

0.0 

GAIN  Q22  *  9.17 

Q22  ( 0  )  «  4.09 

ZERGES  OF  222 
(-0.3UE-C1,  C.O 

) 

POLES  OF  022 
(-0.824E+01,  C.O 

) 

( 

0.6336+01, 

o.c 

(-C.247E-C1,  C.O 

1 

(  A  ,  8  )  MEANS  A  +  J8 


CASE  NO.  2A 

GAIN  PH  »  -217. 


ZEROES  OF  Pll 

(-C.74CE  +  C1,  C.O  )  (  0.6*9E+C1,  O.C 

(  —  0  •  1  o  7  c  —  Clf  C.O  )  (  0.0  v  O.Li 


PC LE  S  CF  PU 

(-0.341E+C1,  C.O  )  (  C.ISEE+Ci,  C.C 

( -0. 105c-01,-0. 579E-0 1 )  (-0.105E-01.  Q.579E-C1 


GAIN  P 12  *  -174. 


ZEROES  OF  PI 2 
<-C.*93S  +  Cl , 
(-C. 171E-C1. 

PCLES  CF  PI 2 
(-0.341E+C1 , 
(-0. 1056-01,- 

GAIN  P21  3 

ZEROES  OF  P21 

(-0.110E+01, 

PCLES  CF  P? 1 

(-0.B41E+01, 

(-0.1056-01,- 

GA IN  P22  * 

ZEROES  OF  P22 
(-0.1406+01, 

POLES  OF  P22 
(-0.341E+C1, 

(  -0. 1 05c -01  ,  - 


0.0  ) 

(  0.4C3E+C1, 

O.C 

C.  C  ) 

<0.0  , 

0.0 

C.  0  1 

(  0.155E+01, 

0.0 

C.  57SE-C1) 

(-0.105E-C1, 

0. 579E-C1 

-11.6 

0.0  ) 

(-0.215E-01, 

O 

• 

o 

C.C  ) 

(  C. 159E+C1  , 

0.0 

C. 579E-01  ) 

(-0.105E-01, 

0. 579E-01 

-3. CO 

0.0  ) 

(-0.192E-C1, 

0.0 

0.0  I 

(  0. 15SE  +  C1 , 

O.C 

0.  579E-01) 

(-0.105E-C1 , 

0.579E-01 

-116- 


I 


GAIN  Cll  *  *55. 

^IKCI  *  C. 

ZEROES  OF  Oil 

(-c.uas-ci,  c.  c 

> 

( 

J.O  , 

0,0 

) 

POLES  0 F  Qil 
(-C. 1406+01,  0,0 

) 

( 

-0.192E-C1, 

o.c 

) 

GAIN  C12  »  -7.S9 

012(0)  a  -14.4 

ZEROES  OF  Q12 
(-C.1O0E-C1,  C.C 

) 

POLES  OF  Q12 
(-0.493E  +  C1,  C.O 

) 

( 

0. 403E  +  C1  , 

0.0 

) 

(-0, 1 71S-C1 ,  0.0 

) 

GAIN  021  a  -113. 

021(0)  a  c. 

ZEROES  OF  Q ? 1 
(-C.1C86-C1,  C.C 

) 

( 

► 

o 

• 

o 

o 

• 

O 

» 

PCLES  OF  Q21 
(-0,1106  +  01,  C.O 

) 

( 

-0.215E-01, 

o.c 

) 

GAIN  022  «  6.34 

022(0)  »  4.20 

ZEROES  OF  022 
(-0.108E-01,  C.C 

) 

PCLES  OF  022 
(-C.  74SE+C1,  C.O 

) 

( 

0.6496+C1  , 

O 

• 

O 

) 

(-0.187E-C1,  0.0 

) 

M 


u 


CASE  NO.  28 


(4,8)  MEANS  A+J8 


GAIN  Pll  *  -330, 


ZEROES  OF  Pll 


i-c*7*aE*ci,  c.o  ) 

(-0,1  3  IE -Cl,  C.O  > 

(  0.64SE+C1, 
(0.0  , 

0.0 

0  .  w 

POLES  OF  Pll 
(  -0,  3<*  IE  +C1 »  C.C  ) 

(-0. 105c -01,-0. 579E-C1  ) 

(  0. 1 59E+C1 , 
(-0.105E-C1 , 

C.C 

0. 5796-C1 

GAIN  P12  =  -174. 

ZEROES  OF  PI 2 
(-0.493E  +  C1,  C.O  ) 

(-0.171E-CI,  C.O  ) 

(  0.4Q3E+C1* 

(  0.0  , 

C.C 

0.0 

POLES  OF  PI 2 

(-C.  3416+C1  ,  C.O  ) 

(-0.105E-C1,-C. 579E-01) 

(  0. 159E+01 , 
(-0.1C5E-C1 , 

0.0 

0.579E-C1 

GAIN  P21  =  -17.6 

ZEROES  OF  P21 
(-C.110E+C1,  C.C  ) 

(-0.208E-C1, 

0.0 

POLES  OF  P21 
(-0.341E  +  C1,  C.O  ) 

(-0. 105E-01 , -0. 579E-Q1  ) 

(  0. 159E+C1 » 
(-0.105E-C1, 

o.c 

0.579E-C1 

GAIN  P22  *  -3. CO 

ZEROES  OF  P22 
(-0.140E+C1,  0.0  ) 

(-0.192E-C1 , 

O.C 

POLES  OF  P22 
(-0.341E+C1,  C.O  1 

(-0.105E -01,-0. 579E-01) 

(  0.159E+CI, 
(-0. 1C5E-01 , 

C.C 

0.579E-01 

-118- 


GAIN  <311  *  691. 

Gilt  0  1  *  C. 

ZEROES  of  Oil 
(-0.  1335-01.  C.  C 

) 

(0.0  , 

0.0 

ROLES  CF  Oil 
(-0.1-OE+C1.  C.O 

) 

(-0.192E-C1  , 

C.C 

GAIN  Q12  *  -11.9 

012(0)  =  -20,9 

ZERCiS  OP  012 
(-0. 1336-C1,  C.C 

) 

POLE  S  QF  012 
(-0.493F  +  C1,  C.O 

1 

(  0.402S+C1, 

0.0 

(-0.1715-C1.  C.O 

) 

GAIN  021  =  -118. 

32110)  »  C. 

ZEROES  QF  021 
(-0.133E-CI.  C.O 

) 

(  0.0  , 

0.0 

POLES  OF  021 
(-C.110E  +  01.  C.O 

) 

(-0.208E-C1 , 

o.c 

GAIN  Q22  *  6.29 

Q22  (0  )  »  4.28 

Z5RCES  QF  022 
( -0.  133E-C1.  C.C 

) 

POLES  OF  022 
(-0.  74dE  +  Cl ,  C.C 

1 

(  0.64SE+C1, 

0.0 

(-0.181S-C1,  C.O 

) 

) 

) 

) 

> 

) 


(A, 3)  ME  AN  S  A  +  J3 


CASE  NC.  2C 

GAIN  Pli  =  -372, 

ZSPCES  OF  PH 


(-C.  75^8+01  ,  C«  G  ) 

(  0.65<*E  +  C1, 

O.C 

1-0.172E-C1,  C.G  ) 

(0.0  , 

O.C 

POLES  GF  Pli 

(-0.341E  +  GI,  C.G  ) 

(  0.15SE+C1, 

O.C 

(-0. 105E-C1,-0.575E-01 ) 

GAIN  P12  =  -198. 

(-0. 1C5E-C1 , 

0. 579E-C1 

ZEROES  GF  PI  2 

(-0.<t50E  +  Cl ,  0.0  ) 

(  0.360E+C1, 

O.C 

(  -0.  1  33E-C1 »  C.O  I 

(  O.C  , 

0.0 

PCLcS  OF  P12 

(-0.341E  +  01,  C.O  ) 

(  0.15SE+01, 

0.0 

(-0.  105E-01 » -0.  5  79E-01) 

GAIN  P21  =  -20.2 

(-0. 105E-C1  , 

0. 579E-C1 

ZEROES  CF  P21 

(-C.110E+C1,  0.0  ) 

(-0.2C0E-G1, 

O.C 

POLES  CF  P21 

(-0.341E+01,  C.C  ) 

(  0. 1 55E  +  C1  , 

0.0 

(-0. 105E-C1.-0. 5  79E-01 ) 

GAIN  P22  *  -2.5 

(-0.105E-C1. 

0. 579E-C1 

ZEROES  OF  P22 

(-0.1558  +  01,  C.O  ) 

(-0.2C3E-C1, 

0.0 

POLES  OF  P22 

(-0.341E+C1,  C.O  ) 

(  0. 15SE+C 1 , 

0  •  G 

( -0. 105E -01 ,-0. 5796-01 ) 

(-0.105E-C1, 

0. 579E-G1 

GAIN  ill  * 

giiio)  * 


0.123E+C4 


C. 


ZEROES  CF  ill 
(-C.210E-01,  C.C 

) 

( 

0.0  , 

0.0 

POLES  OF  311 
( -C. 155c«-Cl  ,  0.3 

) 

( 

-0.202E-C1, 

C.C 

GAIN  012  »  -15.5 

012(0)  *  -22.5 

ZEROES  OF  Cl 2 
(  -0.  210E-C1  *  C.C 

) 

PCLES  OF  012 

I  -C. 450E ♦Cl  »  C.O 

) 

( 

0. 36QE  +  C1 , 

o.c 

(-0.1826-01,  C.C 

1 

GAIN  C21  *  -152. 

021(C)  *  C. 

ZEROES  OF  Q21 
(-0.210E-01,  0.0 

) 

( 

0.0  , 

o.c 

POLES  OF  Q21 
( -0. 113c  +01 ,  C.O 

) 

( 

-0.200E-C1, 

O.C 

GAIN  <322  »  8.27 

C22  ( 0  )  *  «.1«5 

ZEROES  OF  <3>2 
(-C.  210F-C1,  C.C 

) 

PCLES  OF  022 
(  -0.  754E  ♦Cl ,  C.O 

> 

( 

0.654E+C1, 

0.0 

(-0.172E-C1,  0.0 

) 

-121- 


CASE  NC.  2D 


(A, 8)  V£ANS  A>J8 


GAIN  Pll  *  -277. 

ZEROES  OF  Pll 
(-0.  757E  +  C1, 


POLES  OF  Pll 
( -0.  341E  +01 » 


GAIN  P12  * 
ZEROES  OF  PI 2 


POLES  OF  P12 
(-0. 341S+C1, 


GAIN  P21  * 

ZERCES  OF  P21 
( -0# 1 09E  +01 1  ( 

PCLES  OF  P21 
(-0.341E  +  C] 
(-0.105E-C) 

GAIN  P22  * 

ZEROES  OF  P22 
(-0.155E+01,  l 

POLES  OF  P22 


C.  0  ) 

(  0.65SE+C1, 

o.c 

C.O  ) 

(0.0  . 

c.c 

C.  0  ) 

(  0.159E+C1, 

0.0 

C.  5 7<5E-0 1 ) 

( -0. 1 C5E-01 t 

0.579E-C1 

-198. 

0.0  1 

(  0.36CE+C1. 

C.C 

C.O  > 

(  0.  C  r 

0.0 

C.O  ) 

(  0.159E+01. 

0.0 

C.  5796-01) 

(-0. 105E-01  , 

C.579E-C1 

-15.2 

o 

• 

o 

(-0. 195E-C1, 

o 

• 

o 

C.  0  ) 

(  0. 159E+C1 * 

0.0 

0. 579E-01) 

(-0. 105E-C1 « 

0.579E-C1 

-2.5 

o 

• 

o 

(-0. 2036-C1. 

O.C 

0.0  ) 

I  0.159E+C1, 

O.C 

0.5796-01) 

(-0. 105E-Q1 . 

0. 579E-01 

-122- 


GAIN  QH  *  S3C. 

<311(01  *  C. 

ZEROES  OF  QH 
(-0.2215-01*  C.O 

) 

( 

0.0  f 

0.0 

) 

POLES  OF  Cll 
(-C.155E+C1,  0.0 

) 

( 

-0, 2C2E-C1 , 

C.C 

) 

GAIN  012  *  -11.7 

SI 2  (  0  )  *  -16.4 

ZEROES  OF  Cl  2 
(-0.221E-C1,  C.C 

) 

POLES  OF  C12 
(  —  0.  45  toE  ♦  Cl  t  C.O 

) 

( 

0.360E+01 , 

0.0 

) 

(-0.183E-C1,  C.O 

) 

GAIN  021  *  -153. 

C21(CI  »  C. 

ZEROES  OF  C21 
(-0.221E-C1,  C.O 

) 

( 

0.0  , 

O.C 

) 

POLES  OF  Q21 
( -0. 109E+C1 »  0.0 

1 

( 

-0.155E-01, 

o.c 

) 

GAIN  Q22  *  e.38 

022(0)  »  4.15 

ZEROES  OF  Q22 
(-0.22  IE-  Cl »  C.C 

) 

POLES  OF  C 22 
(-C. 757E+C1,  C.O 

) 

( 

0.65SE+01 * 

0.0 

) 

(-0.  lo7E-Cl,  C.C 

) 

! 


I 


(  A  f  3 )  MEANS  A+JB 


CASS  NC.  26 


GAIN  PL1 


-218. 


ZEROES  OF  Pll 
<-C.7*?5  +  Cl,  C.C 

(—  o « i  6  ■*  c — oi »  c.o 

OCLSS  OF  Pil 
(  -G.  341E  +  C1  f  C.O 
( -C. 1J56-C1,-C. 5796-01 


GAIN  P 12  * 


-196. 


ZEP06S  OF  PI 2 
(-0.4506  +  C1*  C.O 
(-C.183S-C1,  C.C 

PCLES  OF  P12 
(  -C.  3416+C1  *  o.  g  _ 
(-0.  105E-01 » -C.  5796-01 


GAIN  P21 


-11.6 


Z=RC5S  OF  P21 
(-C. 1106+01,  0.0 

PCLcS  OF  P21 
(-0.3M1E+C1  ,  0.0 
(-0. 1056-01,-0. 5796-01 


GAIN  922  * 


-2.5 


ZEROES  OF  P22 
( -0.  1556+C1 ,  0.0 

PCLES  OF  P22 
(-0.3416+C1,  C.O 
(-0.1056-01 ,-C. 5 79E-01 


(  0. 6496+ Cl ,  0.0 
(  0.0  ,  O.G 


(  0. 1 596+C1 ,  O.C 
■(-0. 105E-C1*  0. 5796-C1 


(  0.3606+01,  C.C 
(  O.C  ,  0.0 


(  0.15SE+C1,  0.0 
(-0.1C5E-C1,  C. 5796-C1 


(-0. 1926-C1,  C.O 


(  0. 1 596+C1 ,  O.C 
(-0.1056-C1,  0.5/5E-C1 


(-0.2C36-C1 *  C.C 


(  0.159E+C1,  O.C 
(-0.1056-01,  0.579E-C1 


CASE  NO.  3A 


(  A  ,  9  )  MEANS  A  +  J9 


gain  Pll  *  -567. 

ZEROES  OF  Pll 
l -0. 13+5  *02 1  0.0  ) 

(  -0.  <*42:  -01 1  0.0  ) 

(  0.1C96+G2, 
(0.0  , 

o.c 

C.O 

POLES  OF  PH 
(-0.612E+C1,  C, 0  ) 

(-0.2196-01 ,-C. 53CE-01 ) 

(  0. 146E  +  C1 , 
(-0.2196-01 , 

ol5306-01 

GAIN  P12  =  -366. 

ZEROES  OF  PI  2 
(-C.S96E  +  C1,  G.O  ) 

(-0.4276-01,  0.  0  J 

(  0.749E+01, 
(0.0  , 

0.0 

0.0 

POLES  OF  P12 
(-0.612E+C1,  C.C  ) 

(-0.2196-01,-0. 5306-01) 

(  0. 146E+01, 
(-0.2196-01, 

o.c 

0. 5306-01 

GAIN  P21  *  -31.2 

ZEROES  OF  P21 
(-0.2645  +  01,  0.0  ) 

(-0.45CE-C1  , 

0.0 

POLES  OF  P21 
(-0.6126+C1,  G.O  ) 

(-0. 2196-C1,-C. 53CE-Q1 ) 

(  0.146E+C1, 
(-0.21SE-01  , 

O.C 

0.530E-01 

GAIN  P22  *  -9.15 

ZEROES  OF  P22 
(-0.2966+01,  0.  0  ) 

(-0.4346—01  , 

0.0 

PCL6S  CF  P22 
(-0.612E  +  C1,  C.O  ) 

( -0. 2196-01, - C. 5 3CE-01) 

(  0.146E+01, 
(-0.219E-C1, 

O.C 

C. 530E-C1 

-126- 


GAIN  Oil  *  691. 

Cl  1(0)  *  C. 

ZEROES  OF  0*1 
(-0.300E-01,  C.O 

) 

( 

o 

• 

o 

*• 

O.C 

) 

PCLES  CF  Cl  1 
(-0.296E4.C1,  0.0 

) 

( 

-0.434E-C1 , 

O.C 

) 

GAIN  Q12  *  -17.2 

012(0)  =  -127. 

ZEROES  OF  Q12 
(  -0.  300E-01 ,  C.O 

) 

POLES  OF  012 
(-0.986E+C1,  C.O 

) 

( 

0.749E+01* 

o 

• 

O 

) 

(-0.427E-C1,  C.O 

) 

GAIN  021  *  -202. 

021(0)  -  0. 

ZEROES  OF  Q21 
(-0. 300E-C1,  C.O 

) 

( 

O.C  , 

O 

• 

o 

) 

POLES  OF  Q21 
(-0.264E  +  G1,  C.O 

1 

( 

-0. 450E-C1 , 

o 

• 

o 

) 

GAIN  022  -  11.2 

022 ( 0 )  *  90.1 

ZEROES  OF  022 
(  -0.  300E-01 ,  0.0 

) 

POLES  OF  Q22 
(-0.139E+02,  0.0 

) 

( 

0.  109E+02, 

O.C 

> 

(-0.492E-C1,  C.O 

) 

-127- 


CASE  NO.  3b 

(A, 3)  MEANS 

A  +  JB 

GAIN  Pll  *  -SIS. 

ZEROES  OF  Pil 

I-0.133E  +  C2,  C.O  1 

<-0.*33E-Cl»  C.  C  ) 

(  0.  109E  +  02  ♦ 
(0.0  , 

O.C 

0.0 

PCLES  CF  Pll 
(-C.612E  +  C1,  0.0  ) 

( -C. 2196-01 »-C. 53CE-01  1 

(  0. 14&E  +  C1 1 
(-0. 219E-C1 f 

O.C 

0. 530E-C1 

GAIN  P12  *  -368. 

ZEROES  OF  PI 2 
(-C.986E+C1,  0.0  ) 

(  -0.  427E-C1 «  C.O  ) 

(  0.749E+C1, 

(  0.0  , 

C.O 

0.0 

PCLES  OF  P12 
(  -0.  612E  +  01 «  C.O  ) 

( -0. 219E-01 »-C.  530E-01  ) 

(  0. 146E+01 t 
(-0.219E-01  , 

0.0 

C.530E-C1 

GAIN  P21  *  -44.4 

ZEROES  OF  P21 
(-C.265E  +  C1  ,  0.0  J 

( -0. 446E-C1 t 

0.0 

PCLES  OF  P21 

I-0.612E  +  C1,  0.0  ) 

(-0.  219E-C1  »-C.  53CE-01  ) 

(  0.146E+C1, 
(-0.219E-C1, 

O.C 

0.530E-01 

GAIN  P22  »  -9.15 

ZEROES  OF  P22 
( -0.  296E  +  01 »  0.0  ) 

(-0.4346-C1 , 

C.O 

PCLES  OF  P22 

(  -C.  612E  +  C1 »  C.G  I 

<-0.219E-0l,-0.53CE-01) 

(  0.146E+CI, 
(-0.2196-01, 

0.0 

0. 520E-C1 

-128- 


GAIN  Qll  =  966, 

Cl 1 ( 0  )  =  c. 

ZEROES  OF  Qll 
(-0.3345-01 »  C30 

) 

(0.0  , 

o.c 

PCL5S  OF  Gil 
(-G.  29fct*Clt  C.  3 

) 

(-0.434E-C1 , 

0.0 

GAIN  Q12  *  -24.1 

Q 1 2  (  0 )  *  -173. 

ZEROES  OF  Q12 
(-0.334E-C1,  C.O 

) 

POLES  OF  Q12 
( -O. 986c  +  C1 »  C.O 

) 

(  0.749E+C1, 

c.c 

(-0.427E-C1,  C.O 

) 

GAIN  Q21  »  -2CC. 

C2K0I  =*  C. 

Z5RGES  OF  021 
(-0.334E-01,  C.O 

1 

(  0.0  t 

0.0 

POLES  OF  C21 
(-0.263E+01*  C.O 

> 

(-0. 446E-C1 t 

o.c 

GAIN  Q22  *  10.8 

Q22(0)  »  40.1 

ZEROES  OF  Q2 2 
(-C.334E-C1,  C.O 

) 

POLES  OF  022 
<  -0.  133E+C2  »  C.O 

) 

(  0. 1 08E  +  C2  » 

c.c 

< -0. 43  eE-Cl  *  C.O 

) 

-129- 


CASE  NC.  3C 

(A, 3)  MEANS 

A*  J0 

GAIN  Pll  *  -976. 

ZERCES  OF  ®11 
(  -  C.  i  i  ♦  C2  *  C.  0  ) 

(  0. 1106*02, 

0.0 

(-0.432E-C1,  C.O  ) 

(  J.C  . 

0.0 

PCLhS  OF  PH 
(-0.612E  +  C1,  C.C  ) 

(  0.146E+C1. 

0.0 

(-0.219E-C1.-C. 53CE-01 ) 

(-0. 219E-C1, 

C. 530E-C1 

GAIN  P12  *  -485. 

ZEROES  OF  PI 2 
( -O. 916c ♦01 »  0.0  ) 

(  0.691E+C1. 

C.C 

(-0.436E-C1,  C.O  ) 

(  0.0  , 

0.0 

POLES  OF  P12 
(  -0.  6 12E+C1  *  C.O  1 

(  0. 1466+01. 

o.c 

(-0. 2 19c -Cl. -C. 53CE-01) 

(-0.2196-01 , 

0. 530E-C1 

GAIN  P21  =  -54.3 

ZERCES  OF  P21 
( -C.  2a*tE  +  Cl »  C.O  ) 

(-0.44CE-C1 , 

O 

• 

o 

PCLES  OF  P21 
(  -0.  6 12E  ♦  Cl  •  C.C  ) 

(  0.146E+C1, 

o.c 

(-0.219E-01.-C.53CE- 01) 

(-0.219E-C1, 

C. 530E-01 

GAIN  P22  *  -9.73 

ZERCES  JF  P22 
( -0.  3096*01 1  C.O  ) 

(-0.4426-01. 

O.C 

PCLES  OF  P22 
(  -0.  6 12E  *C1 .  C.C  ) 

(  0. 146E+C1, 

O.C 

(-0.2i9E-01.-C,  530E-01  ) 

<-0.2196-01, 

0. 520E-C1 

GAIN  Qll 
011(01  = 


0# 173E+C4 


C. 


ZEROES  2F  Gil 
(-0.-59E-01,  G*  0 

) 

( 

0.0  ♦ 

o.c 

t»CLES  OF  Gil 
(-0.309E+C1,  C.O 

) 

( 

-0. 443E—  C 1 » 

o.c 

GAIN  G12  *  -34.6 

012(0)  *  -215. 

ZEROES  OF  012 
(-C.459E-01,  C.C 

) 

POLES  3F  012 
(  -0.  916c  +01  »  0.0 

) 

( 

0.  631E  +  C1 1 

o.c 

(-C.436E-C1,  C.O 

) 

GAIN  021  =  -309. 

021(0)  *  0. 

ZEROES  OF  Q21 
( -0*  459E-C1 »  C.O 

) 

( 

0.0  , 

C.C 

POLES  QF  Q21 
(-0.264E+C1,  C.O 

) 

< 

-0.44CE-C1, 

C.C 

GAIN  022  *  17.2 

Q22(0«  =  45.5 

ZEROES  OF  022 
(  -C.  459E-01  *  C.C 

) 

POLES  OF  Q22 
(-0.134E+02.  C.O 

) 

( 

0. 110E+C2 > 

C.C 

(  -0.  4325-C1 »  C.C 

) 

CASE  NC.  2D 


U,8>  PE  AN  S  A  +  JB 


GAIN  PII  *  -72S. 

ZERCES  0 F  PII 
(-0.13^*5+02,  C.  0  ) 

( -  0. 429E- C 1 »  C.O  ) 

t  0. 1C9E+C2  , 

(  0.  c  , 

o.c 

C.C 

°CLE  S  OF  Pll 

(-0. 6125+C1 ,  C.  C  1 

(-0. 2196-01,-0.  530E-01  ) 

(  0.146E+C1, 
(-0.219E-C1, 

C.O 

0. 520E-C1 

GAIN  P12  *  -435. 

ZERCES  OF  P 1 2 
(-0.916E+C1,  0.0  1 

(-0.  436E-C1,  C.O  ) 

(  0.631E+C1, 
(0.0  , 

0.0 

0.0 

PCLES  OF  P12 
(-0.612E  +  01,  C.O  ) 

(-0.  21<?E-C1,-C.  530E-01 1 

(  0.146E+01, 
(-0. 2196-01 t 

o.c 

C.530E-01 

GAIN  P21  =  -40.3 

ZERCES  OF  P21 

(-0.  264P  +  C1 ,  0.0  ) 

(-0.437E-C1 , 

C.C 

PCLES  CF  P21 
(-0.612E  +  C1,  0.0  1 

(-0. 2 19E-C1, -C. 53C6- 01 1 

(  0.14AE+G1, 
(-0. 21SE-C 1 , 

O.C 

0. 530E-C1 

GAIN  P22  *  -9.73 

ZERCES  OF  P22 
(-C.309E+G1,  0.0  ) 

(-0.443E-01, 

0.0 

PCLcS  OF  P22 
(-0.612E  +  C1,  C.C  1 

(-0.219c-01f-C.530E  — 01  ) 

(  0. 146E+C1, 
(-0.2 19E-C1 » 

0.0 

C. 520E-C1 

GAIN  Qll  a 
Q 1 1  (  0  )  = 


0.1286+04 


C. 


ZEROES  OF  Qll 
(-0.4796-C1,  C.O 

) 

( 

o 

• 

o 

O.C 

POLES  OF  Qll 
(-0.309E+C1,  0.0 

) 

( 

—0,44  26—  Cl, 

O.C 

GAIN  Q12  *  -25.7 

Q12I0)  =  -156. 

ZEROES  OF  Q12 
(  -0,  479E-C1 ,  C.O 

) 

POLES  OF  Q12 
(  -0.  916c  +  Cl .  0.0 

) 

( 

0.6816+01 , 

0.0 

(  -0.  426E-C1 »  C.C 

) 

GAIN  Q21  *  -209. 

Q2 1  (  0  )  =  C. 

ZEROES  OF  Q21 
(-0.4796-01,  0.0 

) 

( 

0.0  , 

O.C 

PCLES  OF  Q21 
(-0.264E+01,  C.O 

) 

( 

-0.437E-C1, 

o 

• 

o 

GAIN  Q22  »  17.1 

Q22(  0  )  *  45.5 

ZEROES  OF  Q22 
(-0.4796-01,  C.C 

) 

POLES  OF  Q22 
( -0. 1 34fc  +  C2  *  C.Q 
(-0.429E-C1,  C.C 


)  (  0. 1096+C2 »  O.C 

» 


GAIN  311  =  G . 

311(0)  =  C. 

1C2E+C4 

ZEROES  jF  3 4  1 

(-0.491E-C1,  0.0 

) 

(  0.0 

POLES  OF  311 

(-0.  3G9E  +  G1 ,  C.O 

) 

(-0.443E-C1 

GAIN  312  =  -20.5 

Q  i  2  C  0  >  =  -128. 

zeroes  of  312 

(-C.491E-C1,  C.O  ) 

POLES  OF  Q12 

(-0.916E+CI,  C.O  )  (  0.681E+CI 

(-C.436E-C1,  C.C  ) 

GAIN  Q2i  =  -3C4. 

321(0)  =*  C. 

ZEROES  OF  Q21 

(-0.491E-C1,  0.0  )  (  0.0 

POLES  OF  321 

(-0.265E+G1,  C.O  )  (-0.435E-C1 

GAIN  322  *  16.3 

322(0)  »  45.5 

ZEROES  OF  322 

<-C. 491E-C1,  C.O  ) 


OGLES  JF  322 
(-Q.133c+02,  C.C 
(  -C.  <*27E-C1  *  C.O 


)  (  0. 1QSE+C2 

) 


-135- 


CASE  NO.  4A 

(A, 3)  MEANS 

A+JB 

GAIN  Pll  *  -S3. 2 

ZERCES  OF  Pll 

(-0. 4455+ Cl,  C.C  ) 

(-0. 149E-C1.  C.C  ) 

(  0  •  3946+0 1  . 

(  0  •  C  t 

0.0 

O.C 

PCIES  OF  Pll 
(-0.214E  +  C1.  C.C  ) 

(-0.959S -02,-0. 731E-01  ) 

(  0. 124E  +  01  . 
(-0.959E-C2. 

C.O 

0. 781E-C1 

GAIN  P12  *  -52.9 

ZERCES  OF  P 1 2 
(-0.234E+01,  C.O  ) 

(-0.9696-C2,  C.O  ) 

<  0.191E+C1, 
(0.0  , 

O.C 

0.0 

POLES  Or  P12 
(-0.214E  +  C1,  C.O  ) 

(-0.  959E-C2.-0.  781E-01 ) 

(  0.124E+C1, 
(-0. 958E-C2 , 

C.C 

0.7  81E-01 

GAIN  P21  =  -4.86 

ZEROES  OF  P21 

(-0.  5  sfcE  +C  J  »  C.O  ) 

(-0, 255E-C1 , 

O.C 

POLES  OF  P'2 1 
(-0.214e-*-Cl,  C.O  ) 

(-0.959F-C2.-0. 7 81E-01 ) 

(  0.1246+C1, 
(-0.955E-C2, 

O.C 

0.7  61E-C1 

GAIN  P22  =  -0.3C1 

ZERGES  OF  P22 

(-0.  132c+Cl .  0.0  » 

<-0. 1416-C1, 

O.C 

POLES  OF  P22 
(-0.214E  +  C1,  C.C  » 

(-0.959E-C2.-0. 781E-01 ) 

(  0.124E+C1, 
(-0. 95GE-C2 . 

C.C 

0. 761E-C1 

-136- 


GAIN  Oil  *  761. 

Cl  1 ( 0  )  =  0. 

ZEROES  OF  Oil 
(-0.  576E-C2,  0.0 

) 

(0.0  , 

O.C 

POLES  OF  Cll 
( -0. 132E  +  C1 T  C. 0 

) 

(-0. i 41E-C1  . 

O.C 

GAIN  Q12  *  -■4. 33 

012 1  0  I  =  -2.6 

ZEROES  OF  012 
(-0.576E-C2.  C.O 

) 

POLES  OF  Q12 
( -0. 2 3^c +01 ♦  C.O 

) 

(  0.191E+C1, 

O.C 

(-C.969E-02.  C.O 

) 

GAIN  Q21  *  -47,1 

021(0)  =  0. 

ZEROES  OF  02 1 
(-0.576E-C2,  C.O 

) 

(  O.C  t 

O.C 

POLES  OF  Q21 
(-0.556E  +  C0,  0.0 

) 

(-0.255S-C1, 

0.0 

GAIN  022  *  2.46 

Q22 ( 0 )  *  0.337 

ZEROES  OF  022 
(-0.576E-C2,  C.O 

) 

POLES  OF  022 
( -0.  4*»5E  •♦•01 1  C.O 

) 

(  0.294E+C1, 

O.C 

( -0.  1<*9E-Ci  ,  C.O 

) 

i 


I 


-137- 


CASE  NC.  A3 

(A, 8)  *EANS 

A  + JB 

GAIN  Pli  =  -121. 

ZEROES  OF  Pil 
(-G.^96+01  ,  0,0  ) 

<-C. 1-05-01,  C.O  ) 

(  0.398E+C1, 
(0.0  , 

0  .  c 

0.0 

PCLES  OF  Pil 
J-C.Zl^E+Cl,  C.  0  ) 

<-C.95  9E-C2.-C.  7  81E-01) 

(  0.12AE+C1, 
(-0. S5SE-C2  , 

O.C 

0.781E-Q1 

GAIN  P12  *  -52.9 

ZSRCES  OF  P 1 2 
(-0.23AE+C1.  C.O  ) 

(-C.969E-02.  C.O  ) 

<  0.1916+C1, 

(  O.C  , 

O.C 

O.C 

PCLES  OF  P12 
(-0.21^6  +  01,  C.O  I 

(-0.959E-02.-0. 781E-Q1 ) 

(  0.12AE+01, 
(-0.  9596-C2  , 

0.0 

0.7S1E-01 

G1IN  P21  *  -6.46 

Z5»06S  OF  P21 

(-0. 555E+C0,  0.0  ) 

(-0. 2A6E-C1, 

C.O 

PCLES  OF  P21 
(-0.214E+01,  C.O  i 

( -0. 959E-C2 » -0. 781E-01 ) 

(  0.12AE+C1, 
(-0. 959E-C2 , 

C.O 

C.761E-01) 

GAIN  P22  »  -0.301 

ZEROES  OF  P22 
(-0.1326+01,  C.O  ) 

(-0.1A1E-Q1 , 

0.0  ) 

POLES  OF  P22 
|-0.21*E+C1,  C.O  ) 

(-C.959E -02,-0.7816-01) 

(  0.12A6+C1, 
(-0.959E-02 , 

O.C  ) 

0.781E-01 ) 

GAIN  Gil  *  0.1C1E+04 

C1K0)  «  C. 

ZEROES  -JF  Gil 

( -C. 6*66-02 t  G,  0  ) 


PCLES  OF  qj.1 
(-0.  132E+C1,  G.  0 


GAIN  Q12 
012(0)  » 


-5.77 
-A.  49 


ZEROES  OF  Q12 
(-C.646E-C2*  C.  C 

PCLcS  OF  Cl 2 
<-C.  234E  +  Q1  ,  c.C 
(-0.969E-C2,  C.C 


GAIN  C21 
321(G)  = 


-47.2 

C. 


ZEROES  OF  Q21 
(-0.646E-C2,  C.O 

POLES  OF  321 
(-0.555E+CC,  0.0 


GAIN  022 
Q22 ( 0  )  * 


2. Si 

0.337 


ZEROES  OF  Q22 
(-0.646E-C2,  C.C 

PCLES  OF  022 
(-C.448E  +  C1,  0.0 
(-0.  1406-01,  C.O 


(  G.o  ,  Q.c 

(-0.141E-C1,  Q.C 


(  0.191E+01,  O.C 


<  O.C  r  Q.C 

(-0.246E-C1,  O.C 


J  (  0.398E+C1 »  0.0 


CASS  NO.  AC 


(A, 8)  WEANS  A+JB 


GAIN  PI!  3  -136. 

ZEROES  GF  Pll 


(-C.45IS  +  C1,  C.  0  > 

(-0.127S-C1,  C.O  ) 

(  0.4C1E+C1, 
(0.0  . 

o.c 

o.c 

PCLES  CF  Pll 
(-0.214E+C1,  C.  C  ) 

(-0.959t-C2.-C. 7  31 6-01  ) 

(  0.  124E  +  C1 » 
(-0. 959E-C2 . 

C.C 

0. 7  filc-Cl 

GAIN  P12  3  -66.5 

ZERCES  GF  P12 
(  -C,  2206  +  Gl  .  C.O  ) 

(-C.  1286-Cl.f  C.O  ) 

(  0.176E+C1. 
(0.0  . 

O.C 

0.0 

PCLES  GF  PI2 
(-0.214E  +  01,  C.O  ) 

(-0.9596-02,-0. 781E-01 ) 

(  0.124E+C1, 
(-0.  959E-C2  » 

0.0 

C.761E-01 

GAIN  P21  3  -7.38 

ZEROES  OF  P21 
(-0.555E+C0,  0.0  ) 

(-0.232E-C1, 

O.C 

PCLES  CF  P21 
(-0.214E+CI,  C.C  ) 

(-0.959E-02. -0. 781E-01) 

(  0. 124E+C1, 
(-0.959E-C2. 

0.0 

0.7E1E-CI 

GAIN  P22  3  -0.242 

ZERCES  OF  P22 
(-0.179E+G1,  0.0  1 

(-0.160E-CI, 

0.0 

POLES  GF  P22 
( -0.  2 14c  +01 ,  C.O  ) 

(-Q.959E-C2.-0. 781E-01) 

(  0.124E+C1, 
(-0.95SE-C2  t 

O.C 

0.781E-01 

G- IN  Gil  » 
011(0)  * 


0.189E+04 


Z£°CE  S  OF  qil 
(  -0.  12 3c  -Cl »  C.  0 

PCLES  OF  Oil 
(-0.179E  +  C1,  0,0 


GAIN  012 
G12 ( 0  )  = 


-6.88 

-2,99 


ZEROES  OF  G12 
( -0.  128E-C1 »  C.  0 

PCLES  CF  Cl 2 
(-0.  220c ♦Cl ,  C.  0 
( -0. 128E— Cl »  C,C 


-62.0 


gain  C21  * 
021(0)  = 


ZEROES  OF  02 1 
(-0.  1286-01,  C.O 

PCLES  OF  G21 
(-0.555E+CC,  C.O 

GAIN  Q22  »  3.26 

022 ( 0 )  *  0.419 

ZEROES  OF  Q22 
(-0.1236-01,  C.O 

ROLES  OF  022 
(-0.4515+01,  0.0 
(-0.127E-C1,  C.O 


)  (0.0  ,  0 . ( 

'  (-0.16C6-01,  0.< 

) 

j  (  0. 1766+01,  O.C 

)  (  O.C  ,  0.0 

»  (-0.232E-01,  0.0 

) 

j  (  0.4C1E+C1,  O.C 


-141- 


CASE  NC.  40 


( A  ,  0  )  MEANS  A  +  JB 


GAIN  PH  *  -1C1. 

ZEROES  OF  Pll 
( -0.452=+Cl ,  0.0 
(-0. 117E-C1 ,  C.  C 

POLES  CF  Pll 


GAIN  P 12  * 

ZERCES  OF  PI  2 
( -O.  2206  +  C1 .  ( 
(-0.1286-01,  ( 

POLES  JF  PI  2 
(-0.2146+C1,  ( 
(-C.  9596-02, -( 

GAIN  P21  = 

ZEROES  OF  P21 
(-C.  5556  +  CC,  ( 

POLES  OF  P21 
(-G.2146+C1,  < 
( -0.  9596-C2  »“( 

GAIN  P22  * 

ZEROES  OF  P22 
(-0.179E  +  C1,  < 

POLES  OF  P22 


0.  0  ) 

(  0.4C3E+C1. 

o.c 

C.  C  ) 

(0.0  . 

0.0 

C.  0  ) 

(  0.124E+C1* 

o.c 

0.  7 8 1E-0 1 ) 

(-0.  9596-C2, 

0. 761E-01 

-66.5 

C.O  ) 

(  0.116E+C1, 

O.C 

C.O  ) 

(0.0  . 

0.0 

0.0  ) 

(  0.124E+C1, 

O.C 

0. 781E-01  ) 

(-0.9596-02, 

0.7  61E-01 

-5.  51 

C.  J  ) 

( -0.222E-01  , 

0.0 

0.0  ) 

(  0.124E+C1, 

C.C 

0.  781E-01) 

(-0.  959E-C2  , 

0. 7S1E-01 

-0.242 

0.  0  1 

(-o.i6ce-eit 

O.C 

C.  C  ) 

(  0. 1246+01  . 

0  •  c 

0,  7816-CU 

(-0.959E-C2, 

0#  761E-01 

-142- 


GA  IN  Qll  *  0.141E+04 

G 1 1  (  0  >  *  C. 

ZEROES  OF  Gil 

(-0.132E-01,  0.0  )  (  0.0  ,  0.C  ) 

POLES  OF  Gil 

( -0. 1 79E ♦Cl t  C.O  I  (-0.16JE-01.  O.Q  ) 

GAIN  C12  *  -5.14 

012(C)  *  -2.75 

ZEROES  OF  012 
(-0.132E-C1,  C.O  ) 

POLES  OF  012 

( -0. 220c >01 >  C.O  )  (  0.176E+C1,  C.C  ) 

(-0.128E-C1,  C.C  ) 

GAIN  021  *  -62.0 

G 2 1  (  0  )  *  C. 

ZEROES  OF  021 

( -0. 132E-01 »  C.O  )  (  0.0  ,  O.C  ) 

POLES  OF  C21 

( -0.  555c  +00 »  G.O  J  (-0.222S-C1,  O.C  ) 

GAIN  022  *  3.38 

022(0)  =»  0.<tlS 

ZEROES  OF  022 
( -0. 1 325 -01 .  C.O  ) 

POLES  JF  Q22 

(-0.452S+01,  0.0  )  (  0.403E+C1,  O.C  > 

(-0.1172-C1,  C.O  ) 


CASE  NO.  4? 

( A , S  )  MEANS 

A*  J8 

GAIN  Pll  *  -79.3 

ZEROES  OF  PU 
t,“C#^ol.r  +  LL  9  C«  0  ) 

(  0.402E  +  0U 

O.C 

(-0.UJE-C1,  C.  C  ) 

(  o.c  , 

0.0 

PCLES  OF  PU 

(-0.  21<MTfCl  ,  C.C  > 

(  0.1245+C1, 

O.C 

(-0. 9596— C2 » -C. 7  816-01 1 

(-0. 959E-C2 , 

0.7  61E-C1 

GAIN  P12  -  -66.5 

ZEROES  uF  P12 
(  -0.  220E+C1 «  C.O  ) 

(  0.176E+C1, 

O.C 

(  -0.  1 2  dc-Cl  *  C.O  I 

(  0.0  , 

0.0 

PCLES  OF  P 1 2 

(-0. 214E+C1 .  C.O  > 

(  0.124E+C1, 

o.c 

(  -C.  959E-C2.-0.  781E-01  ) 

(-0. 959E-02  » 

0. 781E-01 

GAIN  P21  =  -4.30 

ZEROES  OF  P21 

(-0.  556E+CC,  C.O  ) 

(-0. 215E-C1 , 

O.C 

PCLE  S  OF  P21 

(-0. 214E  +  01  ,  C.C  ) 

(  0. 124E+C1 * 

o.c 

( -0.9596-02,-0. 781E-01 I 

(-0. 959E-C2 . 

0.781E-C1 

GAIN  P22  *  -0.242 

ZEROES  OF  P22 

(-0.1 796  +  01 »  C.O  ) 

(-0.160E-C1, 

O.C 

PCLES  OF  P22 
(-0.214E+C1,  C.C  ) 

(  0.124E+C1, 

O.C 

(-0.959t-Q2,-G.7aiE-01) 

( -0. 959E-C2 . 

0.7E1E-C1 

GAIN  Qll  *  0.110E+C4 

cm  a »  =  c. 

ZEROES  of  Til 

(-0. 136E-01 ,  C.O  )  (  0.0  ,  O.C 

POLES  OF  Til 

(-0. 179E+C1,  0.0  )  (-0.16CE-01*  O.C 

GAIN  Q12  *  -4.C2 

012(0)  =  -2.01 

ZEROES  OF  Q12 

(-0.  136E-C1,  C.  C  ) 

POLES  OF  Q12 

(-C.220E+C1.  C.O  )  (  0.116E+C1,  0.0 

(-0.12SE-G1,  C.C  ) 

GAIN  Q21  *  -62.0 

021(0)  =  C. 

ZEROES  OF  021 

( -0. 136c -01 >  C.O  )  (  0.0  >  0.0 

POLES  OF  Q21 

(-0. 556E+C0t  0.0  )  (-0. 2 15E-C1 »  O.C 

GAIN  Q22  *  3.37 

022 ( 0  )  *  0.419 

ZEROES  OF  022 
(-C.136E-C1,  C.O  ) 


POLES  OF  Q22 
(-0.451E  +  C1,  0.0 
(  -0.  1 10E-01 »  O.C 


)  (  0.4C2E+C1'  C.C 

) 


(A, 3)  MEANS  A*JB 


CASE  NC.  5A 


CAIN  Pll  =  -134, 


ZERCES  ZF  Pll 
(-C.715E+C1,  C.  C  ) 

(-0,1345-01,  C.O  ) 

PCLES  CF  Pll 

(-0. 3056*01,  C.O  ) 

(-0.  3556-02,-0.  536E-01  ) 

GAIN  P12  *  -145. 

ZEROES  OF  P12 
(-0.5C6E  +  C1,  0.0  ) 

(  -0.  120E-C1 ,  C.O  ) 

PCLES  OF  PI 2 
(-0.3056  +  C1,  C.O  ) 

(-0. 8555-C2, -0.  536E-01 ) 

GAIN  P21  -  -10.0 

ZEROES  JF  P21 
( -0. 9156  *C  0 ,  0.0  ) 

PCLES  OF  P21 
(-0.3056*01,  C.O  ) 

<-0.855E-C2,-0. 536E-01 » 


GAIN  P22  *  -3.37 


ZERCES  OF  P22 
( -0. 1 065*C1 , 


0.0  ) 


PCLES  OF  P22 
(-0.3056*01,  0.0  l 

(-0. 8555-02,-0. 5366-01 ) 


(  0.627E+C1,  O.C 

(  O.C  ,  C.C 


(  0.1456+C1,  C.C 
(-0.855E-C2,  0. 5266-C1 


(  0.424E*C1,  0.0 

(  0.0  ,  0.0 


l  0. 145E+01  ,  0.0 
(-0.8556-C2,  0.536E-01 


( -0, 163E-0 1 ,  0.0 


(  0.  145E  +  C1 ,  C.O 
(-0. 855E-C2 ,  0. 536E-01 


(-0. 1446-01,  0.0 


(  0.1456+C1,  O.C 
(-0. 855E-C2 ,  0.536E-01 
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1 

i 

i 

lw 

9 

i 

i 


GAIN  Cll  =  259. 

ClllCI  =  c. 

ZEROES  OF  Cll 

(-C.  279E-C2,  C.O  )  (  0.0  ,  0.0  ) 

POLES  OF  Cll 

(-0. 106E+01 ,  C.O  )  (-0.l44E-01t  O.C  ) 

GAIN  C12  *  -5.86 

<312(0!  =  -11.4 

ZEROES  OF  Q12 
(-0.279E-C2.  C.O  ) 

POLES  OF  Q12 

(-0.5C6E+C1,  C.O  )  (  0. 424E+C1 »  O.C  ) 

(  -0.  120E-C1  *  C.C  ) 

GAIN  Q21  *  -87.4 

C2HCI  *  C. 

ZEROES  OF  Q21 

(-0.279E-C2,  C.O  >  (  0.0  »  0.0  ) 

POLES  OF  C21 

(-0. 915E+CC,  C.O  I  (-0.163E-01,  0.0  ) 

GAIN  C22  *  4.76 

Q22  (  0  )  *  3.S6 

ZEROES  OF  Q2 2 
(-Q.279E-C2.  C.O  ! 

POLES  OF  Q22 

( -0. 713E+01 i  0.0  )  (  0. 627E+C1 *  C.C  ) 

(-0.134E-C1,  C.O  > 


t 


1 


CASc  NC.  5? 


(  A  ,  8  )  MEANS  A  +  J8 


GAIN  P 11  =  -272. 

ZEROES  Jr  Pll 
(  -0 .  7  L  +01  »  C.O  ) 

(-C.129E-C1,  C.C  ) 

(  0. 633E+C1 t 

(0.0  t 

- 

V-  •  U 

0.0 

PCLES  CF  Pll 
(-0.3J5E+C1,  C.O  ) 

(-C.  6  5  5E-C2  t-C.  5366-01 ) 

(  0.145E+C1. 
(-0. 655E-C2 , 

o.c 

0. 536E-01 

GAIN  P 12  *  -149. 

ZERCES  OF  PI 2 

(-0.  506E  +  C1,  C.O  ) 

(-0.120E-C1,  C.O  ) 

(  0.424E+C1, 
(0.0  , 

C.C 

0.0 

PCLES  OF  PI  2 

(-C.  205E+C1,  C.O  ) 

(-0. 6  5  5E-C2.-0. 5366-01  ) 

(  0.145E+01, 
(-0. 855E-C2 ♦ 

0.0 

0.536E-01 

C-AIN  P 21  =  -15.1 

ZEPCES  OF  P21 
( -C.  9 1 4c  +0  0  ,  C.O  ) 

(-0. 1576-01, 

0.0 

PCLES  OF  P21 
(  -0.  305E  +  01 «  C.O  ) 

(-0. 655E-C2  ,-C. 536E-01) 

(  0.1456+Clt 
(-0. 6556-02, 

o.c 

0.536E-01 

GAIN  P22  *  -3.37 

ZEROES  OF  P22 
(-0.106E+01,  C.O  J 

(-0.1446-01, 

C.C 

PCLES  OF  P22 
(-G.3C5E  +  CI,  C.O  ) 

(-0.8556-02,-0. 536E-01) 

(  0.145E+C1, 
(-0. 8556-02 , 

0.0 

0. 536E-01 

Cl  1 ( 0  ) 


c 


i  r?ces  o f  on 

( —  0 . fc 2 6 c — C 2 »  C. 0 

) 

(0.0  » 

c.c 

POLES  OF  Cl  1 
(-0.106E+C1,  C.O 

) 

(-0. 144E-C1, 

C.C 

GAIN  £12  =  -8.99 

012(0)  *  -16.5 

ZEPGES  CF  Cl 2 
(—  0.626E— 02  »  C.  C 

) 

PCLES  OF  Q12 
(-0.506E+C1,  0.0 

) 

(  0.424E+C1, 

o 

• 

O 

(  -C.  12CS-C1 »  C.O 

) 

GAIN  021  a  -88.7 

021(0)  =  C. 

ZEROES  OF  Q21 
( -0. 626E-C2  t  C.O 

) 

(  0.0  t 

C.C 

PCLES  OF  021 
(-0.914E+CC,  C.O 

) 

(-0. 157E-C1, 

O 

• 

O 

GAIN  022  *  4.93 

022(0)  *  3.  96 

ZEROES  OF  022 
(  -0.  626E-C2  f  C.O 

) 

PCLES  OF  022 
(-0.719E+01,  0.0 

) 

(  0.633E+C1 » 

o 

• 

o 

(-0.129E-C1,  C.O 

) 

L  4  5C 


NO. 


(A, 8)  MEANS  A+J8 


GAIN  PH  = 


-314. 


2 1 5 C E  S  JF  PH 
• -C. 713E+C1 ,  C.C  ) 

f  '«»  lZit  - Cl,  C  »  G  ) 

POLES  OF  PH 
(-0.305E+C1,  C.O  ) 

( -C. 855e-C2 , -0. 536E-01 ) 


GAIN  P 12 


-  156. 


ZERCES  OF  P12 
( -C.  465E  +  C1  *  0.0  ) 

(  -0.  132E-C1 1  C.C  I 

PCLES  OF  P12 
(-C.3056  +  C1,  C.O  ) 

(-C.  855E-C2  » -C. 536E-01 ) 


GAIN  P2l  r 


-17.5 


ZEROES  OF  P21 

(-G. 914E+C0,  C.O  | 

PCLES  OF  P21 

(-0 .3056  +  01*  C.O  ) 

(-0.855E-C2.-0. 5366-01 ) 


GAIN  P22  * 


-3.52 


ZERCES  OF  P22 
(-0.113E+C1,  0.0  ) 

POLES  OF  P22 
(-0.305E  +  C1*  C.O  ) 

(-0.8  5  55-02,-0.  5366-01 1 


(  0.634E+C1,  C.C 
(0.0  ,  C.O 


(  0.145E+C1,  C.C 
(-0.855E-C2,  C.536E-01 


(  0.387E+01,  O.C 
(  O.C  ,  C.O 


(  0.145E+C1,  o.o 
(-0. 855E-C2  *  C.526E-C1 


(-0.150E-C1,  C.C 

<  0.145E+C1,  C.C 
(-0.855E-C2,  0.536E-C1 


(-0.155E-C1,  0.0 


(  0.145E+C1,  C.C 
(-0.855E-C2,  0 • 536E-C1 


POLES  OF  Oil 
(-0*1136  +  01,  C.O 


GAIN  Q12 
212(0)  * 


-11.  <3 
-17.6 


ZEROES  OF  212 
( -0.  178E-C1,  C.C 


POLES  OF  Q12 
(-0.  <*696  +  01  ,  C.O 
(-0.1326-01,  C.O 


GAIN  Q21  * 
221(0)  * 


-132. 


ZEROES  OF  Q2 1 
(-0.1786-01,  C.O 


POLES  OF  Q21 
(-0.9146+CQ,  C.O 


GAIN  Q22  * 
222(0)  » 


7.38 
4.  73 


ZEROES  OF  222 
(-0.178E-01,  C.O 


POLES  OF  Q22 
(-0.7196+01,  0.0 
-0.  122E-C1,  C.O 


)  (-0.1S5E-C1,  0. ( 

) 

j  (  0. 3876+G1 ,  Q.C 
\ 

)  (  0.0  ,  Q.C 

)  (-0. 15C6-C1 ,  C.O 

) 

J  (  0.6346+C1,  O.C 


CASE  NO.  50 


(A, 8)  M£4\S  A  +  JB 


GAIN  «li  =  -23C. 

ZEROES  CF  =>■! 


(  -C.  72-E+C1  *  C.C  ) 

( -0.  1136  — Cl*  C.C  ) 

(  0.635E  +  01 , 

(0.0  , 

0.0 

0.0 

PCLtS  OF  PU 
(-0.3056+01,  C.O  ) 

(-0. 855E-C2.-0.  536E-01  ) 

(  0. 145E+C1 » 
(-0.855E-C2* 

C.O 

0. 5366-01 

GAIN  P12  *  -186. 

ZEROES  OF  P12 
(-0.^*696  +  01,  0.0  ) 

(-0.  132E-C1,  C.O  ) 

(  0.3876+01, 
(0.0  , 

0.0 

0.0 

POLES  OF  P12 
(  -0.  3  056  +C1  ,  C.O  I 

( -0. 8  55E-C2  »"Q.  5366-01 1 

(  0. 145E+C1 » 
(-0. 655E-C2 , 

0.0 

0. 5366-01 

GAIN  P21  *  -13.0 

ZERGES  OF  P21 
(-0.913E+CC,  C.O  ) 

(-0. 1466-01, 

0.0 

POLES  OF  P21 
(-0.305E+C1,  C.O  ) 

(-0.8556-02 ,-0.5366-01) 

(  0. 145E+C1, 
(—0, 8556-02 ♦ 

C.C 

0.536E-C1 

GAIN  P22  »  -3.52 

ZERGES  OF  P22 

(-0.  113E  +  C1  »  0.0  1 

(-0. 1556-C1, 

C.O 

POLES  OF  P22 
(-0.305E+C1,  C.O  > 

(-0. 855E-C2 ,-0. 536E-01 ) 

<  0. 145E+C1 , 
(-0. 8556-02, 

0.0 

0.5366-01 
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GAIN  Oil  * 

311(01  =  C. 

ZEROES  OF  311 
(-0.193E-C1,  C. 


POLE S  IF  Qll 
(-0. 113E+01 ,  C.C 


494. 


GAIN  012 
012(0)  * 


-8.  68 
-12.7 


ZEROES  OF  Q12 
l -0. 193E-01 r  C.C 

PCLES  OF  .012  „ 

(-C.469E  +  01*  C.O 
(-0.132E-01,  C.O 


-133. 


GAIN  021  * 
021(0)  = 


ZEROES  OF  02 1  „ 

(-C.  193E-C1,  C.O 

PCLES  OF  Q21 
(-C.913E+C0,  C.O 

GAIN  Q22  *  7.54 

022(0)  *  4.73 

ZEROES  OF  022 
(-0.193E-01,  C.C 


PCLES  OF  022 
( — C. 724E+C1, 
(-0.118E-Q1, 


C.O 


C.O 


1  (0.0  ,  C.O 

1  (-0.155E-CI.  C.C 

) 

)  (  0. 387E+C1 »  0.0 

) 

)  (  0.0  r  0.0 

)  (-0. 146E-C1  i  C.O 

I 

)  (  0. 638E+C1 »  O.C 

1 


I 

) 

1 

) 

) 


CASc  NC.  5E 


(A, 8)  »*EANS  A  +  J8 


GAIN  Pi:  =  -136. 

Z  E  R  C  6  S  'JF  01  1 

(-0.7J9E+01,  C.O  )  {  0.62*c+Clt  0.0 

(  —  Oi  i  Cl  r  0j  C  )  (0.0  »  C  •  C 

POLES  OF  Oil 

(-0.305S+C1,  0,0  )  (  0,  145E  +  C1 .  O.C 

(-0.855c-02.-0. 536E-01)  (-0.S55E-02,  0.5366-01 

GAIN  P12  *  -15c. 

ZEROES  OF  PI 2 

(  -0.  4695 ♦ Cl .  0.0  >  (  0. 387E  +  C1  .  0.0 

(-0.132E-C1,  C.C  >  (  0.0  ,  0.0 

POLES  OF  PI 2 

( -0.  3 05E ♦  Cl »  C.O  >  (  0.145E+C1,  0.0 

(-0.855E-C2.-0.536E-01)  (-0.855E-C2»  0.526E-C1 

GAIN  P21  *  -10.0 

ZEROES  OF  P21 

(-0.91SE+C0,  0.0  J  (-0.144E-C1,  O.C 

POLES  OF  P21 

(-0.305E+C1,  0.0  J  (  0.145E+C1,  O.C 

( -0. 855E-02  » -0. 5  36E-01 )  (-0.855E-C2,  0.536E-C1 

GAIN  P22  *  -3.52 

ZEROES  OF  P22 

( -0. 1 13c ♦01 .  0.0  )  (-0.1556-Clt  O.C 

POLES  OF  P22 

( -0. 305c  *01 »  0.0  I  (  0. 145E+C1 1  O.C 

(-0.855E-02.-Q. 536E-01)  (-0.855E-C2.  0.536E-C1 
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'i 


GAIN  Qll  =  371. 

auto  *  c. 

ZEROES  C F  Cll 
( -0. 2G5E-CI «  C. 0 

) 

(  0.0  , 

3.0 

POLES  OF  Qll 
(-C.112S+C1,  C,C 

) 

(-0. 1556-C1 , 

C.C 

GAIN  Q12  *  -6.68 

01 2  I  0  J  *  -9.61 

ZEROES  QF  C12 
(-0.  2056-01,  C.C 

J 

POLES  OF  C12 
(-C.469E  +  C1,  G.O 

) 

(  0 • 387E  +  C1  , 

O 

• 

O 

(-0.  1326-01,  C.O 

) 

GAIN  C21  *  -13C. 

G  2 1  (  0  )  *  C. 

ZEROES  OF  Q21 
(  -0.  205E-C1 ,  0.0 

I 

(  0.0  , 

o 

• 

o 

POLES  OF  Q21 
(-0.913E+00,  C.O 

) 

(-0. 1446-C1, 

o 

• 

o 

GAIN  Q22  *  7.C3 

Q22  ( 0  l  *  4.73 

ZEROES  OF  Q22 
(-C.2C5E-C1,  C.C 

) 

POLES  OF  C22 
(-C.7C9E  +  C1,  C.O 

) 

(  0.624E+C1, 

o 

• 

o 

(-0.1166-01,  C.C 

) 

|Q*il 


design 


DEGREES 


Chapter  4 


design  of  r..v:  translation  <  •  ,)  mode 


-  .  I  Problem  s ta tenen t 

In  thij  mode  chi  plane  inputs  are  listed  as  o'  ,  :  ,  i 

ora 

(canard,  rudder,  aileron),  and  its  outputs  are:  side  velocity 
v  =  / A  dt  ,  yaw  {i  =  /r  dt  ,  roll  ®  (Fig.  4.1).  The  objective  is  to 
achieve  a  commanded  side  velocity  v  with  small  yaw  0  and  small 


roll  o  angles.  For  a  command  of  v^  =  60  ft/sec,  we  shall  aim  for 


,  o  ,l 

S  1  tor  i v  i 


and  5  2°  for 


Hence 


21  (60)57.3 


=  -70  db 


b31  =  2b21  =  -64  db 


(4. la,b) 


As  for  b^  ,  we  take  in  the  meantime  the  model 


2  i-l 


Tu(s) 


r  /\  +  £_W1  +  s  V  1 

L  V1  .3  y  v1  4 ;  j 


(4.2) 


whose  Bode  plot  is  sketched  in  Fig.  4.2.  We  do  not  consider  the  lower 


bound  a^(w)  for  iTn(jw)  I  at  this  point,  because  in  the  previous 
designs  the  corresponding  loop  transmission  was  so  large  that  the 
(T^^  in  Chapter  2,  in  Chapter  3)  was  almost  invariant  over  its 

important  w  range. 


4 . 2  Plant  equations 

While  the  plant  system  here  is  essentially  the  same  as  in  Chapter  2, 
there  the  variables  of  interest  were  sideslip  3  ,  roll  t  ,  side 
force  Ay  whereas  here  they  are  side  velocity,  yaw  and  roll,  so  we  get 
different  plant  equation.  In  (2.1),  yaw  was  not  present  while  side  slip 
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was  in  che  output  vector  £  .  But  the  procedure  for  finding  the  new 
plant  matrix  is  the  same  as  in  Chapter  2,  so  is  not  repeated  here. 

He  re . 


=  l’ 


The  elements  Q..  of  Q=P^=[1/Q..] 
1J  1J 

Chapter  2  are  given  in  Table  4.1  below, 
form 


for  the  ten  cases  liated  in 
in  which  the  entries  are  of  the 


-13  2 

Q..  =  D..s  +  C..s  +  A..  +  B. .  (4.4a) 

ij  ij  LJ  ijs  1J 


except  for 


Qi3 


=  s [D . .s' 
ij 


C.  .s 
1J 


A  .s 


B.  .] 
ij 


(4b) 
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TAHLc  4,1  tL-’^NTS  Q I .)  (SEE  SEC.  4.2) 


CAS-  1 


»'  ■-* 

i 

J 

r 

10 

in 

s 

-t 

-  ■>  J 

/•> 

1  'J 

o.  j  j 

Jo  .  ■* ’Jn*  J  . 

0  .  -  r  =  -  J  - 

-  •  :  -  -v  6  *  0  1 

i 

1  ->  .  *T  T  "  ~ 

wi  J 

.O  ^  +  )  .* 

•_ 

■  , .) 

,  0 

: . .  4  -  _  7 1 

c,  7  "7  v, 

i‘.  - 

;  j 

*7  .  _  • 

J  1  _  -  _  „ 

—  ,  1  3  ' . 1  J  j 

—  .  -  -  ^  _ 

-\  *' 

i.. 

-  -  -  j 

-  J  -  <_  •  * 

2  J  -  .  . 

w  ,  -  - 

rf.  •' 

•»  ; 

3  ^  •  v  i  - 

*4..‘ 

■5 

J  ■>  .5 

-,77---  y. 

■: ,  1 7  v.  :■  -  o  i 

—  t  7  *  i  J 

»  . 

3 
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4.3  Design  of  L, 


The  inequality  to  be  satisfied  is 


^33  ^33  I 

bn  |  |  +  b2i  |  q^-  | 

|i  +  L3|  5 - t  b3 

b31 

For  u  e  [0.1,5]  ,  b,,  I  ^  j  »  b  |  ^33  I 


'11  I  Q- 


21  |  Q. 


(4.5) 


and  >  60  db  ,  so 


we  need 


bn  Q,, 

|l  +  L,|  2  >  60  db 

3  b31  Q31 


(4.6) 


Since  L, 


IQooG,|  >>  1  ,  |l  +  L  |  =  | L  |  (4.6)  becomes 


|g3I  5 


b31  ^31  min 


(4.7) 


The  resulting  bounds  on  | ( juj)  |  are  shown  in  Fig.  4.3,  and 
L30  =  Q330G3  in  Fig‘  4,4‘ 


For  a)  >  8  in  (4.5),  B^wO  db,  so  the  only  requirement  is  stability 
of  (1+L3)  and  a  reasonable  margin  | L 3 / ( H-L 3 )  |  <  u  ;  u  =*  5  db  was 
chosen.  The  results  are  shown  in  Fig.  4.4.  Some  templates  of  Q33  are 
shown  in  Fig.  4.5.  At  low  cj  values,  a  right  half-plane  f  le  of  Q33 
near  the  origin  and  a  zero  at  s  -  0  (recall  (4b))  dominate,  so  we  use 
G3  ~  -8(s+0.1)/s  for  w  <  5  ,  which  handles  both  w  <  . 1  and  the 
requirement  |g|  >5  for  ui  €  [0.1,5]  .  The  final  =  G3  Q330 

chosen,  is  shown  in  Fig.  4.4,  with 


-171- 


a 


-0. 8( l+10s) 


G3(s)  = 


V1  *7*2  ) 


0  +  to)  (‘  +  ) 


(4.8) 


4.4  Design  of 


The  inequality  to  be  satisfied  (see  Chapter  1),  is: 


bH 


<22 


|1  +  L2|  5 


‘21 


+  b31  I  "22 

31  I  Q 
I  ^2  3 


=  B, 


'21 


(4.9) 


For  u>  €  [.1,5]  ,  Oj 

as  in  Eqs .  [4.6,7]  ,  we  need 

b. 


Q22 

>>b 

|  Q22 

’  11 

q21 

31 

Q23 

|g2U 


li 


b21|Q2ll 


>  30  db  . 


,  and  B2  >  70  db  ,  so 


(4.10) 


For  a)  small  ,  b 


31 


Q22 

>>  b 

Q22 

Q23 

11 

Q2] 

because  of  the  zero  of  Q23  at  the  origin,  so  G 2  is  assigned  a  pole 
at  the  origin  and  its  low  frequency  form  is 


jOlOs  +  1) 


K  *  30  db 


(4.11) 


see  Fig.  4.6.  For  u>  >  15  ,  ~  0  so  the  only  demands  are  of 

stability  and  suitable  margin  for  | L ^ / ( 1  +  L2)  j  *  The  template  of 
Q22  determine  the  resulting  bounds  on  L.,q,  see  Figs.  4.7,  4.8. 


L~,Q(j~)  is  sketched  in  Fig.  4.8,  and 


G .(s'  = 


-3 . 1 ( lOs+L) ( . ls+1) 


4.5  Design  of  G 


As  in  Chapter  2,  the  exact  equation  is  used  for  the  design  of 
the  final  loop  involving  .  The  exact  expression  for  T  ^  is 
manipulated  into  the  form 


'ii'11 


’  L1  "  ^11G1  •  Q11  = 


Yl2(l'^3)+Y13(l+L2)~r 

(i+l2)(i+l3)-y23 


(4 . 13a-f ) 


Q. .  Q.  . 

JJ 

Q.  .  Q.  , 

ij  Ji 


r  »  o  0  0  ^  - i _  +  - i _  1 

11  22  33  ^  Q12Q23Q31  Q21Q32Q13  J 


The  design  of  F^  ,  G^  is  therefore  precisely  a  single-loop  design, 

* 

for  which  Ref.  [3]  is  highly  appropriate.  Templates  of  Q  are  obtained 

(Fig,  4.9).  These  are  used,  as  described  in  [3],  to  obtain  bounds  on 
*  -  *  * 

L1q  ,  such  that  A|L^/(1+L^)[  does  not  exceed  the  allowed  range.  The 
bounds  shown  in  Fig.  10  were  used  here.  The  upper  bound  does  not  correspond 
to  that  used  in  Fig.  4.2.  This  may  affect  the  T01  ,  T  results.  If  it 
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does  do  so  too  much  adversely,  we  can  redesign  .  In  any 

resulting  bounds  on  L  are  shown  in  Fig.  11,  as  well  as  L 

r  <  ,  0 . 01 ( 5s+i 

1  s(.ls+l) 

which  involves  very  little  bandwitdth  for  G^(jw)  .  (It  might 
idea  to  make  larger.)  Two  choices' of  F^  were  tried 


case 

10 


,  tiie 
with 


14; 


be  a  good 


F 


la 


0.3 
(s+. 3) 


lb 


0.5 
(s+. 5) 


(4 . 15a, b) 


4 . 6  Simulation  of  Design 

The  two  designs  (F^,  F^)  were  simulated  on  the  digital  computer. 

The  results  are  shown  in  Figs.  4.12,  4.13  respectively  and  are  seen  to  be 

very  good.  The  commanded  v  was  adjusted  so  that  in  most  cases,  one  of 

the  control  surfaces  (it  was  6^  the  canard)  saturated.  The  yaw  and 

roll  values  obtain.  are  very  small. 

Obviously,  tne  v  responses  with  F,,  are  faster  than  those  for 

ib 

Ffa  •  However,  the  canard  then  saturates  for  some  cases  (those  with 

some  overshoot)  at  smaller  commanded  v  values.  One  could  perhaps 

eliminate  the  overshoot  for  these  cases,  by  careful  examination  of 

these  cases,  and  this  would  be  worth  doing  in  a  design  meant  for  actual 

implementation.  If  the  overshoot  in  some  of  these  6  responses  could  be 

eliminated,  the  limiting  factor  in  speeding  up  the  response  of  v  would 

be  the  saturation  values  of  the  5.  rather  that  those  of  some  of  the  .5.  . 

i  i 

In  any  case,  the  results  achieved  nicely  satisfy  the  specifications  of 
Sec.  4.1. 
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Fig.  4.12  (cont.)  Simulation  results,  F  *  0.3/(s+.3) 
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Fig.  4.14  Ccont.)  Simulation  results  with  digital  compensation,  T  -  0.02  seconds,  F 


Expanded  time  scale  In  last  two  figures 


Chapter  5 


DESIGN  BASED  ON  NONLINEAR  MULTIPLE 
INPUT-OUTPUT  PLANT  MODEL 

5. 1  Nonlinear  MIO  Design  Technique 

In  the  preceeding  chapters  a  linear  MIO  (multiple  input-output)  model 
was  used  for  the  plant.  Such  a  model  is  only  an  approximation  of  the  more 
accurate  nonlinear  model  (Eqs.  5. 1-5. 5  below).  Our  design  technique  is  able 
to  cope  precisely  with  a  large  class  of  highly  uncertain  nonlinear  MIO  plants 
It  is  enqphasized  that  the  design  technique  is  exact.  The  rigorous  develop¬ 
ment  of  the  technique  with  two  detailed  numerical  examples,  is  given  in  [  4  ] 
For  the  sake  of  conqpleteness,  a  qualitative  description  of  the  technique  is 
here  included. 

Let  y^  *  wx^  be  the  nonlinear  relation  between  a  n-vector  x^  of  plant 
inputs  and  a  n-vector  y^  of  plant  outputs.  Due  to  uncertainty  w  e  {jJ ,  a  set. 
Let  there  be  n  such  (x^,y J  vector-pairs,  forming  two  nxn  matrices 
X  -(x^, —  x^)  ,  Y  «  (y1,...,yn)  with  the  x^vector  the  i*"*1  column  of  X,  etc, 
and  y_.  »  wx^;  j  *  l,...,n.  It  is  assumed  that  the  n  x_.  are  an  independent 
set,  as  are  the  n  y^.  Let  the  matrix  P  of  transfer  functions  be  defined  by 
Y  *  PX,  where  Y  is  the  matrix  of  Laplace  transforms  of  the  elements  of  Y, 
etc.  The  P  matrix  is  denoted  as  the  LTIE  (linear  time  invariant  equi valent) 
of  the  mio  nonlinear  plant  w,  with  respect  to  the  n-vector  Y  set  of  outputs. 

A  different  X,Y  set  gives  in  general  a  different  P,  when  w  is  nonlinear. 
Repeat  the  process  over  a  set  V  of  n-vectors,  to  obtain  a  set  of  LTIE 
plants.  Repeat  over  the  elements  of  W  *  {w} ,  to  obtain  an  overall  set 
P  *  (P)  of  LTIE  plants.  This  set  P  is  the  LTIE  of  the  set  W,  with  respect 
of  the  set  of  n-vectors  V  -  {Y}  (or  the  set  X  =  {X}). 


The  following  analogy  may  be  helpful  in  comprehending  the  LTIE  idea. 


Imagine  a  barrel  V  =  (v;  of  all  the  desired  n-output  vectors  the  nonlinear 

plant  is  supposed  to  deliver  over  its  liftr.  Also  there  is  a  barrel 

W  =  {w}  of  all  possible  nxn  nonlinear  plants,  due  to  the  parameter  uncertainty 

We  create  by  the  above  technique,  a  barrel  P  =  {P}  of  nxn  linear  time 

invariant  plants,  which  has  therefore  the  following  property.  Suppose  one 

picks  any  w^  from  the  W  barrel  and  any  y  vector  from  the  /  barrel.  Then 

there  is  a  n-vector  x. .  of  inputs,  such  that  y.  =  w.x. ..  Then  in  the 

ij  31  13 

P  barrel,  there  is  guaranteed  to  be  a  P. .,  such  that  Y.  =  P. .x. ..  Thus,  the 

13  3  13  13 

P  barrel  is  equivalent  to  the  W  barrel,  with  respect  to  the  V  barrel. 

The  vital  point  is  that  the  designer  can  forget  about  (Jj .  If  he 
makes  a  proper  linear  time  invariant  design  for  the  P  barrel,  then  this 
same  design  is  guaranteed  to  work  for  the  ft)  barrel.  By  a  proper  design  for 
the  P  barrel,  it  is  meant  that  the  actual  plant  output  is  a  member  of  the 
desired  output  set  V,  no  matter  which  P  in  P  is  used  as  the  plant.  Thus, 
the  nxn  nonlinear  uncertainty  design  problem  has  been  converted  to  a  nxn  mio 
linear  time  invariant  uncertainty  problem.  It  is  important  to  include  in 
y  a  good  sampling  of  all  the  desired  nonlinear  plant  output  vectors, 
because  the  linear  design  is  guaranteed  to  work  only  for  the  P  set  which 
is  generated  in  the  above  manner  by  the  W,  /  sets,  i.e.  P  is  a  function 
of  y  as  well  as  of  W. 

5.2  The  Nonlinear  Plant 

Because  of  the  nonlinear  coupling  between  them,  the  2x2  vertical 
translation  (a  )  mode  of  Chapter  3  and  the  3x3  lateral  translation  (82) 
mode  of  Ch.  4  were  combined  into  a  single  5x5  nonlinear  plant: 


Control  inputs:  Elevons  5  ,  flaps  5^,  canards  5^,  rudder  5^ 
Plant  outputs:  Vertical  velocity  v  ,  pitch  angle  3,  lateral 
yaw  angle  i< .  roll  angle  i. 


ailerons  6 
velocity  v 


Coxjnand  inputs  and  performance  colr-ranc-s:  T:i-.re  are  conar.and-c  simul  tar.-.-eus  iy 

vertical  velocity  of  25  ft. /sec.  and  lateral  velocity  of  60  ft. /sec.,  with 
tolerances  on  these  outpus  as  in  Chs.  3,4.  The  other  tolerances  are: 
pitch  angle  0  _<  0.165,  yaw  < J>  <_  1.0,  roll  <p  <_  2 ,  all  in  degrees. 

Nonlinear  Equation  of  Motion 


I  -  I  I  I  M 

.  y  z  xz  .  xz  x 

P  =  -  q*r  +  - —  r  +  - —  pq  +  — 

X  XXX 
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: —  r 
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(5. la-c) 
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•  .  x 

u  =  v*  r  -  w*  q  +  — 
m 


F 

•  y 

v  =  w*p  -  u-r  +  — 1 *- 


(5. 2a-c) 
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V  =  Total  velocity 


aQ  =  Trim  AOA  (0  =  aQ) 


u„  =  V  cos  a  .  =  V  sin  a 

0  0  0  0 


u  =  u  dt  +  Uq,  v  =  v  dt,  w  =  w  dt  + 


(5.6a-e) 


^  -1  w  .  -1  v 

a  =  tan  —  ,  8  =  sxn  — 
u  V 


Strictly  the  aerodynamic  coefficients  C  are  functions  of  a  and  6. 
However,  here  they  are  so  restricted  in  magnitude  by  the  specifications  that- 
they  are  effectively  constants.  The  non linearities  are  due  to  the 
multiplication,  etc,  of  state  variables  in  the  above  equations  (5. 1-5. 5). 

5.3  The  Nonlinear  Plant  Input-Output  Sets 


Formally ,  the  LTIE  P  set  is  obtained  by  picking  any  y  vector, 
solving  for  x  =  w  ^y  and  repeating  over  ^  and  W.  This  is  one  way,  which 
was  used  in  [4  ] .  Another  way  is  to  simulate  the  nonlinear  plant  and  try 
to  drive  it  in  such  a  way  as  to  produce  the  kind  of  y  vectors  required  by 
the  specifications.  This  is  normally  very  difficult,  but  in  this  problem 
we  already  have  the  linear  designs  of  Chs.  3,4.  So  the  same  f,g  elements 
were  used  as  compensations  around  the  5x5  nonlinear  plant,  the  entire  system 
simulated  on  the  digital  computer  and  subjected  to  the  inputs  listed  in 
the  specifications.  It  was  hoped  that  the  nonlinear  5x5  plant  outputs 
would  not  differ  very  greatly  from  those  achieved  separately  in  the  linear 
2x2  and  3x3  designs.  If  they  did  so  differ,  one  could  experiment  with  other 
command  inputs  and  hopefully  find  some  which  would  result  in  the  desired 


output  vectors.  It  was  an  experiment  worth  trying. 
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The  results  achieved  were  very  interesting,  because  they  were 

closely  similar  to  those  achieved  in  the  two  separate  linear  designs.  This 

is  seen  by  comparing  the  two.  Figs.  5.1  a,b  (nonlinear  results)  should  be 

compared  with  6f.  5^  curves  respectively,  in  Fig.  3.13  (p.  165), 

Figs.  5.2a,b  with  the  6,  v  curves  of  Fig.  3.13.  Each  of  these  gives  the 

z 

results  for  the  five  cases  listed  in  Ch.  3.  It  is  seen  that  case  by  case 

the  curves  are  very  similar.  Figs.  5.3a-c  present  the  nonlinear  results 

for  v, ((),(>  (cases  1-5),  to  be  compared  with  the  linear  results  in  Fig.  4.12, 

d.  183  for  the  same  variables;  Figs.  5.4a-c  (nonlinear,  cases  6-10)  with 

the  linear  results  of  Fig.  4.12,  p,  184;  Figs.  5.5a-b  for  6,5.6 

era 

(nonlinear,  cases  1-5)  with  Fig.  4.12  (p.  185)  linear  results;  Figs. 

5.6a-c  (nonlinear,  cases  2-10)  with  Fig.  4.12  (p.  186)  linear  results. 

It  is  seen  that  case  by  case,  the  results  are  very  similar  for  both  the 
control  inputs  and  the  plant  outputs.  This  proves  that  the  LTIE  of  the  5x5 
nonlinear  plant  (oombined  lateral  and  longitudinal)  for  the  above  cases,  and 
for  this  class  of  plant  outputs  is  essentially  that  used  in  the  two 
separate ,  decoupled  linear  problems  of  Chapters  3,4.  Hence  straightforward 
application  of  our  technique  to  the  resulting  5x5  LTIE  problem  would 
undoubtedly  result  in  the  same  f,g  compensation  functions  obtained  in  the 
two  separate  2x2,  3x3  designs.  In  fact  our  nonlinear  simulation  with  these 
compensations  has  already  been  done.  There  is  no  need  for  any  more  work. 
Figs.  5. 1-5. 6  are  indeed  the  simulation  results. 

The  above  results  suggest  that  at  least  in  mio  flight  control 
problems,  it  is  worth  doing  first  a  design  based  on  linearized  models,  before 
embarking  on  one  based  on  the  more  accurate  nonlinear  model.  The  compensa¬ 
tions  resulting  from  the  linear  design  may  then  be  tried  for  the  nonlinear 


plant.  Comparison  of  the  results  with  the  linear  system,  indicates  how 
serious  are  the  nonlinearities.  In  the  above  example,  the  results  would 
undoubtedly  be  significantly  different  if  the  1,3  ranges  were  much  larger, 


such  that  the  C  are  then  in  the  strongly  nonlinear  ranges  of  tnese 


variables. 
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CHAPTER  6 

FUSELAGE  YAW  POINTING  MODE  (S  ) 

In  this  mode  the  control  surfaces  and  the  plant  outputs  are  trie  same 

as  in  Chapter  2  (Direct  Side  Force  Mode),  so  the  same  plant  equations  apply, 
except  for  the  numbering  (see  6.2). 

6.1  Specifications 
Part  I  2x2  System 

Ninety  degrees  roll  ($)  is  commanded  and  the  sideslip  (B)  is  to  be  main¬ 
tained  as  small  as  possible  with  +2°  as  extremes.  Time  to  roll  90°  is  to  be  as 
small  as  possible  (0.75  sec  if  possible)  with  little  or  no  overshoot.  After 
steady  state  has  been  achieved,  a  sinusoidal  roll  is  commanded  as  in  fig.  6.1-a 
tentatively  up  to  13.5°  amplitude  and  period  bounds  .5  to  2.5  sec  (.8  ir  to 
4  it  rps).  The  roll  output  is  to  be  within  the  tolerances  given  in  fig.  6.1-b 
where  $c  is  the  commanded  roll,  i.e.  error  <  15®.  The  two  control  inputs  are 
the  ailerons  and  the  rudder;  the  canard  input  is  not  used. 

Part  II  3x3  System 

Repeat  part  I  with  sideslip  step  command  added  after  steady  state  has 
been  achieved,  as  in  fig.  6.1-c.  The  resulting  side-force  is  to  be  less  them 
0.2  g  and  the  canard  is  now  used.  The  sideslip  output  is  to  be  within  the 
tolerances  given  in  fig.  6. Id. 

6.1.1  Compatibility  of  Specifications  with 

Inherent  Capabilities _ 

The  control  input  numbering  is  (1)  6  -ailerons,  (2)  5  -rudder, 

3  IT 

(3)  5c~canards.  The  outputs  are  (1)  $-roll  angle,  (2)  S-sideslips, 

(3)  A  -side force, 

y 
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2 

Part  1:  For  the  roll  step  response ,  the  transfer  model  T. .  =  [  ( s / .>  )  “  + 

-  11  n 

1.6  (s/v>  )  +■  1]  was  used  with  damping  factor  ^  =  .  8  to  achieve  small 
n 

overshoot :  *  =  e  in  order  to  arhi-vv-  final  valu-  rrcssinc?  at  t  '  .  ~~ . 

and  u  =  8  was  tentatively  tried  as  an  utcer  tound  Oi.  T,,  !  ju) .  Consider  the 
n  "  11 

resulting  6  due  to  sinusoidal  commanded  roll  ($  ).  From  6  =  P  ^Y , 

d  S  S 

Y  =  ($,6,A  )  1  ,  P  1  =  [P.  .  ]  ,  6  =  P,  Ji  +  P,  8  +  P,  ,A  ,  and  it  is  found  that 

y  I]  a  11  12  13  y 

at  oi  =  12.5<13.5°  commanded  sinusoid  amplitude,  and  the  assigned  tolerances 

of  6.1,  6  »  P  $ ,  with  |  <M  j  12. 5)  |  =  (13.  5° )  |t  ( j  12. 5)  .  This  calculation 
3  X  X  XX 

gives  a  range  of  27°  minimum  (case  4)  to  190°  maximum  (case  2)  for  5  ,  both 
exceeding  the  ailerons  limit  value.  A  similar  calculation  of  5  gives  a 

cl 

range  of  364  deg/sec  (Case  4)  to  2500°/sec  (Case  2),  both  far  exceeding  the 
ailerons  rate  limit. 

The  inverse  calculation  can  be  made:  at  13.5°  amplitude  command,  given 

|<5almax  *  20°,  I lmav  *  56°/sec,  and  =  8  for  the  second  order  T^ls) 

model,  it  is  found  that  rate  saturation  would  occur  at  id  =  3  (Case  3)  and 

position  saturation  at  id  =  9.  For  case  4,  rate  limit  as  low  as  id  =  2,  position 

at  id  =  3.  The  conclusion  is  that  id  *  8  is  beyond  the  system  capability  for 

n 

13.5°  sinusoidal  input.  In  the  meantime  we  use  id  =  5,8  as  the  bounds  on 

n 

T  <j«),  but  only  for  step  input.  Since  the  sinusoidal  roll  command  input 
generally  is  from  another  source,  it  is  possible  to  use  a  separate  prefilter 
for  the  sinusoidal  signals. 

6.1.2  Specifications  on  T 

Those  on  were  given  above.  For  T^,  roll  due  to  sideslip  step 
command  £  2°  dominates,  giving  |t^|  <  2/5  =  0.4.  For  T 22  (B  response  due  to 
8  command) ,  the  following  transfer  function  models  were  found  suitable. 


Upper  bound  b^:  {  (g-*-l)  (-j-j+1)  (80+1'^ 

Lower  bound  a22:  {  (~+l)  2  <^+l>  <^+l)  >_1 

T  =  3/4>  <  2/90  -  :V0J2  dominates  diving  b.  ,  .  T,  =  A  / J _  =  ( 0 . 13) 

21  c-  2132  - 

(32) (57. 3)/5  gives  b,_  =  55  ft/sec/deg,  while  T,_  =  A  /$  =  (.05) (32) (57.31/13.5 

$2.  31  y  c 

gives  b31  =6.8  ft/sec/deg. 

6.2  Design 

The  detailed  design  steps  are  the  same  as  in  the  other  modes ,  and  there 
is  therefore  no  need  to  repeat  them  here.  They  are  somewhat  more  complicated 
because  (a)  the  same  compensations  are  wanted  for  both  the  2x2  system  and  the  3x3 
system,  and  (b)  both  step  and  sinusoidal  inputs  must  be  considered.  Due  to 
the  former,  the  most  severe  bounds  on  G^  and  G (emanating  from  the  2x2  and 
3x3  problems)  must  be  used.  As  for  the  second,  it  has  already  been  noted  that 
the  sinusoidal  response  tolerances  were  not  compatible  with  the  system  satura¬ 
tion  levels  and  the  step  response  tolerances ,  so  the  approach  was  to  design 
the  loops  (G^ ,  i=l-3)  to  satisfy  the  step  response  tolerances,  and  then  find 
separate  prefilters  for  the  sinusoidal  commands  such  that  saturation  resulted 
only  for  an  extreme  flight  condition.  The  precise  design  details  are  not  given, 
because  the  design  technique  is  the  same  as  has  already  been  greatly  detailed 
in  previous  chapters,  especially  Chapter  2.  However,  the  templates ,  bounds 
etc  and  the  compensation  transfer  functions  are  presented. 

6.2.1  Design  graphs 

The  frequency  domain  bounds  on  J T  f  ju>)  |  ,  [ T22  ( j<u)  |  are  shown  in  Figs. 
6.2a,b.  Templates  of  (2x2)  are  shown  in  Fig.  6.3a,  Q  (3x3)  in  Fig.  6.3b, 

Q22  (3x3)  in  Fig.  6.3c,  (3x3)  in  Fig.  6.3d.  Nichols  chart  bounds  on 

and  sketches  of  L.  (joj)  are  shown  in  Fig.  6.4  a-c  for  i  =  1  to  3  (case  9  is 


nominal) . 
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6.2.2  Compensation  functions 

The  compensation  functions  are; 
g  =  -0.5  ( . 8s+l) / < 100 , 0 -ft) 
g2  =  14  ( . 25s+l)/ (120,0. 5 ) 

g3  =  0.14  (150,1200,250)/ (30) (1200,0.5) (500) 
f  (step)  =  1/(6, 0.8) (1200,0.5) (500) 

f  (sinusoidal)  =  (4. 34)  (4. 34)/ (4. 34, . 707)  (4. 34 ,. 707) (12 . 5) 
f  (step)  =  1. 4/(5) (5) (13) (80) 
f  (sinusoid)  =  f  ]_  (sinusoid) 

2 

(A,B)  indicates  a  complex  pair  [(s/u^l  +  (2J s/uO  +  1]  with  A  =  ,  B  =C; 

(E)  indicates  a  factor  [(s/E)+l].  The  first  number  is  the  value  of  the 
transfer  function  at  s  =  0 

6.2.3  Simulation  Results 

In  Figs.  6.5-7,  the  system  is  2x2  for  the  roll  command  #  of  90°  at 
t=0.  At  t=1.8  sec.  there  is  a  sideslip  8C  of  5°,  and  the  system  is  3x3. 

The  roll  responses  are  shown  in  Fig.  6.5a,  sideslip  in  6.5b,  sideforce  in 
6.5c.  The  results  are  very  good,  except  that  in  6.5a,  the  response  time  is 
not  the  desired  0.75  sec.  But  this  is  because  the  system  is  physically 
incapable  of  such  fast  response.  This  is  seen  in  Fig.  6.6a,  6.7a  where  the 
ailerons  limit  in  rate  in  most  cases  and  in  amplitude  in  two  extreme  cases. 

The  rates  are  shown  in  Figs.  6.7a-c. 

In  Figs.  6.8-6.10,  the  system  is  3x3.  A  step  roll  command  is  given  at 
t*0  and  a  sinusoidal  roll  command  with  amplitude  of  15°  and  frequency  12.5  rps 
is  given  at  t=1.8  sec.  The  three  responses  are  shown  in  Figs.  6.8a-c.  The 
control  surface  time  histories  are  shown  in  Figs.  6.9a-c  and  their  rates  in 
6. lOa-c.  From  Fig.  6.8a,  6.9a  it  is  seen  that  the  best  possible  response  is 
being  achieved,  because  of  the  ailerons  amplitude  and  rate  limiting. 
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Figs.  6.11-6.13  are  similar  to  Figs.  6.8-10,  i.e.  the  system  is  3x3, 

there  is  a  step  roll  command  of  90°  at  t=0  and  a  sinusoidal  command  at  t=1.8, 

of  15°  amplitude.  But  now  the  frequency  is  the  lower  bound  one  of  4  rps , 

with  much  better  roll  sinusoidal  response  as  seen  in  Fig.  6  Lla  (note  that 

saturations  of  6  ,5  occur  in  the  extreme  cases  -  Figs.  6  12a.  6  13a)  The 
a  a 

actuator  responses  are  shown  in  Figs.  6.12a-c.  their  rates  ■_  i  Figs.  6.13a-c. 

Figs.  6.14-6.16  are  the  responses  in  the  3x3  system  dt-sign  due  to  a 
2.5°  sideslip  command  at  t=0 ,  followed  by  a  sideslip  sinusoidal 

command  at  t=l  of  1.5°  (see  Fig.  6.1c)  amplitude,  frequency  15  rps  ; 

6.14a-c:  0 ,  3,  A^;  6.14d-f:  ailerons,  rudder,  canard;  6.14g-i  :  rates  of 

these  control  surfaces.  In  the  next  group.  Figs.  6.l5a-i,  the  command 
is  again  sideslip,  step  of  2.5°  at  t=0,  sinusoid  at  t=l  of  1.5°  amplitude, 
but  now  the  frequency  is  4  rps. 

In  all  of  the  above  the  interactive  specifications  are  excellently 
satisfied.  The  responses  to  commands  which  do  not  satisfy  the  desired  tolerances 
are  due  to  the  physical  inability  to  do  so,  i.e.  the  response  specifications 
are  not  reasonable  in  view  of  the  ailerons  amplitude  and  rate  limits. 


L  ANGLE(D 


NARDS  RATE(DD 


CASES  1  -  5 


Figure  6.8b,c  Sideslip  and  lateral  acceleration,  3x3  (*  =90°, 

c 

t=l.  8,.  sinusoid  ui  =  12.5  rps,  amplitude  15°). 
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CASES  1  -  5 
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Figure  6.14a  Boll 

Figure  6.14  a-i  Responses  due  to  slides lip  command  8c:  2.5°  step  at  t»0, 
sinusoid  at  t»1.0,  amplitude  1.58,  frequency  12.5  rps. 
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Figure  6.14b  Sideslip 
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Figure  6.14c  Acceleration 
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Figure  6.14d  Ailerons 
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Figure  6.15a  Roll 


Figure  6.15a-i  Responses  due  to  sideslip  command  ;  2.5°  step  at  t=0, 
sinusoid  at  t=1.0,  amplitude  1.5°,  frequency  4  rps. 
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CASES  1  -  5 
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Figure  6.15g  Aileron  rate 
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Figure  6.15h  Rudder  rate 
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DISCUSSION  AND  CONCLUS IONS 


The  outstanding  feature  of  the  design  procedure  presented  in  Chapter  1, 
is  its  conversion  of  the  multiple  input-output  (MIO)  design  problem  with 
uncertainty,  into  a  number  of  single  input-output  (SIO)  single  loop  design 
problems  with  uncertainty  and  disturbances.  For  a  large  problem  class,  it 
is  guaranteed  that  the  designs  which  satisfy  the  SIO  problems,  are  also 
satisfactory  for  the  original  MIO  problems.  To  use  this  method,  it  is  very 
helpful  to  assign  apriori  performance  tolerances  on  the  system,  to  be  achieved 
despite  the  parameter  uncertainty.  These  can  be  tentative  tolerances,  changed 
in  the  course  of  the  design  because  of  the  transparency  of  the  design  technique, 
which  clearly  reveals  the  trade-off  between  the  severity  of  the  tolerances  and 
the  price  paid  in  loop  gain  and  bandwidth. 

Heretofore  this  design  procedure  was  successfully  applied  to  academically 
contrived  prcuxems  [1-3,6].  They  all  had  very  large  parameter  uncertainty, 
and  the  most  complex  involved  a  3x3  plant  [2]  with  a  range  of  plant  instability, 
and  some  nonmin imum- phase  elements.  In  its  closed-loop  performance  specifica¬ 
tions,  some  off-diagonal  channels  were  specified  to  be  "basically  noninteracting", 
while  some  others  were  to  have  significant  interaction. 

The  aircraft  problems  treated  in  this  report  are  the  first  practical  real- 
life  ones  to  be  handled  (at  least  by  the  authors)  by  this  new  technique,  which 
revealed  its  flexibility  in  several  ways.  A  high-frequency  condition  (Sec.  1.3.1) 
existed  for  this  technique,  which  while  not  absolutely  necessary,  was  still 
extremely  useful  to  simplify  the  design.  This  high-frequency  condition  was  not 
satisfied  in  Chapter  2.  A  method  was  found  for  overcoming  this  obstacle  and 
which,  in  fact,  significantly  improved  the  technique.  This  method  was  explained 
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in  Sec.  1.3.1  and  thereafter  used  successfully  in  the  design  of  each  mode 
(Chapters  2-4) .  It  is  a  good  example  of  the  importance  of  application  to 
real  problems,  by  the  academic  researcher.  (See  Appendix) 

Also,  "scheduling"  is  a  well-known,  long-standing  very  useful  technique 
in  flight  control  for  economizing  on  the  feedback  loop  bandwidths.  Therefore, 
the  authors  considered  it  important  to  see  how  it  could  be  incorporated  in  a 
simple  manner,  within  the  design  technique.  This  was  treated  in  Sec.  2.8  and 
it  was  seen  that  the  technique  goes  well  with  scheduling,  because  it  reveals 
the  frequency  ranges  in  which  scheduling  is  most  useful,  and  those  in  which 
it  is  of  lesser  usefulness. 

An  area  which  deserves  more  theoretical  treatment  and  practical  experience 
involves  the  use  of  nondiagonal  compensation.  This  was  considered  and  applied 
to  some  extent  in  Sec.  2.9,  but  insufficiently.  This  topic  definitely  requires 
more  research. 

There  is  considerable  interest  now  in  digital  realization  of  compensation 
networks.  A  digital  compensator  always  contributes  effectively  a  right  half- 
plane  zero  in  its  -ransfer  function,  with  its  accompanying  phase  lag,  which 
is  not  good  in  feedback  design.  The  larger  the  sampling  rate,  the  smaller  this 
phase  lag.  The  technique  lends  itself  very  easily  uo  finding  the  smallest 
sampling  rates  commensurate  with  specific  performance  degradations,  thus  per¬ 
mitting  the  designer  to  make  intelligent  trade-offs.  The  minimum  sampling  rate 
needed  for  the  feedback  properties  of  the  system  (parameter  uncertainty,  dis¬ 
turbance  attenuation) ,  need  not  be  the  same  sampling  rate  needed  for  its  filter 
properties,  i.e.  the  command  inputs  nay  not  need  to  be  sampled  anywhere  as 
often  as  the  feedback  sensors.  And  some  of  the  sensors  may  not  need  to  be 
sampled  as  often  as  some  others.  The  feedback  technique  reveals  this  informa¬ 
tion  very  readily. 


Finally,  it  Is  worth  emphasizing  that  the  authors  for  the  most  part 
have  hardly  any  experience  In  flight  control  design.  They  are  specialists 
in  feedback  theory.  Yet  the  synthesis  technique  permitted  them  to  zero  in 
relatively  quickly  and  in  a  systematic  manner,  onto  the  vital  design  conflicts 
a« A  trade-offs.  This  is  the  purpose  of  a  scientific  theory  -  to  enable  a 
non-specialist  in  the  specific  application  to  get  to  the  heart  of  the 
matter,  without  the  need  for  a  long  apprenticeship.  But  it  is  important 
for  the  theorists  to  relatively  early  apply  their  theory  to  realistic 
problems.  In  this  case  we  were  confronted  with  a  situation  which  wreaked 
some  havoc  with  our  theory  (see  p.  13,  and  Appendix),  which  forced  us  to  a 
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IMPROVED  DESIGN  TECHNIQUE  FOR  UNCERTAIN 
MULTIPLE  INPUT-OUTPUT  FEEDBACK  SYSTEMS 

Isaac  Horowitz* 


ABSTRACT 

This  paper  presents  a  synthesis  technique  for  linear  time  invariant 

nxn  multiple  input-output  (mio)  feedback  systems  with  constrained  "plant" 

P,  and  output  feedback.  Due  to  uncertainty,  P  is  known  only  to  be  a  member 

2 

of  a  set  P  ■  {p}.  It  is  required  for  all  Pe  P,  the  n  system  transfer 
functions  t  be  members  of  specified  sets  of  acceptable  outputs. 

A^;  u,  v  *  1,  ...,  n.  The  problem  is  rigorously  converted  into  a  number 
of  single  input-output  (sio)'  uncertainty  problems,  whose  solutions  are 
guaranteed  to  solve  the  original  mio  problem.  The  technique  has  severed, 
advantages  over  a  previous  one,  which  also  converted  the  mio  system  to 
equivalent  sio  systems:  (1)  Fixed  point  theory  is  not  needed  to  rigourously 
justify  the  theory.  The  justification  is  very  simple .  (2)  There  is 

significantly  lass  overdesign  inherent  in  the  method.  (3)  If  arbitrary 
small  sensitivity  is  desired  over  arbitrary  large  bandwidth,  then  the  set  P 
must  satisfy  certain  constraints  u  »*■.  It  is  shown  that  these  constraints 
are  inherent  in  all  linear  time  invariant  compensation  techniques.  In  the 
old  technique,  these  constraints  were  always  present,  even  in  non- arbitrarily 
small  sensitivity  specifications. 
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ABBREVIATIONS 


PS 
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Adj 

adjoint 

a«  X* 

arbitrarily  large 

a*  s  • 

arbitrarily  small 

BM 

bandwidth 

lhp 

left  half  plane 

lti 

linear  time  invariant 

mio 

multiple  input-output 

»P 

minimum-phase 

rhp 

right  half-plane 

sio  single  input-output 


IMPROVED  DESIGN  TECHNIQUE  FOR  UNCERTAIN 
MULTIPLE  INPUT-OUTPUT  FEEDBACK  SYSTEMS 

1.  INTRODUCTION 

A  number  of  synthesis  techniques  have  recently  appeared  for  linear 

tit  invariant  (lti)  n*r>  multiple  input-output  (mio)  feedback  systems 

with  significant  plant  uncertainty  (Horowitz  1979,  Cruz  et  al,  1981, 

Doyle  and  Stein  1981,  Postlethwaite  et  al  1981,  Lehtomaki  et  al  1981, 

Safanov  et  al  1981,  Davison  and  Ferguson  1981,  Zames  and  Bensoussdn  1982) . 

Most  of  these  are  based  on  the  singular  value  approach,  and  do  not  readily 
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permit  detailed  control  over  the  Individual  n  closed- loop  responses,  nor 
trade-off  between  the  loop  bandwidths.  The  method  of  (Horowitz  1979)  does 
permit  such  detailed  control.  It  is  also  quite  transparent  in  revealing  the 
price  paid  for  the  benefits  of  feedback,  in  terms  of  the  various  loop 
transmission  bandwidths,  and  permitting  trade-off  between  them.  It  has 
been  applied  successfully  to  a  number  of  highly  uncertain  mio  problems 
(Horowitz  and  Sidi  1980,  Horowitz  and  Loecher  1981),  including  the  YF16CCV 
3x3  flight  control  system  (Horowitz  et  al  1981) . 

The  method  is  based  on  fixed  point  theory  and  replaces  the  uncertain 
nxn  mio  system  by  a  number  of  single  input-output  (sio)  uncertain  systems, 
whose  solutions  are  guaranteed  to  solve  the  original  mio  problem.  However, 
there  is  some  overdesign  inherent  in  the  method.  This  paper  presents 
several  refinements  of  the  method,  with  reduced  overdesign,  and  retaining 
the  sio  equivalence  approach.  One  of  these  does  not  require  fixed  point 
theory  for  its  derivation,  only  simple  matrix  manipulation.  Also, 
constraints  on  the  mio  plant  are  derived  which  permit  this  design  method 
to  achieve  arbitrarily  small  (a.s.)  sensitivity  over  arbitrarily  large 
(a.l.)  bandwidth.  It  is  shown  that  these  constraints  are  inherent  in  every 
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2.  THE  2x2  SYSTEM 

The  general  nxn  system  Is  defined  in  Fig.  1,  with  the  matrices: 
plant  P  -  ] ,  compensations  G  -  [g  , ]  diagonal  ,  F  all  nxn  matrices  of 

transfer  functions.  Preliminary  quasi-diagonalization  of  the  plant  may 
be  desirable,  which  is  here  absorbed  in  P.  This  is  discussed  in  Sec.  6. 

Hence  diagonal  G  does  not  necessarily  mean  diagonal  compensation  inside 
the  loop  in  Fig.  1.  Due  to  uncertainty,  P  e  P  a  defined  set.  Tolerances 
of  the  form 

i  |tw  |  <  buv(w)  ,  V  P  e  P  (la) 

have  been  assigned  on  the  magnitudes  of  the  elements  of  the  closed  loop 

system  transfer  function  matrix  T  »  [t  ( ju) ] .  Such  tolerances  suffice 

uv 

for  minimum-phase  (mp)  elements  (Horowitz  ]979)  and  determine  sets  of 
acceptable  t^tju)  denoted  by 

Aw(ju)  -  {Aw(ju)>  -  set  of  acceptable  tuv(ju).  (lb) 

Often  one  wants  sene  or  all  nondiagonal  elements  of  T  to  be  small  '/  P  e  P  . 

For  such  elements,  denoted  as  "basically  noninteracting",  auv.  ■  0,  and 

their  f  __  are  also  made  zero  (Horowitz  and  Sidi  1980) . 
uv 

In  the  new  technique,  the  design  equations  (for  the  2x2  mio  system) 

for  t^,  t^2  (or  alternatively  t^,  t^)  are  identical  to  Eqn  (4a, b)  in 

(Horowitz  1979) .  They  are  derived  from  the  matrix  equation  T  ■  (I  +  PG)  1  PGF 

of  Fig.  1,  by  manipulating  it  into  (P-1  +  G)  T  ■  GF,  letting  P  1  *  [P^]  “ 

[1/Q^j]  “  +  vn  ^th  Vd,  the  diagonal  and  nondiagonal  parts  (the 

diagonal  elements  of  V  are  zero  and  vice-versa  for  V.) ,  giving  T  *  (V.  +  G)  1 

n  a  a 

(GF  -  V  T)  which  is  used  to  define  (potential)  fixed-point  mappings, 
n 
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*gufuv  *  dUV^UU  .  mr 

tuv  1  +  gQ  '  uv  Puifciv 

uuu  ij<u 


(2a, b) 

(4a ,b  Horowitz  1979) 


So  the  design  equations  used  here  for  t^  are: 


,  £uLi ;  diiQn  . 

*li  ■  1  +  L!  '  Li  "  9iCii ' 


(3a-d) 


*ii  *  5^  •  Si  ‘  A2i  (lb,,  F'1-  ‘V 


U/Oijl 


The  objective  is  to  choose  f^(s),  f12(s),  3Uch  that  the  4^(3) 

have  no  right  half-plane  (rhp)  poles,  and  satisfy  the  tolerances  (1)  on 
|tii(j“)|,  V  P  e  P  and  V  t2i  e  A  spearing  in  d^.  These  are 
precisely  the  sio  design  problems  of  Fig.  2a. 


When  t^  is  substituted  for  4^  in  (3a) ,  the  resulting  equations  are 


exact.  Hence,  if  indeed  4^  e  A^  for  all  t2i  (of  d^)  e  A2i,  P  e  P  , 


then  the  design  objectives  for  t^  have  been  achieved  by  the  f^  '  f12'  *1 


but  only  if  actually  t2i  e  A^.  A  systematic  design  procedure  for  choosing 


£11'  £12'  gi  for  tile  above  purpose,  has  been  detailed  with  examples  in 
(Horowitz  and  Sidi  1980,  Horowitz  and  Loecher  1981).  A  2x2  nonlinear 
uncertain  plant  exanple  is  given  in  (Horowitz  and  Breiner  1981) . 

The  final  step,  in  the  2x2  system,  is  to  choose  f21,  f22»  g2  to 


ensure  that  the  t2v(s)  have  no  rhp  poles  and  that  t2v  e  A2v,  p  e  P. 


In  (Horowitz  1979)  the  design  equations  are  again  (2a, b)  with  u  *  2  and 


the  d2v  containing  t^.  Instead,  here,  there  are  used  equations  independent 


>-1. 


of  the  t^,  by  simply  finding  t2^  from  T  ■  (I  +  PG)  PGF : 


£2iL2e  *  d2i 


21 


1  +  h 


92222(1  +  L1)  g2f22(1+V 


2e 


2e  ^2*22e  1  "  YL2  +  Li  1  + 


,ij  "  piipjj  ' 


(4a-e) 


91£UP21  (1~  Y12} 


d  m 

2i  1  -  yl2  *  t. 


tPi,]  • 


with  the  sio  structures  of  Fig.  2b.  This  design  is  done  after  Fig.  2a  has 

been  completed  by  mews  of  (3a-d) ,  so  that  L^,  f  f^2  are  known.  Zt  is 

then  necessary  to  find  g 2»  f,  ,  f^2  so  that  in  Fig.  2b  the  outputs 

*^21'  t22  are  8table  satisfy  the  tolerances  (1)  on  jt21|,  ^22^ 

respectively,  which  are  single- loop  problems  similar  to  Fig.  2a,  except  that 

only  the  uncertainty  P  e  P  need  be  considered,  as  the  d^  in  (4a-e)  are  not 

functions  of  the  elements  of  A  ,  which  they  we  in  (3a-d) .  At  each 

uv 

step, design  execution  is  that  of  single  input-output  (sio)  single  loop 
systems — which  is  what  makes  this  design  procedure  so  tractable. 

Susmarizing,  the  theoretical  justification  of  the  above  design  procedure 
is  as  follows:  The  design  specifications  (1)  are  satisfied  for  t2  of 
Fig.  2b,  by  proper  choice  of  the  f 2^,  g2  for  the  given  g  .  So  the 

design  is  satisfactory  for  t2^,  t22  because  (4a ,b)  are  exactly  the  expressions 
for  t21,  t22,  (even  if  the  specs,  for  t^2 ,  t^  are  not  satisfied).  Now 
Figs.  2a  have  been  designed  (via  Bqs.  3a-d)  ,  so  that  are  3table 

and  satisfy  the  specs,  on  t^,  t^2  V  P  e  P  ,  if  the  t^  appearing  in  d^ 
are  in  A2i  (which  they  are).  And  the  equations  for  $^2  correspond 

precisely  to  those  for  t^,  t^2*  Thus,  no  fixed  point  theory  is  needed 
to  rigorously  justify  this  design  procedure,  although  the  idea  and  approach 
were  motivated  by  the  fixed  point  method  in  (Horowitz  1979).  There,  the 
design  equations  for  t2_,  t22  we  of  the  same  form  as  Eqs.  (3a-d)  with  the 
d2^  functions  of  t^,  so  fixed  point  theory  is  required  to  justify  the 
method. 


2.1  Seduced  Overdesi 


Pigs.  2a  are  the  same  as  in  (Horowitz  1979)  in  which  t21 ,  t  appear 

in  the  disturbance  4^,  d^.  There  is  inherent  overdesign  in  Fig.  2a, 

because  in  reality  there  is  correlation  between  the  t^^,  t22  and  the 

2^  of  P  *■.  This  correlation  is  not  being  exploited.  The  uncertainties 

in  t21  e  ^2i'  fc22  e  ^22  ***  assumed  independent  of  P  e  P  in  Figs.  2a 

(see  Horowitz  and  Loecher  1981  for  discussion  how  to  reduce  this  overdesign.) 

But  such  overdesign  does  not  exist  in  Fig.  2b  because  d2^,  4... ...  are  not 

functions  of  the  elements  of  any  A  .  in  (Horowitz  1979)  the  design  for 

uv 

t.^,  t22  involved  overdesign  precisely  as  in  Figs.  2a.  The  above  procedure 
can,  of  course,  be  reversed  with  $21,  $22  using  Eqs.  (2a, b)  with  u  ■  2, 
and  t^,  t12  using  Eqs.  (4a-e)  by  exchanging  numbers  1,2. 

3.  3x3  DESIGN  EQUATIONS 


Nine  design  equations  may  be  defined  in  a  manner  similar  to  (3,4)  of 


Sec.  2.  Let  $lv,  v 


1-3,  be  the  same  as  in  (Horowitz  1979,  Eqs.  4a, b  which 


are  Eqs.  (2a, b)  here) : 

♦ 


flvLl  *  dlv°ll 


rlv 


1  +  t. 


?1211  '  dlv 


*•  p  t 
i»<l  rli^iv 


(5a-c) 


giving  sio  problems.  L^  and  the  three  f^v  are  chosen  so  that  are  stable 
and  satisfy  the  tolerances  (1)  for  all  t^v  e  appearing  in  d^v  amd 
all  P  e  P  -  see  (Horowitz  and  Loecher  1981)  for  procedure  amd  detailed 
design  example.  The  equations  for i^aie  obtained  from  Eq.  2  setting 
u  «  2  in  which  t^v  and  t3v  appear,  but  the  t^v  are  here  replaced  by  the 
of  (5),  giving  for  v  ■  1-3,  in  the  notation  of  Eqs.  (3,4): 


vjl* 


d2v  "  iT  {t3v  [  —2±J±-  ,  p 

92  JV  P11(l+L1)  *23 


]  -  lv  1  21  } 
1  +  1*. 


again  sio  problems.  and  the  f^  cure  known  from  the  designs  of  (5). 
Here,  the  f_  and  L_  are  chosen  so  that  are  stable  for  all  t,  e  A 

2V  2e  2V  JV  jV 

appearing  in  d^v  and,  of  course,  for  all  P  e  P.  Although  the  fonts  for 
d2v  and  I.  in  (6)  are  different  from  those  in  (5) ,  (5a, 6a)  are  otherwise 
identical  in  form,  so  the  design  techniques  for  both  are  basically  the 
same,  as  detailed  in  (Horowitz  and  Loecher  1981) . 

Finally,  sio  design  equations  for  tJv  are  gotten  by  finding  t3v  fro 
T  ■  (I  +  PG)  ^PGF,  or  from  Eq.  (2)  for  u  ■  3  and  eliminating  t^,  t2i 
by  means  of  Eqs.  (5,6): 


(1+1^)  (1+L2)  -  y12 

r23  (1+L1) 

*  T13  (ltL2>  *  *  VS1 

?23P31 

P  P 

u  -  32  21 

P21P33  ' 

3  P31P22 

CliLl233ni 

+  f2iL2Q33n2 

22P21P32  "  P31  1  +  V 


UP12P3 1  -  P32(1  *  V 
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Since  f^v»  f2v  are  known,  the  only  unknowns  in  the  above 

are  the  f^  and  g^.  Eqs.  (7a- j)  constitute  single-loop  uncertinty  problems, 
for  which  the  technique  and  examples  of  the  references  (Horowitz,  etc.) 
apply,  i.e.  they  are  chosen  so  that  the  t^  are  stable  and  satisfy  the 
tolerances  in  (1).  Note,  again,  that  at  each  step,  design  execution  is  that 
of  sio  single- loop  systems. 

The  justification  of  the  above  design  approach  is  the  same  as  of  the 
2x2  case:  Suppose  the  nine  and  three  g^  in  (7a- j)  are  such  that  the 
t^  are  stable  and  their  tolerances  are  satisfied,  which  is  so  by  def initio 
here.  Now  the  design  based  on  (6a-c)  guarantees  that  the  t2^  tolerances 
are  satisfied,  providing  the  t3i  appearing  in  d2i  are  c  A3i  -  which 
is  the  case  here.  Hence,  the  t^  tolerances  are  also  satisfied.  Finally, 
the  design  based  on  (5a-c)  guarantees  that  the  t^  tolerances  are  satisfied, 
providing  that  the  t2i»  t3i,  appearing  in  dly  are  e  A2i,  A3i  respectively, 
which  has  been  established.  In  the  above ,  there  is  some  overdesign  in 
(5a-c)  because  the  t^v  e  A^,  (i  ■  2,3),  (v  -  1,2,3)  appear  as  disturbances 
uncorrelated  to  the  plant  uncertainty.  In  (6a-d)  there  is  less  overdesign 
because  only  the  t3  so  appear,  while  there  is  no  such  overdesign  in 
(7a- j).  Of  course,  the  order  can  be  changed  and  equations  of  the  form 
(Sa-c)  used  for  the  second  or  third  channel,  etc. 

4.  nxn  SYSTEM,  n>3 

The  procedure  for  generating  the  design  equations  for  nxn  mio  systems 
with  n  >  3,  should  be  clear  from  the  above.  One  uses  for  any  channel  (say 
the  first)  design  equations  in  which  all  the  t^  (i  /  1)  appear  as 


disturbances.  Denote  these  as  Eqs.  A.  These  equations  can  be  derived 
or  obtained  from  (Horowitz  1979,  Eqs.  (4a, b)  with  u  ■  1) .  For  the  next 


chosen  channel  (say  the  second),  one  starts  with  Eqs.  (4a, b)  of  (Horowitz  1979) 
with  u  »  2,  but  eliminates  all  the  t^v  by  means  of  Eqs.  A  above.  Denote  the 
resulting  design  equations  as  Eqs.  B.  For  the  next  chosen  channel  (say  the 
third),  start  again  with  Eqs.  (4a, b)  u  »  3,  but  eliminate  all  the  t^y,  t^v 
by  means  of  Eqs.  A,  B  above.  The  process  continues  until  the  end,  and  the 
theoretical  justifi cation  is  the  same  as  given  above  for  n  =  2,3. 

Mixtures  of  the  technique  in  (Horowitz  1979)  and  the  one  described  above, 
may  also  be  used.  For  example,  for  the  3x3  system,  Eq.  2  may  be  used  for  both 
channels  1,2  and  (7a- j)  for  channel  3.  The  theoretical  justification  is 
now  as  follows:  Design  for  channel  3  is  correct  by  definition,  but  fixed 
point  theory,  precisely  as  in  (Horowitz  1979)  is  used  to  justify  the  designs 
for  channels  (1,2) :  The  two  sets  of  design  equations  for  t^v,  t2v  are 
taken  as  the  mappings  on  the  acceptable  sets  and  the  third  set  of 

mappings  is  simply  ^  »  aJv  eA^  The  nine  f^,  three  gu  have  been 
chosen  so  that  these  mappings  map  A^  into  themselves,  etc.,  so  a  fixed  point 
exists,  etc.,  as  in  (Horowitz  1979).  For  larger  n,  it  is  clear  that  a 
larger  variety  of  mixtures  is  possible,  giving  the  designer  useful  flexibility. 
However,  he  must  understand  the  sio  design  theory  used  in  the  design 
execution,  which  reveals  the  cost  of  feedback  and  the  available  tradeoffs 
among  the  loops,  in  order  to  be  able  to  exploit  this  flexibility  to  its 
fullest  extent.  A  detailed  3x3  example  based  on  fixed  point  theory  is  given 
in  (Horowitz  and  Loecher  1981} . 
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5.  CONDITIONS  FOR  EXISTENCE  OF  A  SOLUTION 

This  section  considers  the  conditions  required  for  the  applicability 
of  the  above  design  technique.  Also,  it  considers  the  inherent,  irreducible 
conditions  applicable  for  lti  compensations  in  general,  and  compares  the 
two  sets  of  conditions.  This  is  done  for  "arbitrarily  small  (a.s.)  sensitivity 
defined  as  achievement  of  a.s.  sensitivity  of  the  t  over  arbitrarily  large 
(a.l.)  bandwidth.  Such  "a.s.  sensitivity"  also  achieves  a.l.  attenuation 
over  a.l.  bandwidth,  of  external  disturbances  acting  on  the  plaint.  This 
problem  has  also  been  studied  in  an  abstract  setting  by  (Zames  and  Bensoussan 
1982). 


In  Fig.  2a,  it  is  required  that  e  A  ,  $12  e  A^,  for  all 
P  e  P,  t21  e  A2l»  amd  t  e  A22.  In  the  general  n  case,  the  disturbance 
component  in,  for  exasple,  is 


divQn  y  ^  (Adj  p)u/(Ad3  p)lx 

l+g.Q. .  i/1  fciv  .  .  _ 

1  11  1  +  g  det  P 

”*■  /Vj'i  n  \ 


Right  half-plane  (Rhp)  poles  of  (Adj  P)  ^  are  normally  cancelled  by  similar 
poles  in  det  P,  since  (Adj  P)^  is  a  term  in  expansion  of  det  P.  Rhp 
zeros  of  (Adj  P>11  are,  of  course,  normally  cancelled  by  similar  ones  in 
the  denominator.  There  may  be  exceptional  cases  when  in  det  P,  for  example, 
a  rhp  pole  of  (Adj  P)^  Is  cancelled  by  an  identical  zero  of  p^  and 
does  not  appear  in  the  other  terms  of  det  P.  Such  cases  are  excluded. 


5.1  Conditions  for  "a. s.  sensitivity"  in  single-loop  desie 


In  Eqs.  (3),  it  is  seen  that  a.s.  tolerances  over  a.l.  bandwidth 
(i.e.  "a.s.  sensitivity")  for  t^,  t^2  are  achievable  if  ■  gl®ll  can  be 
made  a.l.  over  a.l.  bandwidth.  Indeed  this  is  at  least  theoretically 
possible,  if  satisfies  certain  constraints.  These  have  been  detailed  in 
(Horowitz  1979,  Appendix  1) ,  so  are  only  qualitatively  described  here  by  means 
of  Fig.  3. 

Fig.  3  is  the  extended  logarithmic  complex  plane  (Nichols  chart) . 

Since  P  ranges  over  P,  the  set  *  g^Q^  ( ju») }  is  not  a  single  complex 
number  (at  any  fixed  u)  in  the  Nichols  chart  but  a  region,  denoted  as 
Tp  (L ( jo> ) ) ,  the  template  of  L^,  which  is  the  same  as  Tp  (Q^(ju)  but 
translated  horizontally  by  20  log  [ ( jm) |  db  and  vertically  by  Arg 
g1(joi)  degrees,  because  -  gl^ll‘  In  a  desi9n  with  significant  plant 
uncertainty  Tp  (L^(  joj) )  must  lie  relatively  high  up,  above  the  zero  db  line 
as  shown  in  Fig.  3  for  u>^.  This  is  so  over  the  important  u  range  of  t^, 
b12  bandwidths  generally) ,  in  order  to  achieve  the  desired  sensitivity 

reduction.  Such  large  values  for  Tp(L^(j<u) )  can  be  maintained  theoretically 
for  any  finite  u  range,  if  is  m.p.  For  those  u>  for  which  Tp  (L^jw) 
is  so  located,  the  uncertainty  in  the  magnitude  and  phase  of  Q^(joi) 

(i.e.  the  area  of  Tp  (Q^tjuj)),  can  be  arbitrarily  large.  (Note,  however, 
that  in  order  to  maintain  Tp  (L^ ( jw) )  above  the  zero  db  line,  any  zeros  of 
Qll(joi)  on  the  ju  axis  must  be  known  and  finite  in  number  in  order  for 
g^(s)  to  be  assigned  poles  there.  (Obviously,  transcendental  compensation 
can  be  used  for  special  countable  cases.)  If  the  range  of  such  ju>  axis 
zeros  is  uncertain,  then  the  specifications (1)  must  be  modified  to  permit 
such  zeros  of  t^.) 


Sooner  or  later  L^(ju)  must  decrease  and  -*■  0  as  u»  -*■  00 .  Stability 

over  the  range  of  P requires  that  Tp  (l^  ( jw) )  move  downward  in  between  the 

vertical  lines  V^,  without  the  points  .  .  .  O^O^O^,...  lying  in  any 

Tp  (L^(jw)).  This  appears  to  allow  (360-26^)  degrees  phase  width  for 

Tp  Q, ,  (jot) ,  with  0  phase  margin.  However  /L,  (jw)  must  be  negative  on  the 
li  n  —  l 

average,  in  order  for  1 1,^ ( jti»)  |  to  decrease.  So  in  practice  only  (180-  0^) 
phase  width  is  tolerable  in  this  range.  As  u>  increases  and  the  Tp  (L^(jcu)) 
descend  lower  on  the  chart  below  the  zero  db  line,  clearly  their  width  may 
increase  again,  but  it  is  essential  that  the  points  ...0^,0^,...  never  be  a 
part  of  any  Tp  (L^jw)).  unstable  are  included  in  the  above  discussion 

and  don't  require  separate  treatment.  It  follows  from  the  above  that  the 
manageable  uncertainties  depend  on  the  assigned  t^v  tolerances,  but  two 
important  constraints  are  stated  here  for  the  case  of  "a.s.  sensitivity." 

(5.1a)  *  det  P/Adj.  P)  ^  must  be  mp. 

(5.1b)  Sqapose^w.,  }  for  all  w  >  w.  ,  the  width  of  Tp  (q  (ju>))  exceeds 
n  n  11 

^180-0^),  9^  a  desired  phase  margin,  then  it  is  impossible  to  achieve  a.s. 

tolerances  for  u»  >  w.  .  This  prevents  "a.s.  sensitivity"  if  Q  -  kir(s+z.)/ 

h  uu  l 

it  (s+  P.) ,  with  the  k  uncertainty  including  a  sign  change  which  is  independent 
of  the  signs  of  the  z^,  p^.  Also  excluded  is  a  factor  (1+Ts)  in  the  ntxnerator 
or  denominator  of  Quur  with  the  uncertainty  in  f  including  a  sign  change  which 
is  independent  of  other  parameters. 

5.2  Application  of  5.1  to  New  Design  Technique 

The  above  constraints  therefore  apply  to  the  Qyu  *  det  P/(Adj  P)QU 


of  the  first  channel  u,  used  in  the  new  design  technique  of  Secs.  2,  3.  So 


from  (5.1a),  0^  must  be  m.p.  Vp  e  P  (rhp  poles  are  tolerable),  for  "a.s. 
sensitivity"  design.  Suppose  P  ■  tpuv)  has  each  p^v  -*■  kuv/se  as  s  +  ®  . 
For  n  ■  2 


P11P22  ~  >12*21  ^  *11*22  "  *12*21 


uu 


'  w 


k  se 
w 


(8) 


Let  K  »  [k^] ,  so  (5.1b)  states  that  there  may  be  no  change  in  the  sign  of 


det  K A  &*  P  ranges  over  P.  In  this  paper  it  is  assumed  that  all  k^v  >  0 


for  all  P  e  P,  so  (5.1b)  gives  k^jk  >  *12*21  *,p  e  P'  or  wice  versa*  To 


remove  the  ambiguity  it  is  also  assumed  that  the  plant  terminals  are  numbered 


so  that  for  at  least  one  P  e  P,  ^^22  >  *12*21'  30  giv®3  the  constraint 


*11*22  >  *12*21  ^  P  e  P 


(9) 


This  is  a  diagonal  dominance  condition  as  3  +  ».  (ZameS  and  Bensoussan) 
have  defined  a  diagonal  dominance  condition  as  s  ■+■  <*,  in  more  abstract  form. 

The  above  applies  to  the  first  channel,  say  no.  1,  for  which  (3)  is 
used.  Eqs.  (4)  are  used  for  the  second  channel.  The  mp  condition  of 
(5.1a)  therefore  applies  to  Q^d+L^),  no8t  o£  whi°h  is  not  new  because 
mp  (1+L^)»  det  P  are  already  required.  As  for  (5.1b),  there  are  two  extremes. 
"A.s.  sensitivity"  can  be  achieved  by  L^  bandwidth  >>  that  of  L^  (denoted 
by  BW (L^)  >>  BH(L^) ) ,  and  then  (4b)  implies  (5.1b)  applies  to  Q22-  It  can 
also  be  achieved  by  the  opposite  strong  inequality,  and  then  (5.1b)  applies 
to  P22'  For  the  condition  assumed  with  Eqs.  8,9)  of  no  sign  changes  in  the 
k^,  the  results  are  the  same.  It  may  also  be  so  in  the  general  case  but 
this  would  require  consideration  of  simultaneous  sign  changes  among  the  kuv. 


which  is  not  done  here. 


n-3.  Application  of  (5.1)  to  the  first  channel  makes  it  applicable  to 

Q  .  Application  to  the  second  (eqs.  6)  gives  the  same  results  as  to  the 

second  in  n«2  above ,  because  L^e  has  the  same  form  in  both  (compare  4b,  6b). 

If  "a.s.  sensitivity"  is  achieved  (as  it  may  be)  by  BW (L^ )  >>  BW(L2)  >>  BW(L3) 

the  result  is  that  (5.1)  apply  to  Q^,  ®22'  ®33*  If  tiie  °PP°s^te  is  d0116* 

0 

in  the  notation  of  (3.7),  they  apply  to  22  , 

1  ^12 

_ Q33(1  ***^12) _ 

1  ‘  <Y12  +  Y23  +  Y13}  +  <Y12U2  +  '(13V3) 

It  has  not  been  ascertained  whether  these  two  sets  of  constraints  are 
identical.  However  Sec.  5.3  shows  that  (5.1)  must  always  apply  to  each  Quu* 
u  -  1  to  n.  Constrains  for  n  >  3  may  be  similarly  developed. 

5. 3  Inherent  Constrains 

It  is  important  to  determine  whether  the  above  constrains  are  due  to 

the  specific  design  technique,  or  are  inherent  in  the  problem  itself.  For 

this  purpose  examine  (4a)  for  t22.  How  can  "a.s.  sensitivity"  of  t^2  be 

achieved  despite  large  uncertainty  in  P?  Clearly  by  large  L2e,  the  usual 

feedback  method.  Large  L2e  is  achieved  by  large  92^22'  because  lar<?e 

(needed  likewise  for  small  t, ,  sensitivity)  gives  L.  *  g_Q,,_  *  L_. 

lx  «e  «  44  z 

The  latter  also  attenuates  d22,  which  may  not  be  small  because  of  g^  in  is 
numerator.  This  same  principle  applies  to  all  t  ,  and  is  basically  the 
same  as  that  derived  from  examination  of  (3),  i.e.  there  is  need  for  a. 1. 

L2 ,  over  a. 1.  bandwidth  in  order  to  achieve  "a.s.  sensitivity."  But 
do  the  constraints  of  Sec.  5.1,  in  particular  (5.1a,b)  apply  to  Q22  and  Q^? 
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This  is  indeed  so,  and  proven  by  (3a)  and  its  analog  for  t^  (by 
interchanging  1,2) ,  by  simply  asking  whether  stable  t  is  possible  if  1  + 
has  rhp  zeros?  For  if  not,  and  since  a.l.  over  a.l.  bandwidth  is  needed, 
it  follows  that  (and  L^)  must  satisfy  constraints  (5.1).  Suppose  (1  +  L^) 
has  rhp  zeros.  These  are  rhp  poles  of  t^,  unless  in  (3a  for  i  »  1)  the 
numerator  of  t^  has  these  same  zeros.  Suppose  it  has  them,  and  there  is  a 
small  change  in  f^.  Since  F  is  outside  the  feedback  loops,  system  stability 
is  unaffected.  The  zeros  of  1  +  g^Q^  in  (3a)  are  thereby  unaffected,  so 
neither  should  the  zeros  of  the  numerator  of  (3a  for  i  «  1) .  Hie  term  f^g^Q^ 
is  unaffected,  but  t21  is  affected— see  (4)  with  i  ■  1.  Hence,  the  hypothesis 
(1  +  L^)  has  rhp  zeros  is  untenable,  so  (5.1)  applies  to  (and  Q22) . 

Comparing  these  results  with  Sec.  5.2,  the  conclusion  is  that  "a.s. 
sensitivity"  may  be  achieved  by  the  new  design  technique,  with  plant  constraints 
which  are  inherent  and  irreducible,  i.e.  not  more  severe  than  inherently 
necessary.  This  is  achieved  by  letting  BHlL^)  >>  BW(L2>  >>  BW(L3),  wherein 
the  only  constraints  are  (5.1)  applicable  to  Q^,  Q22#  Q33»  which  have  been 
shown  to  be  inherent.  This  is  associated  with  the  following  design  order: 
first  (Eqs.  5),  second  L2 (6) ,  third  L3(7).  The  design  procedure  is 

facilitated  by  such  inequalities ,  because  then  in  (6b)  ,  I>2e  -*■  I<2  -  g2Q22 
over  its  important  design  range,  and  ■*  «  ^3^33'  If  other  factors  are 

equal  (Tp  (Q  ) ,  tolerances  on  t^,  i  ■  1-3  for  each  u) ,  there  is  a  natural 
tendency  for  this  order  of  the  inequalities  because  of  the  inherent  greater 
overdesign  of  L^,  lesser  of  I>2  and  least  of  L3  (recall  Secs.  2.1,  and 
last  paragraph  in  Sec.  3. ) . 
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However,  the  above  is  rather  of  theoretical,  academic  interest,  because 

it  applies  only  for  "a.s.  sensitivity,"  defined  at  the  beginning  of  Sec.  5. 

Thus  it  applies  if  one  is  given  a  plant  with  specific  uncertainty  range  and 

is  challenged  to  achieve  a.s.  performance  tolerances  over  a.l.  bandwidths. 

However,  given  a  mio  plant  set  P  and  sets  of  tolerances  A^v  (u,v  »  n) ,  it  is 

conceivable  that  the  latter  are  not  achievable  by  the  new  technique,  but  are 

inherently  achievable.  The  reason  (e.g.  for  n  -  2)  is  that  in  Eq.  (2),  the 

demand  on  to  achieve  the  desired  sensitivity  reduction,  may  be  greater  than 

is  inherently  needed,  because  of  the  overdesign,  discussed  in  Sec.  2.1,  and  is 

therefore  unachievable,  because  of  the  nature  of  P,  say  it  has  some  non-rap 

elements.  Thus,  a  specific  problem  may  be  inconpatible  with  (5.1)  only  due 

to  this  overdesign.  There  may  exist  as  yet  undiscovered  better  methods  with 

reduced  demand  on  L. ,  which  renders  them  conq>atible  with  (5.1).  From  Sec.  2.1, 

it  is  clear  that  the  best  method  will  be  achieved  with  Eqs.  (3)  by  maximum 

use  of  the  correlation  which  exists  between  the  P  uncertainty  and  the  t^  e 

in  (3a-d) .  A  suggestion  for  this  purpose  has  been  given  in  (Horowitz  and 

Loecher  1981) .  One  subdivides  the  plant  set  into  subsets  P^  which  are 

correlated  with  subsets  A^^  of  A^.  Eqs.  (3)  (similarly  5,6)  are  now  applied 

to  these  pairs  P. ,  A  separately  for  each  i.  This  approach  has  not  as  yet 
i  uvi 

been  attempted  in  any  numerical  problem. 

It  is  worth  noting  that  constraints  of  the  diagonal  dominance  type  as 
appear  in  the  design  technique  of  (Horowitz  1979) .  However,  they  are 
always  present  there,  even  if  "a.s.  sensitivity"  is  not  attempted.  In  the 
new  technique  they  are  in  effect  only  for  "a.s.  sensitivity."  Hence,  it  is 
possible  that  a  specific  synthesis  problem  with  given  A^v,  P  sets  may  not  be 
solvable  by  the  older  method,  but  is  solvable  by  the  new  method.  This  was  the 
case  in  (Horowitz  et  al  1981) . 


6.  NOKD1AGONAL  G 


The  constraints  on  P  in  Sec.  5  were  deduced  on  the  assumption  G  is  diagonal 
Are  these  constraints  eased  if  nondiagonal  G  is  used?  To  answer  this  question, 
let  H  be  a  fixed  lti  pre-coopensator  matrix  inserted  ahead  of  the  plant  and  let 

V  «  ph  be  the  new  effective  plant  in  the  set  V  »  {PH,  P  e  P} .  The  design 
techniques  with  diagonal  G,  are  now  applied  to  set  V  instead  of  set  P.  If 
H  is  helpful  in  overcoming  some  constraint,  then  it  is  necessary  that  the 
constraint  violated  by  P,  is  not  violated  by  I/.  The  constraint  det  P  is  mp. 

V  P  e  P  is  not  eased  at  all,  because  det  v  -  (det  P)  (det  H) ,  and  obviously 
cancellation  of  rhp  zeros  of  det  P  by  det  H  cannot  be  done  for  many  reasons. 

The  other  important  constraint  involving  diagonal  dominance  as  s  -*■  ®,  is  also 
not  eased,  because  it  applies  to  the  sign  of  det  V  not  changing,  u  s  +  •  . 

Thus  the  constraints  on  P  for  "a.s.  sensitivity"  are  not  eased  by  nondiagonal 
G. 

However,  H  may  be  very  helpful  in  reducing  the  amount  of  feedback  needed 
to  achieve  specified  tolerance  sets  (1) ,  for  a  given  plant  set  P,  so  that  a 
design  unachievable  by  diagonal  G  (say,  because  of  non-mp  P  or  sensor  noise 
problems)  may  be  achievable  via  H.  For  example  in  Eq.  (4a) ,  L2q  must  handle 
the  uncertainties  due  to  itself  and  attenuate  the  effective  distub an ce 
set  {d2i>.  For  basically  noninteracting  tolerances  on  tuv(ujrfv),  is  made 
zero,  so  only  the  latter  need  exist.  It  may  be  possible  to  considerably 
reduce  I  d2i  Lax  by  088118  of  H,  by  making  V  -  PH  quasidiagonal,  even  though 
P  has  large  nondiagonal  components. 

Off  diagonal  plant  elements  appear  in  all  the  design  equations  (2-7)  in 
the  'disturbance'  coofKxients,  so  their  reduction  via  H  is  desirable.  How  is 
this  systematically  done  in  the  case  of  significant  P  uncertainty?  For  n  *  2, 


let  the  normalized  H  have  1  for  its  diagonal  elements  and  h 


*  v. 


Then  v12  •  pp^  +  p12,  v21  -  p21  +  VP22‘  The  ob3ective  is  minimize  over  P 

Pl2  P22 

maxlv,  J  ,  I v_ _ |  at  each  U  .  Sketches  of  the  sets  { -  (jw)},  { -  (ju)} 

U  21  P11  P21 
in  the  cosplex  plane,  are  clearly  very  helpful  in  choosing  p(ju>),  v(ju>). 


However,  one  should  check  the  effect  on  the  resulting  sets  of  v^  *  p^  +  vp^ 
V22  *  ^p21  +  p22'  h®08118®  of  the  obligations  on  the  loop  transmissions  due 
to  their  uncertainties .  The  fined,  choice  depends  on  the  relative  isqportance 


of  the  two  terms  in  the  numerators  of  ♦  (See  Horowitz  and 

Loecher  1981  Secs.  3.2,4  for  discussion  relevant  to  this  topic). 


If  the  elements  of  P  have  a  rhp  pole  in  common,  i.e.  P  *  P^/(s-p), 
then  one  should  not  try  to  diagonalize  P  by  means  of  PH  *  A  diagonal,  because 
in  practice  H  *  P  ^  A  with  P  1  j4  P  1  exactly,  giving  PH  «  P  P  1  A  with  rhp 

&  &  3L 

dipoles.  Instead,  one  tries  to  diagonalize  P1  by  means  of  P^H  -  A  ,  giving 
H  -  P”J  A,  and  PH  -  PjP”2  A/(s  -  p) . 


7.  CONCLUSIONS 

This  paper  has  presented  new  synthesis  techniques  for  highly  uncertain 
nxn  mio  lti  feedback  systems  with  output  feedback,  with  the  following 
features: 

2 

(a)  There  is  detailed  control  over  the  n  individual  system  transfer 
functions . 

(b)  The  mio  uncertainty  problem  is  rigorously  converted  into  a  number 
of  sio  uncertainty  problems.  Solutions  of  the  latter  are 
guaranteed  to  be  satisfactory  for  the  former.  Relatively 
siaqple  sio, single  loop  feedback  techniques  can  be  used  to  solve 
the  sio  problems. 


(c)  For  "arbitrary  small  sensitivity"  over  arbitrary  large  bandwidth 
the  technique  in  Secs.  2,  3  give  constraints  on  the  plant  which 
cure  inherent  and  irreducible,  i.e.  every  lti  compensation 
technique  has  these  constraints. 

(d)  Part  of  the  constraints  (at  infinite  s)  in  (c)  were  always 
present  in  the  previously  developed  sio  equivalence  techniques 
of  (Horowitz  1979),  i.e.  even  if  "a.s.  sensitivity"  was  not 
required.  They  are  present  in  the  new  techniques  only  for  a.s. 
sensitivity.  Also,  fixed  point  theory  is  not  required  for 
justification  of  the  new  technique. 

(e)  The  overdesign  inherent  in  the  fixed  point  techniques  of 
(Horowitz  1979) ,  has  been  reduced  but  some  overdesign  is  still 


present 
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